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The last decade has witnessed both quantitative and quaUtative progresses in Shell Model stud- 
ies, which have resulted in remarkable gains in our understanding of the structure of the nucleus. 
Indeed, it is now possible to diagonalize matrices in determinantal spaces of dimensionality up 
to 10^ using the Lanczos tridiagonal construction, whose formal and numerical aspects wc will 
analyze. Besides, many new approximation methods have been developed in order to overcome 
the dimensionality limitations. Furthermore, new effective nucleon-nucleon interactions have been 
constructed that contain both two and three-body contributions. The former are derived from 
realistic potentials (i.e., consistent with two nucleon data). The latter incorporate the pure 
monopole terms necessary to correct the bad saturation and shell-formation properties of the real- 
istic two-body forces. This combination appears to solve a number of hitherto puzzling problems. 
In the present review we will concentrate on those results which illustrate the global features of 
the approach: the universality of the effective interaction and the capacity of the Shell Model 
to describe simultaneously all the manifestations of the nuclear dynamics either of single particle 
or collective nature. We will also treat in some detail the problems associated with rotational 
motion, the origin of quenching of the Gamow Teller transitions, the double /3-decays, the effect of 
isospin non conserving nuclear forces, and the specificities of the very neutron rich nuclei. Many 
other calculations — that appear to have "merely" spectroscopic interest — are touched upon briefly, 
although we are fully aware that much of the credibility of the Shell Model rests on them. 



Contents 



IV. The Lanczos basis 



II. 



III. 









A. 


Level densities 


Introduction 


2 




B. 


The ground state 


A. The three pillars of the shell model 


2 






1. The exp(S) method 


B. The competing views of nuclear structure 


3 






2. Other numerical approximation methods 


1. Collective vs. Microscopic 


3 






3. Monte Carlo methods 


2. Mean field vs Diagonalizations 


4 




C. 


Lanczos Strength Functions 


3. Realistic vs Phenomenological 


4 








C. The valence space 


5 


V. 


The Ohu} calculations 


D. About this review: the unified view 


6 




A. 


The monopole problem and the three-body 
interaction 


The interaction 


7 




B. 


The pf shell 


A. Effective interactions 


8 






1. Spectroscopic factors 


1. Theory and calculations 


8 






2. Isospin non conserving forces 


B. The monopole Hamiltonian 


9 






3. Pure spectroscopy 


1. Bulk properties. Factorable forms 


10 




C. 


Gamow Teller and magnetic dipole strength. 


2. Shell formation 


11 






1. The meaning of the valence space 


C. The multipole Hamiltonian 


13 






2. Quenching 


1. Collapse avoided 


14 








D. Universality of the realistic interactions 


15 


VI. 


Spherical shell model description of nuclear 






rotations 


The solution of the secular problem in a finite 






A. 


Rotors in the pf shell 


space 


16 




B. 


Quasi-SU(3). 


A. The Lanczos method 


16 




C. 


Heavier nuclei: quasi-|-pseudo SU(3) 


B. The choice of the basis 


17 




D. 


The ^^Ar and '^'^Ca super-deformed bands 


C. The Glasgow m-scheme code 


18 




E. 


Rotational bands of unnatural parity 


D. The m-scheme code ANTOINE 


18 








E. The coupled code NATHAN. 


19 


VII. 


Description of very neutron rich nuclei 


F. No core shell model 


20 




A. 


N=8; "Li: halos 


G. Present possibilities 


20 




B. 


N=20; "^^Mg, deformed intruders 



c. 

D. 



:28: Vulnerability 
:40; from "magic" 



^Ni to deformed ''^Gr 



21 
21 
22 
23 

24 
24 
25 

26 

26 
29 
29 

30 

30 
32 
32 
33 



34 
35 
36 
38 
39 
40 

42 
42 
43 

44 

45 



'Electronic address: etienne . caurierOires . in2p3 . f r 
^Electronic address: martinezOieec. f cr . es 
tElectronic address: frederic.nowacki@ires.in2p3.fr 
^Electronic address: alfredo.poves9uajii.es 
^Electronic address: andres .zukerOires . in2p3 . f r 



VIII. Other regions and themes 47 

A. Astrophysical applications 47 

B. /3/3-decays 50 
G. Radii isotope shifts in the Calcium isotopes 51 

D. Shell model calculations in heavier nuclei 52 

E. Random Hamiltonians 53 



2 



IX. Conclusion 53 

Acknowledgments 54 

A. Basic definitions and results 54 

B. Full form of Tim 55 

1. Separation of Hmnp and HmO 56 

2. Diagonal forms of Tim 57 

C. The center of mass problem 57 
References 58 



I. INTRODUCTION 

In the early days of nuclear physics, the nucleus, com- 
posed of strongly interacting neutrons and protons con- 
fined in a very small volume, didn't appear as a system to 
which the shell model, so successful in the atoms, could 
be of much relevance. Other descriptions — based on the 
analogy with a charged liquid drop — seemed more natu- 
ral. However, experimental evidence of independent par- 
ticle behavior in nuclei soon began to accumulate, such 
as the extra binding related to some precise values of the 
number of neutrons and protons (magic numbers) and 
the systematics of spins and parities. 

The existence of shell structure in nuclei had already 
been noticed in the thirties but it took more than a 
decade and numerous papers (see Elliott and Lane, 1957, 
for the early history) before the correct prescription was 
found by Mayer (1949) and Axel, Jensen, and Suess 
(1949). To explain the regularities of the nuclear proper- 
ties associated to "magic numbers" — i.e., specific values 
of the number of protons Z and neutrons the authors 
proposed a model of independent nucleons confined by a 
surface corrected, isotropic harmonic oscillator, plus a 
strong attractive spin-orbit term^, 

U{r) = ^hior^ +dP + Cl- s. (1) 

In modern language this proposal amounts to assume 
that the main effect of the two-body nucleon-nucleon 
interactions is to generate a spherical mean field. The 
wave function of the ground state of a given nucleus is 
then the product of one Slater determinant for the pro- 
tons and another for the neutrons, obtained by filling 
the lowest subshells (or "orbits") of the potential. This 
primordial shell model is nowadays called "independent 
particle model" (IPM) or, "naive" shell model. Its foun- 
dation was pioneered by Brueckner (1954) who showed 
how the short range nucleon-nucleon repulsion combined 
with the Pauli principle could lead to nearly independent 
particle motion. 



^ We assume throughout adimensional oscillator coordinates, i.e., 

r — > (mw/hy/'^ r. 



As the number of protons and neutrons depart from 
the magic numbers it becomes indispensable to include in 
some way the "residual" two-body interaction, to break 
the degeneracies inherent to the filling of orbits with 
two or more nucleons. At this point difficulties accu- 
mulate: "Jensen himself never lost his skeptical attitude 
towards the extension of the single-particle model to in- 
clude the dynamics of several nucleons outside closed 
shells in terms of a residual interaction" (Weidenmiiller, 
1990). Nonetheless, some physicists chose to persist. One 
of our purposes is to explain and illustrate why it was 
worth persisting. The kcypoint is that passage from the 
IPM to the interacting shell model, the shell model (SM) 
for short, is conceptually simple but diflicult in practice. 
In the next section we succinctly review the steps in- 
volved. Our aim is to give the reader an overall view of 
the Shell Model as a sub-discipline of the Many Body 
problem. "Inverted commas" are used for terms of the 
nuclear jargon when they appear for the first time. The 
rest of the introduction will sketch the competing views 
of nuclear structure and establish their connections with 
the shell model. Throughout, the reader will be directed 
to the sections of this review where specific topics are 
discussed. 

A. The three pillars of the shell model 

The strict validity of the IPM may be limited to closed 
shells (and single particle -or hole- states built on them), 

but it provides a framework with two important compo- 
nents that help in dealing with more complex situations. 
One is of mathematical nature: the oscillator orbits de- 
fine a basis of Slater determinants, the m-scheme, in 
which to formulate the Schrodinger problem in the oc- 
cupation number representation (Fock space). It is im- 
portant to realize that the oscillator basis is relevant, 
not so much because it provides an approximation to the 
individual nucleon wave-functions, but because it pro- 
vides the natural quantization condition for self-bound 
systems. 

Then, the many-body problem becomes one of diago- 
nalizing a simple matrix. We have a set of determinantal 
states for A particles, ci|i • • • aj^|0) = \<pi), and a Hamil- 
tonian containing kinetic /C and potential energies V 

U = ^lCija\aj- ^ Vijkia\a\akai (2) 

ij i<j k<l 

that adds one or two particles in orbits i, j and removes 
one or two from orbits k, I, subject to the Pauli prin- 
ciple ({ajoj} = Sij). The eigensolutions of the prob- 
lem |$o,) = J2i C7,Q I (/>/), are the result of diagonalizing 
the matrix {(j)j\'H\(l)i') whose off-diagonal elements are 
either or ±Vijki ■ However, the dimensionalities of the 
matrices — though not infinite, since a cutoff is inherent 
to a non-relativistic approach — are so large as to make 
the problem intractable, except for the lightest nuclei. 
Stated in these terms, the nuclear many body problem 
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docs not differ much from other many fcrmion problems 
in condensed matter, quantum liquids or cluster physics; 
with which it shares quite often concepts and techniques. 
The differences come from the interactions, which in the 
nuclear case are particularly complicated, but paradox- 
ically quite weak, in the sense that they produce suf- 
ficiently little mixing of basic states so as to make the 
zeroth order approximations bear already a significant 
resemblance with reality. 

This is the second far-reaching physical component 
of the IPM: the basis can be taken to be small enough 
as to be, often, tractable. Here some elementary def- 
initions are needed: The (neutron and proton) ma- 
jor oscillator shells of principal quantum number p = 
0, 1, 2, 3 . . ., called s, p, sd, pf ■ ■ ■ respectively, of energy 
!iu){p + 3/2), contain orbits of total angular momentum 
j = 1/2, 3/2 . . .p + 1/2, each with its possible jz pro- 
jections, for a total degeneracy Dj = 2j + 1 for each 
subshell, and Dp = {p + l){p + 2) for the major shells 
(One should not confuse the p = 1 shell — ^p-shell in the 
old spectroscopic notation — with the generic harmonic 
oscillator shell of energy fiuj{p + 3/2)). 

When for a given nucleus the particles are restricted 
to have the lowest possible values of p compatible with 
the Pauli principle, we speak of a Ofio) space. When the 
many-body states are allowed to have components in- 
volving basis states with up to N oscillator quanta more 
than those pertaining to the Ohuj space we speak of a 
Nhu) space (often refered as "no core") 

For nuclei up to A « 60, the major oscillator shells pro- 
vide physically meaningful Ofko spaces. The simplest pos- 
sible example will suggest why this is so. Start from the 
first of the magic numbers N = Z = 2; ^He is no doubt a 
closed shell. Adding one particle should produce a pair of 
single-particle levels: and Pi/2- Indeed, this is what 
is found in ^He and ^Li. What about ^Li? The lowest 
configuration, should produce four states of angular 
momentum and isospin JT = 01, 10, 21, 30. They are ex- 
perimentaly observed. There is also a JT = 20 level that 
requires the P3/2P1/2 configuration. Hence the idea of 
choosing as basis p"^ [p stands generically for p^/ 2 Pi 12) 
for nuclei up to the next closure N = Z = 8, i.e., ^^O. 
Obviously, the general Hamiltonian in (2) must be trans- 
formed into an "effective" one adapted to the restricted 
basis (the "valence space"). When this is done, the re- 
sults are very satisfactory (Cohen and Kurath, 1965). 
The argument extends to the sd and pf shells. 

By now, we have identified two of the three "pillars" 
of the Shell Model: a good valence space and an effec- 
tive interaction adapted to it. The third is a shell model 
code capable of coping with the secular problem. In Sec- 
tion I.C we shall examine the reasons for the success of 
the classical Ofko SM spaces and propose extensions ca- 
pable of dealing with more general cases. The interaction 
will touched upon in Section LB. 3 and discussed at length 
in Section II. The codes will be the subject of Section HI. 



B. The competing views of nuclear structure 

Because Shell Model work is so computer-intensive, it 

is instructive to compare its history and recent devel- 
opments with the competing — or alternative — views of 
nuclear structure that demand less (or no) computing 
power. 



1. Collective vs. Microscopic 

The early SM was hard to reconcile with the idea of 
the compound nucleus and the success of the liquid drop 
model. With the discovery of rotational motion (Bohr, 
1952; Bohr and Mottelson, 1953) which was at first as 
surprising as the IPM, the reconciliation became ap- 
parently harder. It came with the realization that 
collective rotors are associated with "intrinsic states" 
very well approximated by deformed mean field determi- 
nants (Nilsson, 1955), from which the exact eigenstates 
can be extracted by projection to good angular momen- 
tum (Peierls and Yoccoz, 1957); an early and spectacular 
example of spontaneous symmetry breaking. By further 
introducing nuclear superfluidity (Bohr et al., 1958), the 
unified model was born, a basic paradigm that remains 
valid. 

Compared with the impressive architecture of the uni- 
fied model, what the SM could offer were some striking 
but isolated examples which pointed to the soundness of 
a many-particle description. Among them, let us men- 
tion the elegant work of Talmi and Unna (1960), the 
fy2 model of McCullen, Bayman, and Zamick (1964) — 
probably the first successful diagonalization involving 
both neutrons and protons — and the Cohen and Kurath 
(1965) fit to the p shell, the first of the classical OfkJ re- 
gions. It is worth noting that this calculation involved 
spaces of m-scheme dimensionalities, dm, of the order of 
100, while at fixed total angular momentum and isospin 
djT ~ 10. 

The microscopic origin of rotational motion was found 
by Elliott (1958a,b). The interest of this contribution 
was immediately recognized but it took quite a few years 
to realize that Elliott's quadrupolc force and the under- 
lying SU(3) symmetry were the foundation, rather than 
an example, of rotational motion (see Section VI). 

It is fair to say that for almost twenty years after its 
inception, in the mind of many physicists, the shell model 
still suffered from an implicit separation of roles, which 
assigned it the task of accurately describing a few spe- 
cially important nuclei, while the overall coverage of the 
nuclear chart was understood to be the domain of the 
unified model. 

A somehow intermediate path was opened by the well 
known Interacting Boson Model of Arima and lachello 
(1975) and its developments, that we do not touch here. 
The interested reader can find the details in (lachello and 
Arima, 1987). 
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2. Mean field vs Diagonalizations 

The first realistic^ matrix elements of Kuo and Brown 
(1966) and the first modern shell model code by French. 
Halbert, McGrory, and Wong (1969) came almost simul- 
taneously, and opened the way for the first generation 
of "large scale" calculations, which at the time meant 
djT ~ 100-600. They made it possible to describe the 
neighborhood of ^^^O (Zuker, Buck, and McGrory, 1968) 
and the lower part of the sd shell (Halbert et ai, 1971). 
However, the increases in tractable dimensionalities were 
insufficient to promote the shell model to the status of 
a general description, and the role-separation mentioned 
above persisted. Moreover the work done exhibited seri- 
ous problems which could be traced to the realistic ma- 
trix elements themselves (sec next Section LB. 3). 

However, a fundamental idea emerged at the time: 
the existence of an underlying universal two-body in- 
teraction, which could allow a replacement of the uni- 
fied model description by a fully microscopic one, based 
on mean-field theory. The first breakthrough came 
when Baranger and Kumar (1968) proposed a new form 
of the unified model by showing that Elliott's quadrupole 
force could be somehow derived from the Kuo-Brown ma- 
trix elements. Adding a pairing interaction and a spher- 
ical mean field they proceeded to perform Hartree Fock 
Bogoliubov calculations in the first of a successful se- 
ries of papers (Kumar and Baranger, 1968). Their work 
could be described as shell model by other means, as it 
was restricted to valence spaces of two contiguous major 
shells. 

Tliis limitation was overcome when Vauthcrin and 
Brink (1972) and Decharge and Gogny (1980) initi- 
ated the two families of Hartree Fock (HF) calculations 
(Skyrme and Gogny for short) that remain to this day 
the only tools capable of giving microscopic descriptions 
throughout the periodic table. They were later joined by 
the relativistic HF approach (Serot and Walecka, 1986). 
(For a review of the three variants sec Bender, Hce- 
nen, and Rcinhard, 2003b). See also Peru, Girod, and 
Berger (2000) and Rodriguez-Guzman, Egido, and Rob- 
ledo (2000) and references therein for recent work with 
the Gogny force, the one for which closest contact with 
realistic interactions and the shell model can be estab- 
lished (Sections VI. A and II. D). 

Since single determinantal states can hardly be ex- 
pected accurately to describe many body solutions, ev- 
erybody admits the need of going beyond the mean field. 
Nevertheless, as it provides such a good approximation 
to the wavefunctions (typically about 50%) it suggests 
efficient truncation schemes, as will be explained in Sec- 
tions IV.B.l and IV.B.2. 



3. Realistic vs Phenomenological 

Nowadays, many regions remain outside the direct 
roach of the Shell Model, but enough has happened in 
the last decade to transform it into a unified view. Many 
steps — outlined in the next Section I.D — are needed to 
substantiate this claim. Here we introduce the first: A 
unified view requires a unique interaction. Its free pa- 
rameters must be few and well defined, so as to make 
the calculations independent of the quantities they are 
meant to explain. Here, because we touch upon a dif- 
ferent competing view, of special interest for shell model 
experts, we have to recall the remarks at the end of the 
first paragraph of the preceding Section I.B.2: 

The exciting prospect that realistic matrix elements 

could lead to parameter-free spectroscopy did not mate- 
rialize. As the growing sophistication of numerical meth- 
ods allowed to treat an increasing number of particles 
in the sd shell, the results became disastrous. Then, two 
schools of thought on the status of phenomenological cor- 
rections emerged: In one of them, all matrix elements 
were considered to be free parameters. In the other, 
only average matrix elements ( "centroids" ) — related to 
the bad saturation and shell formation properties of the 
two-body potentials — needed to be fitted (the monopole 
way). The former lead to the "Universal srf" (USD) in- 
teraction (Wildcnthal, 1984), that enjoyed an immense 
success and for ten years set the standard for shell model 
calculations in a large valence space (Brown and Wildcn- 
thal, 1988), (see Brown, 2001, for a recent review). The 
second, which we adopt here, was initiated in Eduardo 
Pasquini PhD thesis (1976), where the first calculations 
in the full pf shell, involving both neutrons and pro- 
tons were done^. Twenty years later, we know that the 
minimal monopole corrections proposed in his work are 
sufficient to provide results of a quality comparable to 
those of USD for nuclei in fj /2 region. 

However, there wore still problems with the monopole 
way: they showed around and beyond ^^Ni, as well as in 
the impossibility to be competitive with USD in the sd 
shell (note that USD contains non-monopole two-body 
corrections to the realistic matrix elements). The solu- 
tion came about only recently with the introduction of 
three-body forces: This far reaching development will 
be explained in detail in Sections II, II. B and V.A. We 
will only anticipate on the conclusions of these sections 
through two syllogisms. The first is: The case for real- 
istic two-body interactions is so strong that we have to 
accept them as they are. Little is known about three- 
body forces except that they exist. Therefore problems 
with calculations involving only two body forces must be 
blamed on the absence of three body forces. This argu- 



^ Realistic interactions are those consistent with data obtained in 
two (and nowadays three) nucleon systems. 



^ Two months after his PhD, Pasquini and his wife "disappeared" 
in Argentina. This dramatic event explains why his work is only 
publicly available in condensed form (Pasquini and Zuker, 1978). 
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ment raises (at least) one question: What to make of all 
calculations, using only two-body forces, that give satis- 
factory results? The answer lies in the second syllogism: 
All clearly identifiable problems are of monopole origin. 
They can be solved reasonably well by phenomenologi- 
cal changes of the monopole matrix elements. Therefore, 
fitted interactions can differ from the realistic ones basi- 
cally through monopole matrix elements. (We speak of 
R-compatibility in this case) . The two syllogisms become 
fully consistent by noting that the inclusion of three-body 
monopole terms always improves the performance of the 
forces that adopt the monopole way. 

For the Ohio spaces, the Cohen Kurath (Cohen and 
Kurath, 1965), the Chung Wildenthal (Wildenthal and 
Chung, 1979) and the FPD6 (Richter et al, 1991) inter- 
actions turn out to be R-compatible. The USD (Wilden- 
thal, 1984) interaction is R-incompatible. It will be of 
interest to analyze the recent set of p/-shell matrix ele- 
ments (GXPFl), proposed by Honma et al. (2002), but 
not yet released. 



C. The valence space 

The choice of the valence space should reflect a basic 
physical fact: that the most significant components of 
the low lying states of nuclei can be accounted for by 
many-body states involving the excitation of particles in 
a few orbitals around the Fermi level. The history of 
the Shell Model is that of the interplay between experi- 
ment and theory to establish the validity of this concept. 
Our present understanding can be roughly summed up by 
saying that "few" mean essentially one or two contiguous 
major shells^. 

For a single major shell, the classical Ohoj spaces, exact 
solutions are now available from which instructive con- 
clusions can be drawn. Consider, for instance, the spec- 
trum of ^^Ca from which we want to extract the single- 
particle states of the pf shell. Remember that in ^Li we 
expected to find two states, and indeed found two. Now, 
in principle, we expect four states. They are certainly in 
the data. However, they must be retrieved from a jun- 
gle of some 140 levels seen below 7 MeV. Moreover, the 
f'j/2 ground state is the only one that is a pure single 
particle state, the other ones are split, even severely in 
the case of the highest (Endt and van der Leun, 

1990). (Section V.B.I contains an interesting example of 
the splitting mechanism). Thus, how can we expect the 
p f shell to be a good valence space? For few particles 
above **°Ca it is certainly not. However, as we report 
in this work, it turns out that above A = 46, the lowest 
{pf)^ configurations ^ are sufficiently detached from all 



* Three major shells are necessary to deal with super-deformation. 
^ A configuration is a set of states having fixed number of particles 
in each orbit. 



others so as to generate wavefunctions that can evolve to 
the exact ones through low order perturbation theory. 

The ultimate ambition of Shell Model theory is to get 
exact solutions. The ones provided by the sole valence 
space are, so to speak, the tip of the iceberg in terms 
of number of basic states involved. The rest may be so 
well hidden as to make us believe in the literal valid- 
ity of the shell model description in a restricted valence 
space. For instance, the magic closed shells are good va- 
lence spaces consisting of a single state. This does not 
mean that a magic nucleus is 100% closed shell: 50 or 
60% should be enough, as we shall see in Section IV.B.l. 
In section V.C.I it will be argued that the Ohio valence 
spaces account for basically the same percentage of the 
full wavefunctions. Conceptually the valence space may 
be thought as defining a representation intermediate be- 
tween the Schrodinger one (the operators are fixed and 
the wavefunction contains all the information) and the 
Heisenberg one (the reverse is true). 

At best, Ohu) spaces can describe only a limited number 
of low-lying states of the same ( "natural" ) parity. Two 
contiguous major shells can most certainly cope with all 
levels of interest but they lead to intractably large spaces 
and suffer from a "Center of Mass problem" (analyzed in 
Appendix C). Here, a physically sound pruning of the 
space is suggested by the IPM. What made the success 
of the model is the explanation of the observed magic 
numbers generated by the spin-orbit term, as opposed 
to the HO (harmonic oscillator) ones (for shell forma- 
tion see Section II. B. 2). If we separate a major shell 
as HO(p)= p> ® Tp, where p> is the largest subshell 
having j = p + 1/2, and Vp is the "rest" of the HO p- 
shell, we can define the EI (extruder-intruder) spaces as 
EI(p)= rp ® (]7 + 1)>: The {p + 1)> orbit is expelled 
from HO(p + 1) by the spin-orbit interaction and in- 
trudes into HO(p). The EI spaces are well established 
standards when only one "fluid" (proton or neutron) is 
active: The Z or N = 28, 50, 82 and 126 isotopes or 
isotones. These nuclei are "spherical", amenable to ex- 
act diagonalizations and fairly well understood (see for 
example Abzouzi et al., 1991). 

As soon as both fluids are active, "deformation" ef- 
fects become appreciable, leading to "coexistence" be- 
tween spherical and deformed states, and eventually to 
dominance of the latter. To cope with this situation 
we propose "Extended EI spaces" defined as EEI(p)= 
rp ©Ap+i, where Ap = p>©(p> — 2)®. . ., i.e., the Aj = 2 
sequence of orbits that contain p> , which are needed to 
account for rotational motion as explained in Section VI. 
The EEI(l) (pi/2, ^5/2, S1/2) space was successful in de- 
scribing the full low lying spectra for A = 15-18 (Zuker 
et al. 1968, 1969). It is only recently that the EEI(2) 
(si/2, a'3/2! ./7/2, P3/2) space (region around *°Ca) has 
become tractable (see Sections VI. D and VIII. C). EEI(3) 
is the natural space for the proton rich region centered 
in ^°Zr. For heavier nuclei exact diagonalizations are not 
possible but the EEI spaces provide simple (and excel- 
lent) estimates of quadrupole moments at the beginning 
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of the well deformed regions (Section VI. C). 

As we have presented it, the choice of valence space 
is primarily a matter of physics. In practice, when 
exact diagonalizations are impossible, trnncations are 
introduced. They may be based on systematic ap- 
proaches, such as approximation schemes discussed in 
Sections IV.B.l and IV. B. 2 or mean field methods men- 
tioned in Section LB. 2, and analyzed in Section VI. A. 



D. About this review: the unified view 

We expect this review to highlight several unifying as- 
pects common to the most recent successful shell model 

calculations: 

a) An effective interaction connected with both the two 
and three nucleon bare forces. 

b) The explanation of the global properties of nuclei 
via the monopole Hamiltonian. 

c) The universality of the multipole Hamiltonian. 

d) The description of the collective behavior in the 
laboratory frame, by means of the spherical shell model. 

e) The description of resonances using the Lanczos 
strength function method. 

Let us see now how this review is organized along these 
lines. References are only given for work that will not be 
mentioned later. 

Section II. The basic tool in analyzing the interac- 
tion is the monopolc-multipolc separation. The for- 
mer is in charge of saturation and shell properties; it 
can be thought as the correct generalization of Eq. (1). 
Monopole theory is scattered in many references. Only 
by the inclusion of three-body forces could a satisfac- 
tory formulation be achieved. As this is a very recent 
and fundamental development which makes possible a 
unified viewpoint of the monopole field concept, this sec- 
tion is largely devoted to it. The "residual" multipole 
force has been extensively described in a single reference 
which will be reviewed briefly and updated. The aim of 
the section is to show how the realistic interactions can 
be characterized by a small number of parameters. 

Section III. The ANTOINE and NATHAN codes have 
made it possible to evolve from dimensionalities d„i ^ 
5 X 10-* in 1989 to drri ~ 10^, dj ~ 10^ nowadays. 
Roughly half of the-order of magnitude-gains are due 
to increases in computing power. The rest comes from 
algorithmic advances in the construction of tridiagonal 
Lanczos matrices that will be described in this section. 

Section IV. The Lanczos construction can be used to 
eliminate the "black box" aspect of the diagonalizations, 
to a large extent. We show how it can be related to 
the notions of partition function, evolution operator and 
level densities. Furthermore, it can be turned into a pow- 
erful truncation method by combining it to coupled clus- 
ter theory. Finally, it describes strength functions with 
maximal efficiency. 

Section V. After describing the three body mechanism 
which solves the monopole problem that had plagued 



the classical Ohuj calculations, some selected examples of 
pf shell spectroscopy arc presented. Special attention 
is given to Gamow Teller transitions, one of the main 
achievements of modern SM work. 

Section VI. Another major recent achievement is the 
shell model description of rotational nuclei. The new 
generation of gamma detectors, Euroball and Gamma- 
sphere has made it possible to access high spin states 
in medium mass nuclei for which full Ofuv calculations 
are available. Their remarkable harvest includes a large 
spectrum of collective manifestations which the spheri- 
cal shell model can predict or explain as, for instance, 
deformed rotors, backbending, band terminations, yrast 
traps, etc. Configurations involving two major oscillator 
shells are also shown to account well for the appearance 
of superdeformed excited bands. 

Section VII. A conjunction of factors make light and 
medium-light nuclei near the neutron drip line specially 
interesting: (a) They have recently come under intense 
experimental scrutiny, (b) They are amenable to shell 
model calculations, sometimes even exact no-core ones 
(c) They exhibit very interesting behavior, such as ha- 
los and sudden onset of deformation, (d) They achieve 
the highest N/Z ratios attained, (e) When all (or most 
of) the valence particles are neutrons, the spherical shell 
model closures, dictated by the isovector channel of the 
nuclear interaction alone, may differ from those at the 
stability valley. These regions propose some exacting 
tests for the theoretical descriptions, that the conven- 
tional shell model calculations have passed satisfactorily. 
For nearly unbound nuclei, the SM description has to 
be supplemented by some refined extensions such as the 
shell model in the continuum, which falls outside the 
scope of this review. References to the subject can be 
found in Bennaceur et al. (1999, 2000), and in Id Bctan 
et al. (2002); Michel et al. (2002, 2003) that deal with 
the Gamow shell model. 

Section VIII. There is a characteristic of the Shell 
Model we have not yet stressed: it is the approach to 
nuclear structure that can give more precise quantita- 
tive information. SM wave functions are, in particular, 
of great use in other disciplines. For example: Weak 
decay rates are crucial for the understanding of several 
astrophysical processes, and neutrinoless /?/? decay is one 
of the main sources of information about the neutrino 
masscis. In both cases, shell model calculations play a 
central role. The last section deal with these subjects, 
and some others... 

Appendix B contains a full derivation of the general 
form of the monopole field. 

Two recent reviews by Brown (2001) and Otsuka et al. 
(2001b) have made it possible to simplify our task and 
avoid redundancies. However we did not feel dispensed 
from quoting and commenting in some detail important 
work that bears directly on the subjects we treat. 

Sections II, IV and the Appendices are based on un- 
published notes by APZ. 
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II. THE INTERACTION 

The fohowing remarks from Abzouzi, Caurier, and 
Zuker (1991) still provide a good introduction to the sub- 
ject: 

"The use of realistic potentials {i.e., consistent with 
NN scattering data) in shcU-modcl calculation was pi- 
oneered by Kuo and Brown (1966). Of the enormous 
body of work that followed wc would like to extract two 
observations. The first is that whatever the forces (hard 
or soft core, ancient or new) and the method of regu- 
larization (Brucckncr G matrix (Kahana et ai, 1969a; 
Kuo and Brown, 1966), Sussex direct extraction (Elliott 
et al., 1968) or Jastrow correlations (Fiase et al, 1988)) 
the effective matrix elements are extraordinarily simi- 
lar (Pasquini and Zuker, 1978; Rutsgi et al, 1971). The 
most recent results (Jiang et al, 1989) amount to a vin- 
dication of the work of Kuo and Brown. We take this 
similarity to be the great strength of the realistic inter- 
actions, since it confers on them a model-independent 
status as direct links to the phase shifts. 

The second observation is that when used in shell- 
model calculations and compared with data these matrix 
elements give results that deteriorate rapidly as the num- 
ber of particle increases (Halbert et al, 1971) and (Brown 
and Wildenthal, 1988). It was found (Pasquini and 
Zuker, 1978) that in the pf shell a phenomenologi- 
cal cure, confirmed by exact diagonalizations up to 
A=48 (Caurier et ai, 1994), amounts to very simple mod- 
ifications of some average matrix elements (centroids) of 
the KB interaction (Kuo and Brown, 1968)." 

In Abzouzi et al. (1991) very good spectroscopy in the 
p and sd shells could be obtained only through more rad- 
ical changes in the centroids, involving substantial three- 
body (3b) terms. In 1991 it was hard to interpret them 
as effective and there were no sufficient grounds to claim 
that they were real. 

Nowadays, the need of true 3b forces has become ir- 
refutable: In the 1990 decade several two-body (2b) po- 
tentials were developed Nijmegen I and II (Stoks et al, 
1993), AV18 (Wiringa et ai, 1995), CD-Bonn- (Machleidt 
et al., 1996)— that fit the w 4300 entries in the Ni- 
jmegen data base (Stoks et al, 1994) with « 1, and 
none of them seemed capable to predict perfectly the 
nucleon vector analyzing power in elastic {N, d) scatter- 
ing (the Ay puzzle). Two recent additions to the fam- 
ily of high-precision 2b potentials the charge-dependent 
"CD-Bonn" (Machleidt, 2001) and the chiral Idaho-A 
and B (Entem and Machleidt, 2002) — have dispelled any 
hopes of solving the Ay puzzle with 2b-only interac- 
tions (Entem et al., 2002). 

Furthermore, quasi-exact 2b Green Function Monte 
Carlo (GFMC) results (Pudliner et al, 1997; Wiringa 
et al, 2000) which provided acceptable spectra for ^ < 8 
(though they had problems with binding energies and 
spin-orbit splittings), now encounter serious trouble in 
the spectrum of -'^"B, as found through the No Core Shell 
Model (NCSM) calculations of Navratil and Ormand 



(2002) and Caurier et al. (2002), confirmed by Pieper, 
Varga, and Wiringa (2002), who also show that the prob- 
lems can be remedied to a large extent by introducing 
the new Illinois 3b potentials developed by Pieper et al 
(2001). 



It is seen that the trouble detected with a 2b-only 
description — with binding energies, spin-orbit splittings 
and spectra — is always related to centroids, which, once 
associated to operators that depend only on the num- 
ber of particles in subshells, determine a "monopole" 
Hamiltonian Hm that is basically in charge of Hartree 
Fock selfconsistency. As we shall show, the full H can be 
separated rigorously as H = Tim + 'Hm- The multipole 
part Hm includes pairing, quadrupole and other forces 
responsible for collective behavior, and — as checked by 
many calculations — is well given by the 2b potentials. 



The preceding paragraph amounts to rephrasing the 
two observations quoted at the beginning of this section 
with the proviso that the blame for discrepancies with 
experiment cannot be due to the — now nearly perfect — 
2b potentials. Hence the necessary "corrections" to Hm 
must have a 3b origin. Given that we have no complaint 
with Hm, the primary problem is the monopole contri- 
bution to the 3b potentials. This is welcome because 
a full 3b treatment would render most shell model cal- 
culations impossible, while the phenomenological study 
of monopole behavior is quite simple. Furthermore it is 
quite justified because there is little ah initio knowledge 
of the 3b potentials®. Therefore, whatever information 
that comes from nuclear data beyond A=3 is welcome. 



In Shell Model calculations, the interaction appears as 
matrix elements, that soon become far more numerous 
than the number of parameters defining the realistic po- 
tentials. Our task will consist in analyzing the Hamilto- 
nian in the oscillator representation (or Fock space) so as 
to understand its workings and simplify its form. In Sec- 
tion II. A we sketch the theory of effective interactions. In 
Section II. B it is explained how to construct from data a 
minimal Hm, while Section II. C will be devoted to extract 
from realistic forces the most important contributions to 
Hm- The basic tools are symmetry and scaling argu- 
ments, the clean separation of bulk and shell effects and 
the reduction of H to sums of factorable terms (Dufour 
and Zuker, 1996). 



Experimentally there will never be enough data to determine 
them. The hope for a few parameter description eomes from 
chiral perturbation theory (Entem and Machleidt, 2002). 
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A. Effective interactions 



The Hamiltonian is written in an oscillator basis as 



= K,+ ^ ^Istu^tsT ■ -^t^r 
r<s, t<«, r 

+ ^rstuvw'^tstT ' ^uvwT, (3) 

r<s<i, ti<D<t(;, r 

where K, is the kinetic energy, V^^/ the interaction matrix 
elements, Z+p ( Z^r) create (annihilate) pairs (r = rs) 
or triples (r = rst) of particles in orbits r, coupled to 
r = JT. Dots stand for scalar products. The basis and 
the matrix elements are large but never infinite^. 

The aim of an effective interaction theory is to reduce 
the secular problem in the large space to a smaller model 
space by treating perturbatively the coupling between 
them, thereby transforming the full potential, and its re- 
pulsive short distance behavior, into a smooth pseudopo- 
tential. 

In what follows, if we have to distinguish between 
large {Nfko) and model spaces we use H, K, and V for the 
pseudo-potential in the former and H, K andV for the 
effective interaction in the latter. 



1. Theory and calculations 

The general procedure to describe an exact eigen- 
state in a restricted space was obtained independently 
by Suzuki and Lee (1980) and Poves and Zuker (1981a)^. 
It consists in dividing the full space into model and 
external (|a)) determinants, and introducing a transfor- 
mation that respects strict orthogonality between the 
spaces and decouples them exactly: 

\T) = \i) + ^Aioc\a) = \a) - ^Aio,\i) (4) 

a i 

ty\a) = 0, Aia is defined through = 0. (5) 

The idea is that the model space can produce one or 
several starting wavefunctions that can evolve to exact 
eigenstates through perturbative or coupled cluster eval- 
uation of the amplitudes Aia, which can be viewed as 
matrix elements of a many body operator A. In cou- 
pled cluster theory (CCT or exp S) (Coester and Kiimmel 
(1960), Kiimmel et al. (1978)) one sets A = exp S, where 



Nowadays, the non-relativistic potentials must be thought of as 
derived from an effective field theory which has a cutoff of about 
1 GeV (Entem and Machleidt, 2002). Therefore, truly hard cores 
are ruled out, and 7i should be understood to act on a sufficiently 
large vector space, not over the whole Hilbert space. 
* The notations in both papers are very different but the pertur- 
bative expansions are probably identical because the Hermitean 
formulation of Suzuki (1982) is identical to the one given in Poves 
and Zukcr (1981a). This paper also deals extensively with the 
coupled cluster formalism. 



S = Si + S2 + ■ ■ ■ + Sk is a sum of fc-body (kb) opera- 
tors. The decoupling condition (z|7i|a) = then leads to 
a set of coupled integral equations for the Si amplitudes. 
When the model space reduces to a single determinant, 
setting S'l = leads to Hartree Fock theory if all other 
amplitudes are neglected. The S2 approximation con- 
tains both low order Brueckner theory (LOBT) and the 
random phase approximation (RPA) . In the presence 
of hard-core potentials, the priority is to screen them 
through LOBT, and the matrix elements contributing 
to the RPA are discarded. An important implementa- 
tion of the theory was due to Zabolitzky, whose calcu- 
lations for -^He, and '*"Ca, included 5.3 (Bcthe Fad- 
deev) and S4 (Day Yacoubovsky) amplitudes (Kiimmel, 
Liihrmann, and Zabolitzky, 1978). This "Bochum trun- 
cation scheme" that retraces the history of nuclear mat- 
ter theory has the drawback that at each level terms that 
one would like to keep are neglected. 

The way out of this problem (Heisenberg and Mi- 
haila, 2000; Mihaila and Heisenberg, 2000) consists on 
a new truncation scheme in which some approximations 
are made, but no terms are neglected, relying on the fact 
that the matrix elements are finite. The calculations of 
these authors for (up to ^3) can be ranked with the 
quasi-exact GFMC and NCSM ones for lighter nuclei. 

In the quasi-degenerate regime (many model states), 
the coupled cluster equations determine an effective in- 
teraction in the model space. The theory is much sim- 
plified if we enforce the decoupling condition for a single 
state whose exact wavefunction is written as 



\ref) = (I + A1+A2 + ...)(! + -B1 
= {l + Ci+C2 + ...)\ref), 



-B2 + ...)\ref} 
C = expS, (6) 



where \ref) is a model determinant. The internal am- 
plitudes associated with the B operators arc those of an 
eigcnstate obtained by diagonalizing Heff — Ti,{l + A) 
in the model space. As it is always possible to eliminate 
the Al amplitude, at the 5*2 level there is no coupling 
i.e., the effective interaction is a state independent G- 
matrix (Zukcr, 1984), which has the advantage of pro- 
viding an initialization for J?e//- Going to the ^3 level 
would be very hard, and we examine what has become 
standard practice. 

The power of CCT is that it provides a unified frame- 
work for the two things we expect from decoupling: to 
smooth the repulsion and to incorporate long range cor- 
relations. The former demands jumps of say, bOfuv, the 
latter much less. Therefore it is convenient to treat them 
separately. Standard practice assumes that G-matrix el- 
ements can provide a smooth pseudopotential in some 
sufficiently large space, and then accounts for long range 
correlations through perturbation theory. The equation 
to be solved is 



Gijkl — Vijkl — ■ 



a/3 



+ e/3 - Ci - Cj + A 



(7) 



where ij and kl now stand for orbits in the model space 
while in the pair aP at least one orbit is out of it, 
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is an unperturbed (usually kinetic) energy, and A a free 
parameter called the "starting" energy. Hjorth-Jensen, 
Kuo, and Osnes (1995) describe in detail a sophisticated 
partition that amounts to having two model spaces, one 
large and one small. 

In the NCSM calculations an Ntuv model space is cho- 
sen with N w 6-10. When initiated by Zheng et al. (1993) 
the pseudo-potential was a G-matrix with starting en- 
ergy. Then the A dependence was eliminated either by 
arcane perturbative maneuvers, or by a truly interest- 
ing proposal: direct decoupling of 2b elements from a 
very large space (Navratil and Barrett, 1996), further im- 
plemented by Navratil, Vary, and Barrett (2000a), and 
extended to 3b effective forces (Navratil et al., 2000), 
(Navratil and Ormand, 2002). It would be of inter- 
est to compare the resulting effective interactions to the 
Brueckner and Bethe Faddeev amplitudes obtained in a 
full cxp S approach (which are also free of arbitrary start- 
ing energies). 

A most valuable contribution of NCSM is the proof 
that it is possible to work with a pseudo-potential in Nhu 
spaces. The method relies on exact diagonalizations. As 
they soon become prohibitive, in the future, CCT may 
become the standard approach: going to 5*3 for = 50 
(as Heisenberg and Mihaila did) is hard. For A' = 10 it 
should be much easier. 

Another important NCSM indication is that the exci- 
tation spectra converge well before the full energy, which 
validates formally the Oftw diagonalizations with rudi- 
mentary potentials^ and second order corrections. 

The Otiu) results are very good for the spectra. How- 
ever, to have a good pseudopotential to describe energies 
is not enough. The transition operators also need dress- 
ing. For some of them, notably E2, the dressing mecha- 
nism (coupling to 2fiw quadrupole excitations) has been 
well understood for years (see Dufour and Zuker, 1996, 
for a detailed analysis), and yields the, abundantly tested 
and confirmed, recipe of using effective charges of « 1.5e 
for the protons and ~ 0.5e for neutrons. For the Gamow 
Teller (GT) transitions, mediated by the spin-isospin 
aT± operator the renormalization mechanism involves an 
overall "quenching" factor of 0.7-0.8 whose origin is far 
subtler. It will be examined in Sections V.C.I and V.C.2. 
The interested reader may wish to consult them right 
away. 



B. The monopole Hamiltonian 

A many-body theory usually starts by separating the 
Hamiltonian into an "unperturbed" and a "residual" 

part, Ti. = + Wr- The traditional approach consists 
in choosing for Ho a one body (lb) single-particle field. 



Since Ti. contains two and three-body (2b and 3b) compo- 
nents, the separation is not mathematically clean. There- 
fore, we propose the following 

Ti = Km + Hmi (8) 

where Hm-, the monopole Hamiltonian, contains K and 
all quadratic and cubic (2b and 3b) forms in the scalar 
products of fcrmion operators a].^ ■ a.sy^'^ ^ while the mul- 
tipole Tiu contains all the rest. 

Our plan is to concentrate on the 2b part, and intro- 
duce 3b elements as the need arises. 

Hm has a diagonal part, W^, written in terms of num- 
ber and isospin operators (oj^ • ar^)- It reproduces the 
average energies of configurations at fixed number of par- 
ticles and isospin in each orbit {jt representation) or, al- 
ternatively, at fixed number of particles in each orbit and 
each fluid (neutron-proton, np, or j representation). In 
jt representation the centroids 

^ EjV/t?t(2J+l)[l-(-)"+^M . 
E,/(2J+l)[l-(-)^+^'5st] ' ^ ^ 

= J(3Vi + V°), 6,* = V,\-V°„ (9b) 

are associated to the 2b quadratics in number (rUg) and 
isospin operators (T^), 

nist = \ ms{mt-Sst), (10a) 

J- + Ost 

= T-^(r«-rt- ^m,t<5,t), (10b) 

to define the diagonal 2b part of the monopole Hamilto- 
nian (Bansal and French (1964), French (1969)) 

nijt = + ^{ast rUst + bst Tst) + V*^, (11) 

The 2b part is a standard result, easily extended to in- 
clude the 3b term 

^stu '^stu ~l~ ^stu ^stu 

stu 

where mstu = ^'^st'nnm or ms{ms — l)(?7is — 2)/6 and 
Tstu = nT'sTtu or (mg — 2)Tsu- The extraction of the full 
Hm is more complicated and it will be given in detail 
in Appendix B. Here we only need to note that Hm is 
closed under unitary transformations of the underlying 
fermion operators, and hence, under spherical Hartree- 
Fock (HF) variation. This property explains the interest 
of the Hm + T^M separation. 



^ Even old potentials fit well the low energy NN phase shifts, the 
only that matter at OHui. 



r and s arc subslicUs of the same parity and angular momentum, 
X and y stand for neutrons or protons 
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1. Bulk properties. Factorable forms 

Hm must contain all the information necessary to pro- 
duce the parameters of the Bethe Weiszacker mass for- 
mula, and we start by extracting the bulk energy. The 
key step involves the reduction to a sum of factorable 
forms valid for any interaction (Dufour and Zuker, 1996). 
Its enormous power derives from the strong dominance 
of a single term in all the cases considered so far. 

For clarity we restrict attention to the full isoscalar 
centroids defined in Eqs. (B18a) and (B18c)^^ 



Vst — ttst — 



4(4i« + l)- 



Diagonalize Vst 



(13) 



(14) 



^L = ^^k^Usk'ms^Utkmt-^Vssms. (15) 

k s t s 

For the 3b interaction, the corresponding centroids Vstu 
are treated as explained in (Dufour and Zuker, 1996, 
Appendix Bib): The st pairs are replaced by a sin- 
gle index x. Let L and M be the dimensions of the x 
and s arrays respectively. Construct and diagonalize an 
{L + M) x (L + M) matrix whose non-zero elements are 
the rectangular matrices Vxs and V^x- Disregarding the 
contractions, the strictly 3b part V^^-* can then be written 
as a sum of factors 



(16) 



Full factorization follows by applying Eqs. (14,15) to the 
^tu matrices for each k. 



In Eq. (15) a single term strongly dominates all oth- 
ers. For the KLS interaction (Kahana et ai, 1969b) and 
including the first eight major shells, the result in Fig. 1 
is roughly approximated by 



Us 



i.-0.5{l-{l)) + {j-l) 



(17) 



where p is the principal quantum number of an oscillator 
shell with degeneracy Dp — J^sC^Js + 1) = (p+ l)(p + 2), 
(0 = K'^l + 1)/ E/(2^ + !)■ The (unitary) U matrices 
have been affected by an arbitrary factor 6 to have num- 
bers of order unity. Operators of the form {Dg = 2js + 1) 



with 



J2 Ds^s = 



(18) 



vanish at closed shells and are responsible for shell effects. 

As I — (l)) and j — I arc of this type, only the the Up part 
contributes to the bulk energy. 

To proceed it is necessary to know how interactions 
depend on the oscillator frequency oj of the basis, related 
to the observed square radius (and hence to the density) 
through the estimate (Bohr and Mottelson, 1969, Eq. (2- 
157)): 



2 3 4 5 6 7 

P 



FIG. 1 The Us terms in Eq. (17). Arbitrary scale. 



They ensure the right average energies for T = closed shells, 
which is not the case in Eq. (11) if one simply drops the bst 
terms. 



35.59 



(^)V3 (r2) 



40 

IV3 



MeV 



(19) 



A 6 force scales as (hw)^^^. A 2b potential of short range 
is essentially linear in hco; for a 3b one we shall tentatively 
assume an {hu)"^ dependence while the Coulomb force 
goes exactly as {fkoY^'^. 

To calculate the bulk energy of nuclear matter we av- 
erage out subshell effects through uniform filling => 
rUp Dg/ Dp. Though the 17-type operators vanish, we have 
kept them for reference in Eqs. (21,22) below. The lat- 
ter is an educated guess for the 3b contribution. The 
eigenvalue Eq for the dominant term in Eq. 15 is replaced 
by hioVQ defined so as to have Up = 1. The subindex 
m is dropped throughout. Then, using Boole's factorial 
powers, e.g., p'^^^ = p{p — \){p — 2), we obtain the fol- 
lowing asymptotic estimates for the leading terms (i.e.. 
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disregarding contractions) 



(p/+3)(3)(p/ + 2) 



(20) 



\ p 



+ n 



?iwVo[p/(p/+4)]2, 



^(M'/3VopKp/+4)2. 
Finally, relate pf to A : 

^mp = ^2(p+l)(p + 2) 

p p=0 




+ n2 



(3) 



(23) 



Note that in Eq. (22) we can replace (fiw)^ by {fuv)^^'^ 
and change the powers of Dp in the denominators accord- 
ingly. 

Assuming that non-diagonal IC^ and V"^ terms cancel 
we can vary with respect to hiv to obtain the saturation 
energy 

E,=)C'^ + {1- /3a;>/)V^ 

dE ^ IC^ + V^ 

Piv^pf = 



dui 



(k + 1) 

^^(1_4Vo)('^aV^'. 
4 K+1 ^ V2 / 



(24) 



The correct saturation properties are obtained by fix- 
ing Vo so that Es/A « 15.5 MeV. The TwJe « AOA''^/^ 
choice ensures the correct density. It is worth noting that 
the same approach leads to Vc « 3e^.^(.Z— l)/5i?'^ (Duflo 
and Zuker, 2002). 

It should be obvious that nuclear matter properties 
derived from finite nuclei — could be calculated with tech- 
niques designed to treat finite nuclei: A successful the- 
oretical mass table must necessarily extrapolate to the 
Bethe Weiszacker formula (Duflo and Zuker, 1995). It 
may be surprising though, that such calculations could 
be conducted so easily in the oscillator basis. The merit 
goes to the separation and factorization properties of the 
forces. 

Clearly, Eq. (24) has no (or trivial) solution for k = 0, 
i.e., without a 3b term. Though 2b forces do saturate, 
they do it at the wrong place and at a heavy price be- 
cause their short-distance repulsion prevents direct HF 
variation. The crucial question is now: Can we use real- 
istic (2-|-3)b potentials soft enough to do HF? Probably 
we can. The reason is that the nucleus is quite dilute, and 



nucleons only "see" the low energy part of the 2b poten- 
tial involving basically s wave scattering, which has been 
traditionally well fitted by realistic potentials. Which in 
turn explains why primitive versions of such potentials 
give results close to the modern interactions for G-matrix 
elements calculated at reasonable fvu) values. This is a re- 
current theme in this review (see especially Section II. D) 
and Fig. 2 provides an example of particular relevance for 
the study of shell effects in next Section II. B. 2. Before 
we move on, we insist on what we have learned so far: It 
may be possible to describe nuclear structure with soft 
(2+3)b potentials consistent with the low energy data 
coming from the A = 2 and 3 systems. 
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FIG. 2 The KLS (Kahana et al., 1969b) shell effects in Fig. 1 
for the first 10 orbits compared to those produced by a Bonn 

potential (Hjorth- Jensen, 1996) (to within an overall factor). 
Points in the same major shell are connected by lines and 
ordered by increasing j. 



2. Shell formation 

Fig. 2 provides a direct reading of the expected sin- 
gle particle sequences produced by realistic interactions 
above T = closed shells. Unfortunately, they do not 
square at all with what is seen experimentally. In partic- 
ular, the splitting between spin-orbit partners appears to 
be nearly constant, instead of being proportional to I. As 
it could no longer be claimed that the 2b realistic forces 
are "wrong" , the solution must come from the 3b oper- 
ators. The "educated guess" in Eq. (22) may prove of 
help in this respect but it has not been implemented yet. 
What we propose instead is to examine an existing, strict 
(l-|-2)b, model of Ti."^ that goes a long way in explaining 
shell formation. Though we know that 3b ingredients are 
necessary, thciv may be often mocked by lower rank ones, 
in analogy with the kinetic energy, usually represented 
as lb, when in reality it is a 2b operator, as explained in 
Appendix C. Something similar seems to happen with 
the spin-orbit force which must have a 3b origin, while 
a lb picture is phenomenologically quite satisfactory, as 
our study will confirm. We shall also find out that some 
mechanisms have an irreducible 3b character. 
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7^„j manifests itself directly in nuclear spectra through 
the doubly magic closures and the particle and hole states 
built on them. The members of this set, which we call 
cs ± 1 are well represented by single determinants. Con- 
figuration mixing may alter somewhat the energies but, 
most often, enough experimental evidence exists to cor- 
rect for the mixing. Therefore, a minimal characteriza- 
tion, TYJ^, can be achieved by demanding that it repro- 
duce all the observed cs ± 1 states. The task was under- 
taken by Duflo and Zuker (1999). With a half dozen pa- 
rameters a fit to the 90 available data achieved an rmsd of 
some 220 keV. The differences in binding energies (gaps) 
2BE{cs) - BE{cs + 1)- BE{cs - 1), not included in the 
fit, also came out quite well and provided a test of the 
reliability of the results. 

The model concentrates on single particle splittings 
which are of order A^^^"^. The valence spaces involve a 
number of particles of the order of their degeneracy, i.e., 
« {pf + 3/2)2 ^ (3^/2)2/3 according to Eq. (23). As a 
consequence, at midshell, the addition of single particle 
splittings generates shell effects of order A-^^^. To avoid 
bulk contributions all expressions are written in terms of 
ri-type operators [Eq. (18)] and the model is defined by 




ni = {W- AK) + n{ls) + Cl{ll) 
+n(ffi) + Clizni) + n(ffc) + Cl(znc), 



(25) 



FIG. 3 Evolution of cs ± 1 states from '"'Ca to "^^Ca 



Where: 



[I] W-AK 



3/2). (26) 



comes from IC^ and V'. Both terms lead to the same 
gaps, and the chosen combination has the advantage of 
canceling exactly to orders A and A"^/^. 

[II] Cl{ls) « -22 l-s/A^^^ (very much the value in (Bohr 
and Mottelson, 1969, Eq. (2-132)) and Cl{ll) w -22[l{l + 
l)-p{p+3)/2]/A'^/^ (in MeV), are Cl-type operators that 
reproduce the single-particle spectra above HO closed 
shells. 

[III] Cl{ffi), Cl{zni), Cl{ffc) and ^(^nc)^^ are fi-type 
2b operators whose effect is illustrated in Fig. 3: In going 
from ''"Ca to '^^Ca the filling of the f'j/2 (/ for short) 
neutron orbit modifies the spectra in ^^Ca, ^^Sc, ^^Ca 
and '^^K. 

As mentioned, realistic 2b forces will miss the particle 
spectrum in "^^Ca. Once it is phenomenologically cor- 
rected hy Is + II such forces account reasonably well for 
the zni and eei mechanism, within one major exception: 
they fail to depress the /7/2 orbit with respect to the 
others. As a consequence, and quite generally, they fail 



to produce the EI closures. Though zni and eei can be 
made to cure this serious defect, it is not clear how well 
they do it, since the correct mechanism necessarily in- 
volves 3b terms, as will be demonstrated in Section V.A. 
Another problem that demands a three body mechanism 
is associated with the znc and ffc operators which al- 
ways depress orbits of larger I [as in the hole spectra of 
''^Ca and ''^K in Fig. 3], which is not the case for light 
nuclei as seen from the spectra of ^^C and ^^Si. Only 3b 
operators can account for this inversion. 

In spite of these caveats, should provide a plausible 
view of shell formation, illustrated in Figure 4 for t = 




40 50 60 



// stands for "same fluid", zn for "proton-neutron", i for "intra- 
shell" and c for "cross-shell" 



FIG. 4 The different contributions to Ht for N=Z nuclei 
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FIG. 5 H„, for even t = N - Z =10. 
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N—Z = nuclei. 7i^„ is calculated by filling the oscillator 
orbits in the order dictated by Is + II. The W — AK 
term produces enormous HO closures (showing as spikes). 
They are practically erased by the lb contributions, and 
it takes the 2b ones to generate the EI-'^^ spikes. Note that 
the HO magicity is strongly attenuated but persists. The 
drift toward smaller binding energies that goes as A^^^ is 
an artifact of the W — AK choice and should be ignored. 

In Figs. 5 we show the situation for even t = N— Z = 
to 18. Plotting along lines of constant t has the advan- 
tage of detecting magicity for both fluids. In other words, 
the spikes appear when either fluid is closed, and are re- 
inforced when both are closed. The spikes are invari- 
ably associated either to the HO magic numbers (8, 20, 
40) or to the EI ones (14, 28, 50), but the latter always 
show magicity while the former only do it at and above 
the double closures (Z, 7V)=(8, 14), (20,28) and (40,50). 
For example, at Z = 20 "^^Ca shows as a weak closure. 
42,44, 46 g^j.g j^Q^ closed (which agrees with experiment), 
^®Ca is definitely magic, and spikes persist for the heavier 
isotopes. At Z = 40, *°Zr shows a nice spike, a bad pre- 



Remember: EI valence spaces consist in the HO orbits shell p, 
except the largest (extruder), plus the largest (intruder) orbit 
from shell p -|- 1. 



diction for a strongly rotational nucleus. However, there 
are no (or weak) spikes for the heavier isotopes except at 
^"Zr and ^^Zr, both definitely doubly magic. There are 
many other interesting cases, and the general trend is the 
following: No known closure fails to be detected. Con- 
versely, not all predicted closures are real. They may be 
erased by deformation — as in the case of ^°Zr — but this 
seldom happens, thus suggesting fairly reliable predictive 
power for 7i^. 



C. The multipole Hamiltonian 

The multipole Hamiltonian is defined as Ti.M—'H-Ti.m- 
As we are no longer interested in the full TC, but its re- 
striction to a finite space, Hm will be more modestly 
called Hm , with monopole-free matrix elements given by 



JT 



JT 



SrsStuVr, 



(27) 



We shall describe succinctly the main results of Du- 
four and Zuker (1996), insisting on points that were not 
stressed sufficiently in that paper (Section II.C.l), and 
adding some new information (Section II. D). 
There are two standard ways of writing Hm- 



E 



r<s.t<u.r 



H. 



M 



rstuT 



r ii/2 f-r (qi c7 \0 

III J rtsuy^rt^ su } ' 



(28) 



(29) 



where /, 



7 

rtsu 
7 



+ Srs){l + Stu)/^. The matrix 
elements w^tsu ^^'^ ^rstu related through Eqs. (A9) 
and (AlO). 

Replacing pairs by single indices rs = x, tu = y in 
eq. (28) and rt = a, su = b in eq. (29), we proceed 
as in Eqs. (14) and (15) to bring the matrices W^y = 
W^gf.^ and = frtsu I to diagonal form through unitary 



transformations U^i^, i 



H 



M 



^j. to obtain factorable expressions 
Y^U^^Zyr, (30) 



fe,r 




1/2 



(31) 



which we call the E (or normal or particle-particle or 
pp) and e (or multipole or particle-hole or ph) represen- 
tations. Since Tim contains all the 7 = 00 and 01 terms, 
for Hm, i^rltu = ^rltu = (see Eq. (Bll). There are 
no one body contractions in the e representation because 
they are all proportional to 

The eigensolutions in eqs. (30) and (31) using the KLS 
interaction (Kahana et ai, 1969b; Lee, 1969), yield the 
density of eigenvalues (their number in a given interval) 
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FIG. 6 i?-eigenvalue density for the KLS interaction in the 

pf+sdg major shells hu = 9. Each eigenvalue has multiplicity 
[r]. The largest ones are shown by arrows. 



in the E representation that is shown in Fig. 6 for a 
typical two-shell case. 

It is skewed, with a tail at negative energies which is 
what we expect from an attractive interaction. 

The e eigenvalues are plotted in fig. 7. They are very 
symmetrically distributed around a narrow central group, 
but a few of them are neatly detached. The strongest 
have 7'^ = 1-0, 1+1, 2+0, 3^0, 4+0. "If the corre- 
sponding eigenvectors arc eliminated from H in" cq. (31) 
and the associated H in eq. (30) is recalculated^'*, the E 
distribution becomes quite symmetric, as expected for a 
random interaction. 

If the diagonalizations are restricted to one major shell, 
negative parity peaks are absent, but for the positive par- 
ity ones the results arc practically identical to those of 
Figs. 6 and 7, except that the energies are halved. This 
point is crucial: 

J/Wpj andup^ are the eigenvectors obtained in shells pi 
and p2, their eigenvalues are approximately equal Cp^ « 
= e. When diagonalizing inpi+p2, the unnormalized 
eigenvector turns out to be Up^ + Up^ with eigenvalue e. 

In the figures the eigenvalues for the two shell case 
are doubled, because they are associated with normalized 
eigenvectors. To make a long story short: The contribu- 
tion to Hm associated to the large F = 01, and -j — 20 



Inverted commas arc meant to call attention on an erratum 
in (Dufour and Zuker, 1996). 
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FIG. 7 e-eigenvalue density for the KLS interaction in the 
pf-l-sdg major shells. Each eigenvalue has multiplicity [7]. 
The largest ones are shown by arrows. 



terms, 



H 



E'''\{Pl + Pl^,).{Pp + Pp+,) (32) 
Hq = -T^\e^°\{qp + qp+i)- {qp + Qp+i), (33) 



turns out to be the usual pairing plus quadrupole Hamil- 
tonians, except that the operators for each major shell 
of principal quantum number p are affected by a nor- 
malization. E'^^ and e^" are the one shell values called 
generically e in the discussion above. To be precise 



^;=E^-ol^^yv^^y^ (34) 

Qp=Yl SrsQrs/K, (35) 

r 



where fir = jr + ^, Q-, 



{r\\r''Y^\\s) (36) 



and np = ^Dp, A^ = Eql^^{p + 3/2)^ (37) 



1. Collapse avoided 

The pairing plus cjuadnipole (P -I- Q) model has a long 
and glorious history (Barangcr and Kumar, 1968; Bcs and 
Sorensen, 1969), and one big problem: as more shells are 
added to a space, the energy grows, eventually leading 
to collapse. The only solution was to stay within lim- 
ited spaces, but then the coupling constants had to be 
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readjusted on a case by case basis. The normalized ver- 
sions of the operators presented above are affected by 
universal coupling constants that do not change with the 
number of shells. Knowing that fuvo = 9 MeV, they are 



\/hwo 



0.216 in 



\E°^\/huJo = g' = 0.32 and je^o 
Eqs. (32) and (33). 

Introducing A^f « |(p/ + 3/2)'^, the total number of 
particles at the middle of the Fermi shell p/, the rela- 
tionship between g' , k', and their conventional counter- 
parts (Baranger and Kumar, 1968) is, for one shell^^ 

0.32fia; 19.51 ^ ^ . i 



Op 
0.216/ia; 



1 216 



2 AV3Al/| 



_ X_ _ Xo^-5/3 



(38) 



To see how collapse occurs assume m = 0{Df) = 

(9 (^2/3 -J j-^ ^Yic Fermi shell, and promote them to a higher 
shell of degeneracy D. The corresponding pairing and 
quadrupole energies can be estimated as 

(39) 
(40) 

respectively, which become for the two possible scalings 
D 



Ep = -\G\4m{D -m + 2) = -\G\0{mD), 
and E, « -Ix'IQ^ = -\x'\0{m^D), 



Ep{old) = 0{- 



> Ep(new) = 0{ 
Eq{new) = 0{ 



m 

2 



) (41) 

MV3^-)- (42) 

If the TO particles stay at the Fermi shell, all energies go 
as A^/^ as they should. If D grows both energies grow 
in the old version. For sufficiently large D the gain will 
become larger than the monopole loss 0{mD^^'^hu]) = 
0{D^/'^A^/^). Therefore the traditional forces lead the 
system to collapse. In the new form there is no collapse; 
Ep stays constant, Eq decreases and the monopole term 
provides the restoring force that guarantees that particles 
will remain predominantly in the Fermi shell. 

As a model for TLm, P + Q is likely to be a reasonable 
first approximation for many studies, provided it is sup- 
plemented by a reasonable to produce the P+Q+m 
model. 



D. Universality of tlie realistic interactions 

To compare sets of matrix elements we define the over- 
laps 

OaB =Va-Vb=Y1 yrstuAV^stuB (43) 



Oal 



rstuT 

Oab 



VOaaObb 



(44) 



According to Dufour and Zukcr (1996) the bare coupling con- 
stants, should be renormalized to 0.32(1-1-0.48) and 0.216(1-1-0.3) 



Similarly, is the overlap for matrix elements with the 
same T. Oaa is a measure of the strength of an inter- 
action. If the interaction is referred to its centroid (and 
a fortiori if it is monopole- free), Oaa = (^a- The upper 



TABLE I The overlaps between the KLS(A) and Bonn(B) 
interactions for the first ten oscillator shells (2308 matrix el- 
ements), followed by those for the pf shell (195 matrix ele- 
ments) for the BonnC (C) and KB (K) interactions. 



O'kA 


= 17.56 


0\a 


= 2.61 


Oaa 


= 6.49 




= 11.84 


Ohs 


= 2.31 


Obb 


= 4.78 


o% 


= 0.98 


0\b 


= 0.99 


Oab 


= 0.98 


0%c 


= 3.71 


Ohc 


= 0.56 


Occ 


= 1.41 




= 3.02 




= 0.46 


Okk 


= 1.15 




= 0.99 


OcK 


= 0.97 


OcK 


= 0.99 



set of numbers in Table I contains what is probably the 
most important single result concerning the interactions: 
a 1969 realistic potential (Kahana et al. (1969b); Lee 
(1969)) and a modern Bonn one (Hjorth-Jensen, 1996) 
differ in total strength(s) but the normalized 0{AB) are 
very close to unity. 

Now: all the modern realistic potentials agree closely 
with one another in their predictions except for the bind- 
ing energies. For nuclear matter the differences are sub- 
stantial (Pieper et al., 2001), and for the BonnA,B,C 
potentials studied by Hjorth-Jensen et al. (1995), they 
become enormous. However, the matrix elements given 
in this reference for the pf shell have normalized cross 
overlaps of better than 0.998. At the moment, the overall 
strengths must be viewed as free parameters. 

At a fundamental level the discrepancies in total 
strength stem for the degree of non-locality in the po- 
tentials and the treatment of the tensor force. In the old 
interactions, uncertainties were also due to the starting 
energies and the renormalization processes, which, again, 
affect mainly the total strength(s), as can be gathered 
from the bottom part of Table I. 

The dominant terms of Hm are central and table II 
collects their strengths (in MeV) for effective interac- 
tions in the j3/-shell: Kuo and Brown (1968) (KB), the 
potential fit of Richter et al. (1991) (FPD6), the Gogny 
force — successfully used in countless mean field studies 
(Decharge and Gogny, 1980) and BonnC (Hjorth-Jensen 
et al, 1995). 

There is not much to choose between the different 
forces which is a nice indication that overall nuclear 
data (used in the FPD6 and Gogny fits) are consistent 
with the A'' A'' data used in the realistic KB and BonnC 
G-matrices. The splendid performance of Gogny de- 
serves further study. The only qualm is with the weak 
quadrupole strength of KB, which can be understood by 
looking again at the bottom part of Table I. In assess- 
ing the JT = 01 and JT =10 pairing terms, it should 
be borne in mind that their renormalization remains an 
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open question (Dufour and Zuker. 1996). 

Conclusion: Some fine tuning may be in order and 3b 
forces may (one day) bring some multipole news, but as 
of now the problem is Hm, not Hm- 

TABLE II Leading terms of tfie multipole Hamiltonian 



Interaction particle-particle 



particle-hole 





JT=01 


JT=10 


At=20 


Ar=40 


Ar=ll 


KB 


-4.75 


-4.46 


-2.79 


-1.39 


+2A6 


FPD6 


-5.06 


-5.08 


-3.11 


-1.67 


-h3.17 


GOGNY 


-4.07 


-5.74 


-3.23 


-1.77 


+2A6 


BonnC 


-4.20 


-5.60 


-3.33 


-1.29 


+2.70 



III. THE SOLUTION OF THE SECULAR PROBLEM IN A 
FINITE SPACE 

Once the interaction and the valence space ready, it is 
time to construct and diagonalize the many body secular 
matrix of the Hamiltonian describing the (effective) in- 
teraction between valence particles in the valence space. 
Two questions need now consideration: which basis to 
take to calculate the non-zero many-body matrix ele- 
ments (NZME) and which method to use for the diag- 
onalization of the matrix. 



A. The Lanczos method 



10 12 14 16 



FIG. 8 m-scheme dimensions and total number of non-zero 
matrix elements in the p/-shell for nuclei with M = Tz = 

In the standard diagonalization methods (Wilkinson, 
1965) the CPU time increases as N^, N being the di- 
mension of the matrix. Therefore, they cannot be used 
in large scale shell model (SM) calculations. Nuclear SM 
calculations have two specific features. The first one is 
that, in the vast majority of the cases, only a few (and 



very often only one) eigenstates of a given angular mo- 
mentum (J) and isospin (T) arc needed. Secondly, the 
matrices arc very sparse. As can be seen in figure 8, the 
number of NZME varies linearly instead of quadratically 
with the size of the matrices. For these reasons, itera- 
tive methods are of general use, in particular the Lanczos 
method (Lanczos, 1950). As an alternative, the Davidson 
method (Davidson, 1975) has the advantage of avoiding 
the storage of a large number of vectors, however, the 
large increase in storage capacity of modern computers 
has somehow minimized these advantages. 

The Lanczos method consists in the construction of an 
orthogonal basis in which the Hamiltonian matrix (H) is 
tridiagonal. 

A normalized starting vector (the pivot) |1) is chosen. 
The vector |ai) = H\l) has necessarily the form 



\ai) 



iJnll) + 12'), with (112') = 0. 



(45) 



Calculate Hu = {1\H\1) through (l|ai) = Hn. Nor- 
malize |2) = |2')/(2'|2')i/2 to find H12 = {1\H\2) = 
(2'|2')^/2 Iterate until state |A;) has been found. The 
vector |ofe) = H\k) has necessarily the form 



\ak) = Hkk-i\k - 1) + Hkk\k) + \(k + 1)'). 



(46) 



Calculate (fcjofc) = Hkk- Now \{k + 1)') is known. Nor- 
malize it to find Hkk+i = ((fc + l)'|(fc -}- l)')i/2. 

The (real symmetric) matrix is diagonalized at each 
iteration and the iterative process continues until all the 
required eigenvalues are converged according to some cri- 
teria. The number of iterations depends little on the di- 
mension of the matrix. Besides, the computing time is 
directly proportional to the number of NZME and for 
this reason it is nearly linear (instead of cubic as in the 
standard methods) in the dimension of the matrix. It de- 
pends on the number of iterations, which in turn depends 
on the number of converged states needed, and also on 
the choice of the starting vector. 

The Lanczos method can be used also as a projection 
method. For example in the m-scheme basis, a Lanc- 
zos calculation with the operator will generate states 
with well defined J. Taking these states as starting vec- 
tors, the Lanczos procedure (with H) will remain inside 
a fixed J subspace and therefore will improve the con- 
vergence properties. When only one converged state is 
required, it is convenient to use as pivot state the solu- 
tion obtained in a previous truncated calculation. For 
instance, starting with a random pivot the calculation of 
the ground state of ^°Cr in the full p f space (dimension 
in m-scheme 14,625,540 Slater determinants) needs twice 
more iterations than if we start with the pivot obtained 
as solution in a model space in which only 4 particles 
are allowed outside the I/7/2 shell (dimension 1,856,720). 
The overlap between these two 0+ states is 0.985. When 
more eigenstates arc needed (ric > 1), the best choice for 
the pivot is a linear combination of the ric lower states 
in the truncated space. Even if the Lanczos method is 
very efficient in the SM framework, there may be nu- 
merical problems. Mathematically the Lanczos vectors 
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are orthogonal, however numerically this is not strictly 
true due to the limited floating point machine precision. 
Hence, small numerical precision errors can, after many 
iterations, produce catastrophes, in particular, the states 
of lowest energy may reappear many times, rendering the 
method inefficient. To solve this problem it is necessary 
to orthogonalize each new Lanczos vector to all the pre- 
ceding ones. The same precision defects can produce the 
appearance of non expected states. For example in a m- 
scheme calculation with a J = 4, M = pivot, when 
many iteration are performed, it may happen that one 
J = and even one J = 2 state, lower in energy than 
the J = 4 states, show up abruptly. This specific prob- 
lem can be solved by projecting on J each new Lanczos 
vector. 



B. The choice of the basis 

Given a valence space, the optimal choice of the basis 
is related to the physics of the particular problem to be 
solved. As we discuss later, depending on what states or 
properties we want to describe (ground state, yrast band, 
strength function,. . . ) and depending on the type of nu- 
cleus (deformed, spherical,. . . ) different choices of the 
basis arc favored. There are essentially three possibilities 
depending on the underlying symmetries: m-scheme, J 
coupled scheme, and JT coupled scheme. 



TABLE III Some dimensions in the pf shell 



A 


4 


8 


12 


16 


20 


M 


= T,= 


-- 


4000 


2 X 10'^ 


1.10 X 10** 


1.09 X lO'* 


2.29 X 10'-* 


J 







156 


41355 


1.78 X 10^ 


1.54 X 10'' 


3.13 X 10^ 


J 


= T = 





66 


9741 


3.32 X 10^ 


2.58 X 10® 


5.05 X 10® 



As the m-scheme basis consists of Slater Determinants 
(SD) the calculation of the NZME is trivial, since they are 
equal to the decoupled two body matrix element (TBME) 
up to a phase. This means that, independently of the size 
of the matrix, the number of possible values of NZME is 
relatively limited. However, the counterpart of the sim- 
plicity of the m-scheme is that only and are good 
quantum numbers, therefore all the possible (J, T) states 
are contained in the basis and as a consequence the di- 
mensions of the matrices are maximal. For a given num- 
ber of valence neutrons, n,,, and protons, Zy the number 
of different Slater determinants that can be built in the 
valence space is : 



d = 




(47) 



where D„ and Dp are the degeneracies of the neutron and 
proton valence spaces. Working at fixed M and T^ the 

bases are smaller (d = X^mt {M,Tz)). The J or JT 
coupled bases, split the full m-scheme matrix in boxes 



whose dimensions are much smaller. This is especially 
spectacular for the J = T = states (sec tabic III). 

It is often convenient to truncate the space. In the par- 
ticular case of the pf shell, calling / the largest subshell 
(/7/2)) and r, generically, any or all of the other subshells 
of the p = 3 shell, the possible f-truncations involve the 
spaces 



■+f 



m—mo—tmo+t 



(48) 



where mo ^ if more than 8 neutrons (or protons) are 
present. For t = m — mo we have the full space {pf)"^ 
for ^ = 40 -I- m. 

In the late 60's, the Rochester group developed the 
algorithms needed for an efficient work in the (J, T) 
coupled basis and implemented them in the Oak-Ridge 
Rochester Multi-Shell code (French et al, 1969). The 
structure of the calculation is as follows: In the first 
place, the states of rii particles in a given ji shell are 
defined: |7i) = \{ji)^^ViJiXi) [vi is the seniority and Xi 
any extra quantum number). Next, the states of N par- 
ticles distributed in several shells are obtained by succes- 
sive angular momentum couplings of the one-shell basic 
states: 



[|7i)l72)r l73) ...\lk) 



(49) 



Compared to the simplicity of the m-scheme, the cal- 
culation of the NZME is much more complicated. It in- 
volves products of 9j symbols and coefficients of fractional 
parentage (cfp), i.e., the single-shell reduced matrix ele- 
ments of operators of the form 



i^n^n)\ 4' (50) 

This complexity explains why in the OXBASH 
code (Brown et al., 1985), the JT-coupled basis states 
are written in the m-scheme basis to calculate the 
NZME. The Oxford-Buenos Aires shell model code, 
widely distributed and used has proven to be an invalu- 
able tool in many calculations. Another recent code that 
works in JT-coupled formalism is the Drexel University 
DUPSM (Novoselsky and Vallieres, 1997) has had a much 
lesser use. 

In the case of only J (without T) coupling, a strong 

simplification in the calculation of the NZME can be 
achieved using the Quasi-Spin formalism (Ichimura and 
Arima, 1966; Kawarada and Arima, 1964; Lawson and 
Macfarlane, 1965), as in the code NATHAN described 
below. The advantages of the coupled scheme decrease 
also when J and T increase. As an example, in ^^Ni the 
ratio dim(M = J)/dim(J) is 70 for J = but only 5.7 
for J = 6. The coupled scheme has another disadvantage 
compared to the m-schcmc. It concerns the percentage 
of NZME. Let us give the example of the state 4+ in ^°Ti 
(full pf space); the percentages of NZME are respectively 
14% in the JT-basis (see Novoselsky et al, 1997), 5% in 
the J-basis and only 0.05% in m-scheme. For all these 
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reasons and with the present computing facilities, we can 
conclude that the m-scheme is the most efficient choice 
for usual SM calculations, albeit with some notable ex- 
ceptions that we shall mention below. 



C. The Glasgow m-scheme code 

The steady and rapid increase of computer power in 
recent years has resulted in a dramatic increase in the 
dimensionality of the SM calculations. A crucial point 
nowadays is to know what are the limits of a given com- 
puter code, their origin and their evolution. As far as the 
NZME can be calculated and stored, the diagonalization 
with the Lanczos method is trivial. It means that the 
fundamental limitation of standard shell model calcula- 
tions is the capacity to store the NZME. This is the origin 
of the term "giant" matrices, that we apply to those for 
which it is necessary to recalculate the NZME during the 
diagonalization process. The first breakthrough in the 
treatment of giant matrices was the SM code developed 
by the Glasgow group (Whitehead et al., 1977). Let us 
recall its basic ideas: It works in the m-scheme and each 
SD is represented in the computer by an integer word. 
Each bit of the word is associated to a given individual 
state \nljmT). Each bit has the value 1 or depending 
on whether the state is occupied or empty. A two-body 
operator ajajafea; will select the words having the bits 
i,j,k,l in the configuration 0011, say, and change them 
to 1100, generating a new word which has to be located 
in the list of all the words using the bi-section method. 



D. The m-scheme code ANTOINE 

The shell model code ANTOINE^s (Caurier and 
Nowacki, 1999) has retained many of these ideas, while 
improving upon the Glasgow code in the several aspects. 
To start with, it takes advantage of the fact that the di- 
mension of the proton and neutron spaces is small com- 
pared with the full space dimension, with the obvious ex- 
ception of semi-magic nuclei. For example, the 1,963,461 
SD with AI = in '^^Cr are generated with only the 4,865 
SD (corresponding to all the possible M values) in '*'*Ca. 
The states of the basis are written as the product of two 
SD, one for protons and one for neutrons: |/) = |«,a). 
We use, /, J, capital letters for states in the full space, 
i,j, small case Latin letters for states of the first subspace 
(protons or neutrons), a, /3, small case Greek letters for 
states of the second subspace (neutrons or protons). The 
Slater determinants i and a can be classified by their M 
values. Ml and M2. The total M being fixed, the SD's of 



A version of the code can be downloaded from the URL http: 
//sbgatl94 . in2p3 . fr/'theory/antoine /main. html 



the 2 subspaces will be associated only if Mi -I- M2 = M. 
A pictorial example is given in fig. 9. 




FIG. 9 Schematic representation of the basis 

It is clear that for each \i) state the allowed \a) states 
run, without discontinuity, between a minimum and a 
maximum value. Building the basis in the total space by 
means of an i-loop and a nested a-loop, it is possible to 
construct numerically an array R{i) that points to the 
|/) state^^: 

I = R{i) + a (51) 

For example, according to figure 9, the numerical values 
of R are: i?(l) = 0, R{2) = 6, R{3) = 12, ...The re- 
lation (51) holds even in the case of truncated spaces, 
provided we define sub-blocks labeled with M and t (the 
truncation index defined in (48)). Before the diagonal- 
ization, all the calculations that involve only the pro- 
ton or the neutron spaces separately are carried out and 
the results stored. For the proton-proton and neutron- 
neutron NZME the numerical values of R{i), R{j),Wij 
and a, /?, Wap, where {i\H\j) = Wij and {a\H\l3) = 
are pre-calculated and stored. Therefore, in the Lanc- 
zos procedure a simple loop on a and i generates all 
the proton-proton and neutron-neutron NZME, Wjj — 
{I\H\J). For the proton-neutron matrix elements the 
situation is slightly more complicated. Let's assume that 
the \i) and \ j) SD are connected by the one-body operator 
a^gUr (that in the list of all possible one-body operators 
appears at position p), with q — nljm and r = n'l'j'm' 
and m! — m = Am. Equivalently, the \a) and |/3) SD 
are connected by a one-body operator whose position is 
denoted by ^. We pre-calculate the numerical values of 
R{i),R{j),p and a, (3, 11. Conservation of the total M 
implies that the proton operators with Am must be as- 
sociated to the neutron operators with —Am. Thus we 
could draw the equivalent to figure 9 for the proton and 
neutron one-body operators. In the same way as we did 
before for / = R(i) + a, we can now define an index 
K = Q{p) + n that labels the different two-body matrix 



We use Dirac's notation for the quantum mechanical state 
the position of this state in the basis is denoted by i. 
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elements (TBME). Then, wc denote V{K) the numeri- 
cal value of the proton-neutron TBME that connects the 
states \i,a) and Once R{i), R{j),Q{p) and a,P,^ 

are known, the non-zero elements of the matrix in the 
full space are generated with three integer additions: 

I = R(i) + a, J = RU)+f3, K = Q{p)+n. (52) 

The NZME of the Hamiltonian between states |7) and 
I J) is then: 

{I\H\J) = {J\H\I) = V{K). (53) 

The performance of the code is optimal when the 
two subspaces have comparable dimensions. It becomes 
less efficient for asymmetric nuclei (for semi-magic nuclei 
all the NZME must be stored) and for large truncated 
spaces. No-core calculations are typical in this respect. 
If wc consider a. N ~ Z nucleus like ^Li, in a valence 
space that comprises up to 14fia; configurations, there 
are 50000 Slater determinants for protons and neutrons 
with M = 1/2 and t = 14. Their respective counterparts 
can only have M — —1/2 and t — and have dimen- 
sionality 1. This situation is the same as in semi-magic 
nuclei. 

As far as two Lanczos vectors can be stored in the RAM 
memory of the computer, the calciilations arc straightfor- 
ward. Until recently this was the fundamental limitation 
of the SM code ANTOINE. It is now possible to overcome 
it dividing the Lanczos vectors in segments: 



(54) 



The Hamiltonian matrix is also divided in blocks so that 
the action of the Hamiltonian during a Lanczos iteration 
can be expressed as: 



(fc) 



(55) 



The k segments correspond to specific values of M (and 
occasionally t) of the first subspace. The price to pay for 
the increase in size is a strong reduction of performance of 
the code. Now. {I\H\J) and {J\H\I) are not generated 
simultaneously when |/) and |J) do not belong at the 
same k segment of the vector and the amount of disk use 
increases. As a counterpart, it gives a natural way to par- 
allelize the code, each processor calculating some specific 
'J' J . This technique allows the diagonalization of matri- 
ces with dimensions in the billion range. All the nuclei of 
the j3/-shell can now be calculated without truncation. 
Other m-scheme codes have recently joined ANTOINE 
in the run, incorporating some of its algorithmic findings 
as MSHELL (Mizusaki, 2000) or REDSTICK (Ormand 
and Johnson, 2002). 



E. The coupled code NATHAN. 

As the mass of the nucleus increases, the possibilities of 

performing SM calculations become more and more lim- 
ited. To give an order of magnitude the dimension of the 



matrix of '''^Fe with 6 protons and 6 neutrons in the pf 
shell is 10®. If now we consider ^'^^Ba with also 6 protons 
and 6 ncutron-holcs in a valence space of 5 shells, four 
of them equivalent to the p f shell plus the lft.11/2 orbit, 
the dimension reaches 2 x 10^". Furthermore, when many 
protons and neutrons are active, nuclei tend to get de- 
formed and their description requires even larger valence 
spaces. For instance, to treat the deformed nuclei in the 
vicinity of ^''Zr, the "normal" valence space, 2p^/2-, I/5/2) 
2pi/2) I59/2) must be complemented with, at least, the 
orbit 1(^5/2. This means that our span will be limited to 
nuclei with few particles outside closed shells (^''"Xe re- 
mains an easy calculation) or to almost semi-magic, as for 
instance the Tellurium isotopes. These nuclei are spheri- 
cal, therefore the seniority can provide a good truncation 
scheme. This explains the interest of a shell model code 
in a coupled basis and quasi-spin formalism. 

In the shell model code NATHAN (Caurier and 
Nowacki, 1999) the fundamental idea of the code AN- 
TOINE is kept, i. e. splitting the valence space in two 
parts and writing the full space basis as the product of 
states belonging to these two parts. Now, |z) and la) are 
states with good angular momentum. They are built with 
the usual techniques of the Oak- Ridge/Rochester group 
(French et al, 1969). Each subspace is now partitioned 
with the labels Ji and J2. The only difference with the 
m-scheme is that instead of having a one-to-one associa- 
tion (Ml -|- M2 = M), for a given Ji we now have all the 

possible J2, Jmin < J2 < Jmax, with J™„ = \Jq - Ji| 

and Jmax = Jq + Ji- The continuity between the first 
state with Jmin and the last with Jmax is maintained and 
consequently the fundamental relation / = R{i) -\- a still 
holds. The generation of the proton-proton and neutron- 
neutron NZME proceeds exactly as in m-scheme. For the 
proton-neutron NZME the one-body operators in each 
space can be written as = (aj^a^j)'^. There exists 
a strict analogy between Am in m-scheme and A in the 
coupled scheme. Hence, we can still establish a relation 
K = Q{p) + uj. The NZME read now: 



{I\H\J) = {J\H\I) = ■ Kp ■ W{K), 
with hij = {i\0^\j), /la/3 = {a\0^\l3), and 



W{K) DC V{K) 




(56) 



(57) 



where V{K) is a TBME. We need to perform -as in the 
m-scheme code- the three integer additions which gener- 
ate I, J, and K, but, in addition, there are two floating 
point multiplications to be done, since hij and hafj, which 
in m-scheme were just phases, are now a product of cfp's 
and 9j symbols (see formula 3.10 in French et al. (1969)). 
Within this formalism we can introduce seniority trun- 
cations, but the problem of the semi-magic nuclei and 
the asymmetry between the protons and neutrons spaces 
remains. To overcome this difficulty in equalizing the 
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dimension of the two subspaccs wc have generalized the 
code as to ahow that one of the two subspaces contains 
a mixing of proton and neutron orbits. For example, for 
heavy Scandium or Titanium isotopes, we can include 
the neutron I/7/2 orbit in the proton space. For the 
N = 126 isotones we have only proton shells. We then 
take the lii3/2 and I/19/2 orbits as the first subspace, 
while the 2/7/2, 2/5/2, 8^3/2 and 3pi/2 form the second 
one. The states previously defined with Jp and J„ are 
now respectively labelled ^1, Ji and A2, J2, Ai and A2 
being the number of particles in each subspace. In hij 
and hap now appear mean values of all the operators in 
Eq.(50) 

The fact that in the coupled scheme the dimensions 
are smaller and that angular momentum is explicitly en- 
forced makes it possible to perform a very large number 
of Lanczos iterations without storage or precision prob- 
lems. This can be essential for the calculation of strength 
functions or when many converged states are needed, for 
example to describe non-yrast deformed bands. The cou- 
pled code has another important advantage. It can be 
perfectly parallelized without restriction on the number 
of processors. The typical working dimension for large 
matrices in NATHAN is 10^ (10^ with ANTOINE). It 
means that there is no problem to define as many fi- 
nal vectors as processors are available. The calculation 
of the NZME is shared between the different proces- 
sors (each processor taking a piece of the Hamiltonian, 
H = iJ^'^)) leading to different vectors that are added 
to obtain the full one: 

= H^k)^,^ */ = ^ *f ^- (58) 

k 

F. No core shell model 

The ab initio no-core shell model (NCSM) (Navratil 
et al., 2000a, b) is a method to solve the nuclear many 
body problem for light nuclei using realistic inter-nucleon 
forces. The calculations are performed using a large but 
finite harmonic-oscillator (HO) basis. Due to the basis 
truncation, it is necessary to derive an effective interac- 
tion from the underlying inter-nucleon interaction that 
is appropriate for the basis size employed. The effective 
interaction contains, in general, up to ^-body compo- 
nents even if the underlying interaction had, e.g. only 
two-body terms. In practice, the effective interaction is 
derived in a sub-cluster approximation retaining just two- 
or three-body terms. A crucial feature of the method is 
that it converges to the exact solution when the basis 
size increases and/or the effective interaction clustering 
increases. 

In a first phase, their applications were limited to the 
use of realistic two-nuclcon interactions, either G- matrix- 
based two-body interactions (Zheng et al., 1994), or de- 
rived by the Lee-Suzuki procedure (Suzuki and Lee, 1980) 
for the NCSM (Navratil and Barrett, 1996). This re- 
sulted in the elimination of the purely phenomenolog- 



ical parameter A used to define G-matrix starting en- 
ergy. A truly ab initio formulation of the formalism was 
presented by Navratil and Barrett (1998), where conver- 
gence to the exact solutions was demonstrated for the 
A = 3 system. The same was later accomplished for the 
A = 4 system (Navratil and Barrett, 1999), where it was 
also shown that a three-body effective interaction can be 
introduced to improve the convergence of the method. 
The capability to derive a three-body effective interac- 
tion and apply it in either relative-coordinate (Navratil 
et al., 2000) or Cartesian-coordinate formalism (Navratil 
and Ormand, 2002) together with the ability to solve a 
three-nucleon system with a genuine three-nuclcon force 
in the NCSM approach (Marsden et al., 2002) now opens 
the possibility to include a realistic three-nucleon force 
in the NCSM Hamiltonian and perform calculations us- 
ing two- and three-nucleon forces for the p-shell nuclei. 
Among successes of the NCSM approach was the first 
published result of the binding energy of *He with the 
CD-Bonn NN potential (Navratil and Barrett, 1999), the 
near-converged results for A = 6 using a non-local Hamil- 
tonian (Navratil et al., 2001), the first observation of the 
incorrect ground- state spin in "'^"B predicted by the realis- 
tic two-body nucleon-nucleon interactions (Cauricr et al., 
2002). This last result, together with the already known 
problems of under-binding, confirms the need of realistic 
three-nucleon forces, a conclusion that also stems from 
the Green Function Monte Carlo calculations of Wiringa 
and Pieper (2002) (see also Pieper and Wiringa, 2001, 
for a review of the results of the Urbana-Argonne collab- 
oration for nuclei A<10). 

Recently, a new version of the shell model code AN- 
TOINE has been developed for the NCSM applications 
(Caurier et ai, 2001b). This new code allows to perform 
calculations in significantly larger basis spaces and makes 
it possible to reach full convergence for the A = 6 nuclei. 
Besides, it opens the possibility to investigate slowly con- 
verging intruder states and states with unnatural parity. 
The largest bases reached with this code so far are, ac- 
cording to the number of HO excitations, the Uhw space 
for ^Li and, according the matrix dimension, the lOfko 
calculations for ^"C. In the latter case, the m-scheme 
matrix dimension exceeds 800 millions. 

G. Present possibilities 

The combination of advances in computer technology 
and in algorithms has enlarged the scope of possible SM 
studies. The rotational band of ^^sy 3^(3^18 to be still very 
far from the reach of SM calculations, but predictions for 
^^^U are already available. Let us list some of the present 
opportunities. 

i) No-corc calculations. One of the major problems 
of the no-core SM is the convergence of the results 
with the size of the valence space; For ^Li we can 
handle excitations until 14fuj and at least SfiaJ for 
all the p-shell nuclei. 
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ii) The p/-shell. This where the shell model has been 
most successful and exact diagonalizations are now 
possible throughout the region Beyond ^^Ni, as the 
1/7/2 orbit becomes more and more bound, trun- 
cated calculations are close to exact, for instance, 
in ^°Zn (Mazzocchi et ai, 2001) the wave-functions 
are fully converged when 6p-6h excitations are in- 
cluded. 

iii) The g valence space. We use the notation Vp 
for the set of the nlj orbits with 2{n — 1) + I = p 
excluding the orbit with maximum total angular 
momentum j — p + 1/2. This space describes nu- 
clei in the region 28 < N, Z < 50; ^^Ni is taken as 
inert core. Most of them are nearly spherical, and 
can be treated without truncations. The f3/3 decay 
(with and without neutrinos) of ^^Ge and ^^Se are a 
prime example. The deformed nuclei [N Z ~ 40) 
are more difficult because they demand the inclu- 
sion of the 2d^/2 orbit to describe prolate states and 
oblate-prolate shape coexistence. 

iv) The pfg space. To enlarge the pf valence space 
to include the 1^9/2 orbit is hopeless nowadays. 
Furthermore, serious center of mass spuriousness is 
expected in the I/7/2I59/2 configurations (see Ap- 
pendix C). For neutron rich nuclei in the Nickel 
region, a tractable valence space that avoids the 
center of mass problem can be defined: A "^^Ca 
core on top of which pf protons and g neutrons 
are active. 

v) Heavy nuclei. All the semi-magic nuclei, for in- 
stance the = 126 isotones, can be easily studied 
and the addition of few particles or holes remains 
tractable. Some long chains of Tellurium and Bis- 
muth isotopes have been recently studied (Caurier 
et ai, 2003). 



IV. THE LANCZOS BASIS 

There exists a strong connection between the 
Lanczos algorithm, the partition function, Z(f3) = 
exp(— /3H)|i), and the evolution operator exp{iHt). 
In the three cases, powers of the Hamiltonian determine 
the properties of the systems. The partition function 
can be written as Z{P) — p(i?) exp(— /3i?) i.e., the 
Laplace transform of the density of states, a quantity 
readily accessible once the Hamiltonian has been fully 
diagonalized. The evolution operator addresses more 
specifically the problem of evolving from a starting vec- 
tor into the exact ground state. We shall discuss both 
questions in turn. 



A. Level densities 

For many, shell model is still synonymous with diag- 
onalizations, in turn synonymous with black box. One 



still hears questions such as: Who is interested in di- 
agonalizing a large matrix? As an answer we propose to 
examine Fig. 10 showing the two point-functions that de- 
fine the diagonal. Ha, and off-diagonal elements, 
in a typical Lanczos construction. The continuous lines 




i/dim 

FIG. 10 Tridiagonal matrix elements for a 6579 dimensional 
matrix, and the logarithmic and inverse binomial approxima- 
tions (nib) (Zuker et ai, 2001). 

are calculated knowing the first four moments of the ma- 
trix {i.e. of the l-|-2-body Hamiltonian H it represents). 
These results hold at fixed quantum numbers, i.e., when 
the matrix admits no block decomposition. 

When the matrix is diagonalized the level density is 
well reproduced by a continuous binomial 



Pb{x,N,p, S) 



r{N + i) 



N 



r{xN +l)T{xN + 1) S 



(59) 



where x is an a-dimensional energy, x — 1 — x, N the 
number of valence particles, p an asymmetry parameter, 
p + q ~ 1, and S the span of the spectrum (distance 
between lowest and highest eigenstates). Introducing the 
energy scale e, S, the centroid Ec, the variance and x 
are given by 



S = Ne, Ec = Npe, = Np{\ ~ p)e^ 



f, (60, 



Note that pi,{x, N,p, S) reduces to a discrete binomial, 
(^) if x = n/N — ne/ S, with integer n. 

N , p and S are calculated using the moments of the 
Hamiltonian TL, i.e., averages given by the traces of 
Ti.^ , to be equated to the corresponding moments of 
Pb{x, N,p, S), which for low K are the same as those of a 
discrete binomial. The necessary definitions and equali- 
ties follow. 



d-hY{n'^) = (7^^), E, - {n^), Mk = {{n - e,) 
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These quantities also define the logarithmic and inverse 
binomial (nib) forms of Ha and Ha^i in Fig. 10. Note 
that the corresponding lines are almost impossible to dis- 
tinguish from those of the exact matrix. The associated 
level densities are found in Fig. 11. 




2 4 6 8 10 12 14 16 

E(IVIeV) 

FIG. 11 Exact and binomial level densities for the matrix in 
Fig. 10. 

The mathematical status of these results is somewhat 
mixed. Mon and French (1975) proved that the total 
density of a Hamiltonian system is to a first approxima- 
tion a Gaussian. Zuker (2001) extended the approxima- 
tion to a binomial, but the result remains valid only in 
some neighborhood of the centroid. Furthermore, one 
does not expect it to hold generally because binomial 
thermodynamics is trivial and precludes the existence of 
phase transitions. In the example given above, we do 
not deal with the total density, which involves all states, 
p = jrp{'2,J + 1)(2T -f l)/5,/T but with a partial pjT at 
fixed quantum numbers. In this case Zuker et al. (2001) 
conjectured that the tridiagonal elements given by the 
logarithmic and nib forms are valid, and hence describe 
the full spectrum. The conjecture breaks down if a dy- 
namical symmetry is so strong as to define new (approx- 
imately) conserved quantum numbers. 

Granted that a single binomial cannot cover all situa- 
tions we may nonetheless explore its validity in nuclear 
physics, where the observed level densities are extremely 
well approximated by the classical formula of Bethe 
(1936) with a shift A, 

PB{E,a,A) = -^ (4^)V4(j, + A)5/4 - (62) 

Obviously, if binomials are to be useful, they must 
reproduce — for some range of energies — Eq. (62). They 
do indeed, as shown in Fig. 12, where the experimen- 
tal points for ^''Ni (Iljinov et al., 1992) are also given. 
Though Bethe and binomial forms are seen to be equiv- 
alent, the latter has the advantage that the necessary 
parameters are well defined and can be calculated, while 
in the former, the precise meaning of the a parameter 
is elusive. The shift A is necessary to adjust the ground 
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FIG. 12 Binomial, Bethe and experimental level densities for 
«°Ni. 



state position. The problem also exists for the binomial^* 
and the result in the figure (Zuker, 2001) solves it phe- 
nomenologically. The Shell Model Monte Carlo method 
provides the only parameter-free approach to level den- 
sities (Dean et ai, 1995; Langanke, 1998; Nakada and 
Alhassid, 1997) whose reliability is now established (Al- 
hassid et ai, 1999). The problem is that the calculations 
are hard. 

As the shape of the level density is well reproduced 
by a binomial except in the neighborhood of the ground 
state, to reconcile simplicity with full rigor we have to 
examine the tridiagonal matrix at the origin. 



B. The ground state 

Obviously, some dependence on the pivot should exist. 
We examine it through the J = IT — 3 pf states in 
''^Sc^^. The matrix is 8590-dimensional, and we calcu- 
late the ground state with two pivots, one random (ho- 
mogeneous sum of all basic states) , the other variational 
(lowest eigenstate in space). A zoom on the first ma- 
trix elements in Fig. 13 reveals that they are very differ- 
ent for the first few iterations, but soon they merge into 
the canonical patterns discussed in the preceding section. 
The ground state wavefunction is unique of course, but 
it takes the different aspects shown in Fig. 14. In both 
cases the convergence is very fast, and it is not difficult to 
show in general that it occurs for a number of iterations 
of order NlogN for dimensionality d = 2^ . However, 
the variational pivot is clearly better if we are interested 
in the ground state: If its overlap with the exact solu- 
tion exceeds 50%, all other contributions are bound to 



The energies are referred to the mean (centroid) of the distribu- 
tion, while they are measured with respect to the ground state. 
They are reached in the **Ca(p, n)(48)Sc reaction, which will be 
our standard example of Gamow Teller transitions. 
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FIG. 13 Tridiagonal matrix elements after 100 iterations for 
a random and a variational pivot. 



be smaller and in general they will decrease uniformly^*^ . 
We shall see how to exploit this property of good pivots 
to simplify the calculations. 



1. The exp(S) method 
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FIG. 14 The ground state wavefunction in the Lanczos basis 
for a random, and a variational pivot. The plotted values are 
the squared overlaps of the successive Lanczos vectors with 
the final wave-function 



In Section II. A. 1 we have sketched the coupled-cluster 
(or exp(5')) formalism. The formulation in the Lanczos 
basis cannot do justice to the general theory, but it is a 
good introduction to the underlying ideas. Furthermore, 
it turns out to be quite useful. 

The construction is as in Eqs. (45,46), but the succeed- 
ing vectors — except the pivot — are not normalized. Then, 



There are some very interesting counter-examples. One is found 
in Fig. 24, where the natural pivot is heavily fragmented. 



the full wavefunction takes the form 

|0) = (1 + ciPi + C2P2 + • • • + c/P/)|0), (63) 

where Pm is a polynomial in H that acting on the pivot 
produces orthogonal unnormalized vectors in the Lanczos 
basis: Pm\0) — Eq. (46) becomes 



n\m)=V^\m-l) 



I) 



(64) 



To relate with the normalized version (m m); Divide 
by (m|m)"'^/^, then multiply and divide |m— 1) and |m-|-l) 
by their norms to recover Eq. (46), and obtain Em — 



Hr, 



Vm 



{m\m) / {m — l\m — 1) 



The 



secular equation (H — E)\0) 



tt2 

m m — l 

leads to the recursion 



-1 + (Em — E) Cm + Vm+l Cm+1 = 0, 



(65) 



whose solution is equivalent to diagonalizing a matrix 
with Vm+i in the upper diagonal and 1 in the lower one, 
which is of course equivalent to the symmetric problem. 
However, here we shall solve the recursion by transform- 
ing it into a set of non-linear coupled equations for the 
Cm amplitudes. The profound reason for using an unnor- 
malized basis is that first term in Eq. (65), E — E^^ViCi 
implies that once ci is known the problem is solved. Call- 
ing Em = Em — Eo and replacing E ~ Eq+ViCi in Eq (65) 
leads to 



VlCi) Cm + Vm+l Cm+1 = 0. 



(66) 



Cm-l + (Sr, 

Now introduce 

CmPm = exp(y^ S in Pin ), (67) 
Expanding the exponential and equating terms 



C3 = S'i 



C2 = 52 + ^5*1, 



5152 + ■^'S'l, 



etc. 



(68) 
(69) 



Note that we have used the formal identification Pm Pn = 
Pn+m as a heuristic way of suggesting Eqs. (68). In- 
serting in Eq. (66) and regrouping, a system of coupled 
equations obtains. The first two are 



=S2V2 + Sf{^V2 - Vi) + SiEi + 1 



=53l4i + 5222 + 5i 52(1^3 - 

lsfe2 + S!il^V,-lvi) 



Vi)+ 
+ Si. 



(70) 



(71) 



Hence, instead of the usual Cm-truncations we can use 
5m-truncations, which have the advantage of providing a 
model for the wavefunction over the full basis. At the first 
iteration Eq. (70) is solved neglecting 52. The resulting 
value for 5i is inserted in Eq. (71) which is solved by 
neglecting 53. The resulting value for 52 is reinserted 
in Eq. (70) and the process is repeated. Once 5i and 
52 are known, the equation for 53 (not shown) may be 
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FIG. 15 Energy convergence for ground and first excited state 
in ^®Ni as a function of the truncation level t (Caurier et al., 
1999c). 



incorporated, and so on. Since the energy depends only 
on Si, convergence is reached when its value remains 
constant from step to step. With a very good pivot, 
Eq. (70) should give a fair approximation, improved by 
incorporating Eq. (71), and checked by the equation for 

If all amplitudes except Si are neglected, the differ- 
ence scheme (65) is the same that would be obtained for 
an harmonic 7i = + ViS'^ + Successive ap- 

proximations amount to introduce anharmonicities. An 
equivalent approach — numerically expedient — consists in 
diagonalizing the matrices at each c„i truncation level. 
Then, after some iteration, the energies converge expo- 
nentially: Introduce 



conv(«,a, eo,«o) = eo 



exp(— ai) — 1 
exp(— aio) — 1 ' 



(72) 



which equals eo at point i = io- Choose a so as to yield 
conv(i = io + 1) = e(io + 1) ie., the correct energy 
at the next point . Check that e{io + 2) is well repro- 
duced. Fig. 15 provides an example of the efficiency of 
the method for the lowest two states in ^^Ni. The index 
i is replaced by the truncation level i = t/2 [Eq. (48, 
with m = 16, mo = 0]: We set io = 1, fix a so as to 
reproduce the energy at the second iteration and check 
that the curve gives indeed the correct value at the third 
point. The results reproduce those of the 5*2 truncation, 
confirming that "exponential convergence" and exp(S') 
are very much the same thing. In this example, the 
closed shell pivot is particularly good, and the exponen- 
tial regime sets in at the first iteration. In general, it will 
always set in for some value of ig that may be quite large 
(Fig. 14 suggests io « 25 for the random pivot). Fortu- 
nately, it is almost always possible to find good pivots, 
and the subject deserves a comment. 

In the case of ^^Ni the good pivot is a closed shell. As 
a consequence, the first iterations are associated to trun- 
cated spaces of much smaller dimensionalities than the 
total one {dm ~ 10^). This also happens for lighter pf 



shell nuclei, for which the fy2 [or eventually (.f7/2P3/2)'"] 
subspaces provide a good pivot. The same argument ap- 
plies for other regions. For well deformed nuclei Hartree 
Fock shoiild provide good determinantal pivots, and 
hence enormous gains in dimensionality, once projection 
to good angular momentum can be tackled efficiently. 

As we shall see next, the Lanczos and exp S procedures 
provide a convenient framework in which to analyze other 
approaches. 



2. Other numerical approximation methods 

The exponential convergence method introduced in the 
previous section was first described by Horoi et al. (1999), 
under a different but equivalent guise. In later work, a 
hierarchy of configurations determined by their average 
energy and width was proposed. Successive diagonaliza- 
tions make it possible to reach the exact energy by expo- 
nential extrapolation. The method has been successfully 
applied to the calculation of the binding energies of the 
pf shell nuclei (Horoi et al., 2002) and to calculation of 
the excitation energies of the deformed 0+ states of ^^Ni 
and 52 Cr (Horoi et al., 2003). 

The work of Mizusaki and Imada (2002, 2003), is based 
on the fact that the width of the total Hamiltonian in a 
truncated space tends to zero as the solution approaches 
the exact one. They have devised different extrapolation 
methods and applied them to some pf-sheW nuclei. 

Andreozzi, Lo ludice, and Porrino (2003) have recently 
proposed a factorization method that allows for impor- 
tance sampling to approximate the exact eigenvalues and 
transition matrix elements. 

A totally different approach which has proven its power 
in other fields, the density matrix renormalization group, 
has been proposed (Dukelsky and Pittel, 2001; Dukelsky 
et al., 2002), and Papcnbrock and Dean (2003) have de- 
veloped a method based in the optimization of product 
wavefunctions of protons and neutrons that seems very 
promising. 

Still, to our knowledge, the only prediction for a truly 
large matrix that has preceded the exact calculation re- 
mains that of '"'''Ni, J = 2 in Fig. 15. It was borne out 
when the diagonalization became feasible two years later. 



3. Monte Carlo methods 

Monte Carlo methods rely on the possibility to deal 
with the imaginary-time evolution operator acting on 
some trial wavcfunction exp(— /?7i)|0) which tends to the 
exact ground state |0) as /3 oo. This is very much what 
the Lanczos algorithm does, but no basis is constructed. 
Instead, the energy (or some observable Ci) is calculated 
through 



(Ole-z^^lO) 



(0|0) 



(73) 
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which is transformed into a quotient of multidimensional 
integrals evaluated through Monte Carlo methods with 
importance sampling. A "sign problem" arises because 
the integrands are not definite positive, leading to enor- 
mous precision problems. The Green Function Monte 
Carlo studies mentioned at the beginning of Section II are 
conducted in coordinate space. The Shell Model Monte 
Carlo (SMMC) variant is formulated in Fock space, and 
hence directly amenable to comparisons with standard 
SM results (see Koonin, Dean, and Langanke, 1997a, 
for a review). The approach relies on the Hubbard- 
Stratonovich transformation, and the sign problem is 
circumvented cither by an extrapolation method or by 
choosing Hamiltonians that are free of it while remain- 
ing quite realistic {e.g. pairing plus quadrupole). At 
present it remains the only approach that can deal with 
much larger valence spaces than the standard shell model. 
It should be understood that SMMC docs not lead to 
detailed spectroscopy, as it only produces ground state 
averages, but it is very well suited for finite tempera- 
ture calculations. The introduction of Monte Carlo tech- 
niques in the Lanczos construction is certainly a tempting 
project 

The quantum Monte Carlo diagonalization method 
(QMCD) of Otsuka, Honma, and Mizusaki (1998) con- 
sists in exploring the mean field structure of the valence 
space by means of Hartree-Fock calculations that break 
the symmetries of the Hamiltonian. Good quantum num- 
bers arc enforced by projection techniques (Peierls and 
Yoccoz, 1957). Then the authors borrow from SMMC to 
select an optimal set of basic states and the full Hamil- 
tonian is explicitly diagonalizcd in this basis. More 
basic states are itoratively added until convergence is 
achieved. For a very recent review of the applications 
of this method see (Otsuka et al., 2001b). A strong con- 
nection between mean field and shell model techniques is 
also at the heart of the Vampir approach (Petrovici et al, 
1999), (Schmid, 2001) and of the projected shell model 
(PSM) (Kara and Sun, 1995). 



C. Lanczos Strength Functions 

The choice of pivot in the Lanczos tridiagonal con- 
struction is arbitrary and it can be adapted to special 
problems. One of the most interesting is the calculation 
of strength functions (Bloom, 1984; Whitehead, 1980): If 
Uij is the unitary matrix that achieves diagonal form, its 
first column Ui q gives the amplitudes of the ground state 
wavcfunction in the tridiagonal basis while the first row 
Uoj determines the amplitude of the pivot in the j-th 
eigenstate. Uqj plotted against the eigenenergies Ej is 
called the "strength function" for that pivot. 

In practice, given a transition operator T, act with it 

on a target state \t) to define a pivot |0') = T\t) that 
exhausts the sum rule (0'|0') for T. Normalize it. Then, 
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FIG. 16 Evolution of the Gamow- Teller strength function of 
''^Ca as the number of Lanczos iterations on the doorway state 
increases. 
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by definition 

T\t) 

|0) = ^j^^n^ = Uoj\3), whose moments, (74) 
{Q\T-&\Q) = ^Ul^E^, are those of the (75) 

3 

Strength function S{E) = ^5{E - EjpQj, (76) 

j 

As the Lanczos vector |/) is obtained by orthogonahz- 



thc result for 50 iterations, and 250 kcV widths for the 
non converged states. The profiles become almost iden- 
tical. 




10 15 20 



Energy(MeV) 

FIG. 17 Strength functions convoluted with gaussians: Upper 
panel; 50 iterations. Bottom panel; 1000 iterations. 

ing H^O) to all previous Lanczos vectors \i), i < I, the 
tridiagonal matrix elements are linear combinations of 
the moments of the strength distribution. Therefore the 
eigensolutions of the 7x7 matrix define an approximate 
strength function Si{E) = Ez=i./ ^(^ - E,) {i\T\0)^ , 
whose first 27—1 moments are the exact ones. The 
eigenstates act as "doorways" whose strength will be split 
until they become exact solutions for 7 large enough, as 
illustrated in Fig. 16 that retraces the fragmentation pro- 
cess of the sum rule pivot; in this case a ^^Sc "doorway" 
obtained by applying the Gamow- Teller operator to the 
■^^Ca ground state. The term doorway applies to vectors 
that have a physical meaning but arc not eigenstates. Af- 
ter full convergence is achieved for all states in the reso- 
nance region, the strength function has the aspect shown 
at the bottom of the figure. In practice, all the spikes 
are affected by an experimental width, and assuming a 
perfect calculation, the observed profile would have the 
aspect shown at the bottom of Fig. 17, after convoluting 
with gaussians of 150 keV width. The upper panel shows 



V. THE OHm calculations 



In this section we first revisit the p, sd, and pf shells 
to explain how a 3b monopole mechanism solves prob- 
lems hitherto intractable. Then we propose a sample of 
pf shell results that will not be discussed elsewhere. Fi- 
nally Gamow Teller transitions and strength functions 
are examined in some detail. 



A. The monopole problem and the three-body interaction 

The determinant influence of the monopole interaction 

was first established through the mechanism of Bansal 
and French (1964). Its efficiency in cross-shell calcula- 
tions was further confirmed by Zamick (1965), and the 
success of the ZBM model (Zuker, Buck, and McGrory, 
1968) is implicitly due to a monopole correction to a real- 
istic force. Zuker (1969) identified the main shortcoming 
of the model as due to what must be now accepted as a 3b 
eff'ect^^. Trouble showed in 0?k<; calculations somewhat 
later simply because it takes larger matrices to detect it 
in the sd shell than in the ZBM space (dimensionalities 
of 600 against 100 for six particles). Up to 5 particles the 
results of Halbert et al. (1971) with a realistic interaction 
were quite good, but at ^^Na they were so bad that they 
became the standard example of the unreliability of the 
realistic forces (Brown and Wildenthal, 1988) and lead 
to the titanic^^ Universal SD (USD) fit of the 63 matrix 
elements in the shell by Wildenthal (1984). 

Though the pf shell demands much larger dimension- 
alities, it has the advantage of containing two doubly 
magic nuclei *^Ca and ^^Ni for which truncated calcula- 
tions proved sufficient to identify very early the core of 
the monopole problem: The failure to produce EI clo- 
sures. 



The ZBM model describes the region around through a 
Pi/2'51/2'^5/2 space. The French Bansal parameters bpd and bps 
(see Eq. (9b)) must change when going from ^^N to ^®0; which 
demands a 3b mechanism 
2^ It took two years on a Vax. Nowadays it would take an afternoon 
on a laptop. 
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FIG. 18 The spectrum of *^Ca for the KB, KB' and KB3 
interactions compared to experiment 

Fig. 18 gives an idea of what happens in "^^Ca with 
the KB interaction (Kuo and Brown, 1968): there are six 
states below 3 MeV, where only one exists. In Pasquini 
(1976) and Pasquini and Zuker (1978) the following mod- 
ifications were proposed (/ = f-j/2, r = P3/2, P1/2), 

y^^(KBl) = V^^iKB) - (-)^ 300 keV, 
V}0^(KB1) = Vj?^. (KB) - 350 keV, (77) 
y/^(KBl) = V//(KB) - 110 keV. 

The first line defines KB' in Fig. 18. The variants KB2 
and KB3 in Poves and Zuker (1981b) keep the KBl cen- 
troids and introduce very minor multipole modifications. 
KB3 was adopted as standard^^ in successful calculations 
in ^ = 47-50 that will be described in the next sec- 
tions (Caurier et ai, 1994; Martmez-Pinedo et ai, 1996b, 
1997). 

For higher masses, there are some problems, but noth- 
ing comparable to the serious ones encountered in the sd 
shell, where modifications such as in Eq. (77) are ben- 
eficial but (apparently) insufficient. The effective single 
particle energies in Fig. 19 are just the monopole values 
of the particle and hole states at the subshell closures. 
They give an idea of what happens in a 2b description. 
At the origin, in ^^Ca, the spectrum is the experimental 
one. At ^^Ni there is a bunching of the upper orbits that 
the realistic BonnC (BC) and KB describe reasonably 
well, but they fail to produce a substantial gap. Hence 
the need of KBl-type corrections. At the end of the shell 
BC and KB reproduce an expanded version of the "^^Ca 
spectrum, which is most certainly incorrect. The indi- 
cation from Tif^ in Eq. (11. B. 2) is that the bunching of 
the upper orbits should persist. By incorporating this 
hint — which involves the Vrr' centroids — and fine-tuning 



The multipole changes — that were beneficial in the perturbative 
treatment of Poves and Zuker (1981b) — had much less influence 
in the exact diagonalizations. 
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FIG. 19 Effective single particle energies in the p/-shell along 
the N = Z line, computed with the BC and KB3 interactions. 



the Vfr ones Poves et ai (2001) defined a KB3G inter- 
action that brings interesting improvements over KB3 in 
A = 50-52, but around ^^Ni there are still some problems, 
though not as severe as the ones encountered in the sd 
and p shells. Compared to KB3, the interaction FPD6 
has a better gap in ^^Ni, however, the orbit 115/2 is defi- 
nitely too low in '^^Ni. This produces problems with the 
description of Gamow- Teller processes (see Borcea et at, 
2001, for a recent experimental check in the beta decay 
of 56 Cu). 

The classic fits of Cohen and Kurath (1965) (CK) de- 
fined state of the art in the p shell for a long time. As they 
preceded the realistic G-matrices, they also contributed 
to hide the fact that the latter produce in ^°B a catastro- 
phe parallel to the one in ^■^Na. The work of Navratil and 
Ormand (2002), and Pieper, Varga, and Wiringa (2002) 
acted as a powerful reminder that brought to the fore the 
3b nature of the discrepancies. 

Once this is understood, the solution follows: Eq. (77) 
is assumed to be basically sound but the corrections are 
taken to be linear in the total number of valence particles 
m (the simplest form that a 3b term can take). 

Using / = (p3/2, 4/2, /7/2) generically in the 
{p, sd, pf) shells respectively, and r = pi/2 and r = 
c^3/2j5i/2 for the p and sd shells, Zuker (2003) proposes 
{k = Ko + {m — mo) Ki) 
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FIG. 20 The excitation spectrum of B for different interac- 
tions. 



VlfiR)- 1.5 kSto 



(78) 



where R stands for any realistic 2b potential. The results 
for ^"B are in Fig. 20. The black squares (NO) are from 
(Navratil and Ormand, 2002, Fig 4, 6hn). 
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FIG. 21 The excitation spectrum of Na for different inter- 
actions. 

The black circles correspond to the bare KLS G- matrix 
{hui =17 MeV), the white squares to the same with k = 
1.1, and the pentagons to the CK fit. NO and KLS give 
quite similar spectra, as expected from the discussion 
in Sections II. A and II. D. The k correction eliminates 
the severe discrepancies with experiment and give values 
close to CK. 

In ^^Na the story repeats itself: BC and KB are very 
close to one another, the ground state spin is again J = 1 
instead of J = 3 and the whole spectrum is awful. The k 
correction restores the levels to nearly correct positions, 
though USD still gives a better fit. This simple cure 
was not discovered earlier because k is not a constant. 
In ^"^Mg we could still do with the value for ^^Na but 



around ^^Si it must be substantially smaller. The spectra 
in Fig. 22 are obtained with K{m) — 0.9 — O.G5(m — 6), 
and now the spectra are as good as those given by USD. 
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FIG. 22 The excitation spectra of ' 
interactions. 



'Mg and ^^Si for different 



To determine a genuine 3b effect (i.e., the linearity 
of k) a sufficiently large span of A values is necessary. 
In the p shell, corrections to Vrr' become indispensable 
very soon and must be fitted simultaneously, as done 
successfully by Abzouzi et al. (1991), so we do not dwell 
on the subject. As we have seen, in the pf shell KB3, 
which is nearly perfect in ^4 = 47-50, must be turned into 
KB3G for A — 50-52, which also does well at the lower 
masses. To find a real problem with a 2b ii-interaction 
{i.e., compatible with NN data)^^ we have to move to 
56,58-^1 In particular the first B{E2){2 — > 0) transition 
in ^^Ni falls short of the observed value (140 e^fm^) by a 
factor w 0.4 with any monopole corrected KB interaction. 
The problem can be traced to weak quadrupole strength 
and it is not serious: as explained in Section II.D it is due 
to a normalization uncertainty, and Table I shows that 
with equal normalizations BC and KB in the pf shell are 



It may be possible to fit the data with a purely 2b set of matrix el- 
ements: USD is the prime example, but it is R-incompatible (Du- 
four and Zuker, 1996). 
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as close as in the sd shell. 

In (Zuker, 2003) BC was adopted. Small modifications 
of the Vfr matrix were made to improve the (already rea- 
sonable) r-spectrum in ^^Ni (see Fig. 19). The particular 
mixture in Eq. (78) was actually chosen to make possi- 
ble, in the simplest way, a good gap in ^^Ca and a good 
single-particle spectrum in *^Sc. It was found that for 
A = 48, K{m = 8) w 0.43. For A = 56, truncated calcu- 
lations yielded /t(16) ^ 0.28. The B(E2)(2 -> 0) in ^SNi 
indicated convergence to the right value. 



16 












14 












12 












10 
8 








If > 




6 












4 
2 

n 








BCk=0.43 - 
BCk=0.28 
KB3 
EXP 


▼ 

-■ - 

--• 



500 1000 1500 2000 2500 3000 3500 

A(Ej) (keV) 



FIG. 23 Backbending in See text. 

With K{m — 8), BC reproduces the yrast spectrum of 
^*Cr almost as well as KB3, but not with K{m = 16) 
(Fig. 23), indicating that the three-body drift is needed, 
though not as urgently as in the sd shell. Another indica- 
tion comes from the T=0 spectrum of "'^V, the counter- 
part of ^^Na and ^"B in the p/-shGll. The realistic inter- 
actions place again the 3=1 and J=3 states in the wrong 
order, the correct one is re-established by the three-body 
monopole correction. 



within 100'--'200 keV and to give a fairly good descrip- 
tion of the highly deformed excited band of this doubly 
magic nucleus (Mizusaki et at, 2002, 1999; Otsuka et al, 
1998). Other applications to the pf shell can be found 
in (Honma et ai, 1996). The existence of excited collec- 
tive bands in the N=28 isotones is studied in (Mizusaki 
et al, 2001). Very recently, a new interaction for the 
pf shell (GXPFl) has been produced by a Tokyo-MSU 
collaboration (Honma et ai, 2002), following the fitting 
procedures that lead to the USD interaction. The fit 
starts with the G-matrix obtained from the Bonn-C nu- 
cleon nucleon potential (Hjorth- Jensen et al., 1995). The 
fit privileges the upper part of the pf shell, as seen by 
the very large difference in single particle energies be- 
tween the lf7/2 and 12p3/2 orbits (3 MeV instead of the 
standard 2 MeV). In the calcium isotopes i.e., when only 
the T=l neutron neutron interaction is active) this new 
interaction retains the tendency of the bare G-matrices 
(KB, Bonn-C, KLS) to produce large gaps at N=32 and 
N=34, contrary to FPD6 or KB3G that only predict a 
large gap at N=32. Some early spectroscopic applications 
of GXPFl to the heavy isotopes of Titanium, Vanadium 
and Chromium can be found in (Janssens et ai, 2002) 
and (Mantica et ai, 2003), where the N=32-34 gaps are 
explored. 

Each nucleus has its interest, sometimes anecdotic, 
sometimes fundamental. The agreement with experiment 
for the energies, and for the quadrupole and magnetic 
moments and transitions is consistently good, often ex- 
cellent. These results have been conclusive in establish- 
ing the soundness of the minimally-monopole-modified 
realistic interaction(s). There is no point in reproducing 
them here, and we only present a few typical examples 
concerning spectroscopic factors, isospin non-conserving 
forces and "pure spectroscopy" . 



B. The pf shell 

Systematic calculations of the A=47-52 isobars can be 
found in (Caurier et ai, 1994, A = 48, KB3), (Martfnez- 
Pinedo et al, 1997, A = 47 and 49, KB3) and (Poves 
et al, 2001, A = 50-52, KB3G). Among the other 
full calculations let's highlight the following: The 
SMMC studies using either the FPD6 interaction (Al- 
hassid et al., 1994) or the KB3 interaction (Langanke 
et al, 1995). A comparison of the exact results with 
the SMMC can be found in (Caurier et al, 1999b). For 
a recollection of the SMMC results in the pf shell see 
also (Koonin et ai, 1997b). The recent applications of 
the exponential extrapolation method (Horoi et ai, 2002, 
2003). The calculations of Novoselsky et al. (1997) for 
^^Sc and ^^Ti using the DUPSM code (see also the erra- 
tum in Novoselsky et al, 1998b) and for ^^Sc and ^^Ti 
(Novoselsky et al, 1998a). The extensive QMCD calcu- 
lation of the spectrum of ^^Ni have been able to repro- 
duce the exact result for the ground state binding energy 



1. Spectroscopic factors 

The basic tenet of the independent particle model is 
that addition of a particle to a closed shell al\cs) pro- 
duces an eigenstate of the |cs -f 1) system. Nowadays we 
know better: it produces a doorway that will be frag- 
mented. If we choose \cs) = |^®Ca) and \cs + 1) = ["'^Sc), 
the four pf orbits provide the doorways. The lowest, 
fj/2, leads to an almost pure eigenstate. The middle 
ones, P3/2, i/2j are more fragmented, but the lowest level 
still has most of the strength and the fragments are scat- 
tered at higher energies. Fig. 24 shows what happens to 
the /5/2 strength: it remains concentrated on the door- 
way but splits locally. The same evolution with energy 
of the quasi-particles (Landau's term for single particle 
doorways) was later shown to occur generally in finite 
systems (Altshuler et al, 1997). 
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FIG. 24 Spectroscopic factors, {2j + l)S{j,tz), corresponding 
to stripping of a particle in the orbit I/5/2 (Marti'nez-Pinedo 
et al, 1997). 



2. Isospin non conserving forces 

Recent experiments have identified several yrast bands 
in mirror pf nuclei (Bentley et at, 1998, 1999, 2000; 
Brandolini et al., 2002b; Lenzi et al., 2001; O'Leary et al., 
1997, 2002) for A = 47, 49, 50 and 51. The naive view 
that the Coulomb energy should account for the MDE 
(mirror energy differences) turns out to be untenable. 
The four pairs were analyzed in (Zuker et al., 2002) where 
it was shown that three effects should be taken into ac- 
count. A typical result is proposed in Fig. 25, where Vcm 
and Vbm stand for Coulomb an nuclear isospin break- 
ing multipole contributions, while Vcm is a monopole 
Coulomb term generated by small differences of radii be- 
tween the members of the yrast band. The way these 
disparate contributions add to reproduce the observed 
pattern is striking. 
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FIG. 25 Experimental (O'Leary et al., 1997) and calcu- 
lated (Zuker et al, 2002) MED for the pair *^Cr-*^Mn. 

Another calculation of the isospin non conserving ef- 
fects, and their influence in the location of the proton drip 
line, is due to (Ormand, 1997). He has also analyzed the 



A=46 isospin triplet in (Garrett et al., 2001). 



3. Pure spectroscopy 

The first nucleus we have chosen is also one that has 
been measured to complete the mirror band (Bentley 
et al., 2000) in ^4 = 51, whose MED are as well described 
as those of A = 49 in Fig. 25. Such calculations need very 
good wavefunctions, and the standard test they have to 
pass is the "purely spectroscopic" one. 

Work on the subject usually starts with a litany. Such 
as: Quadrupole effective charges for neutrons q,y=0.5 and 
protons q^=1.5 and bare 5-factors g^=5.5857 /xjv, g^—- 
3.3826 /ijv, ffi=l-0 and ^[,=0.0 jiN in Ml transitions 
and moments are used. Except when the Ml transitions 
are fully dominated by the spin term, the use of effective 
g-factors does not modify the results very much due to 
the compensation between the spin and orbital modifica- 
tions. 

Then some spectrum: The yrast band of ^^Mn, calcu- 
lated in the full p/-shell space, is compared in Fig. 26 
with the experimental data. The first part of the test 
is passed. At most the examiner will complain about 
shghtly too high, high spin states. 
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FIG. 26 Yrast band of ^^Mn; experiment vs shell model cal- 
culation in the full p/-shell space. 

Finally the transitions: The transitions in Table IV 
are equally satisfactory. Note the abrupt change in both 
B[Ml) and B{E2) for J = (17/2)", beautifully repro- 
duced by the calculation. The origin of this isomerism is 
in the sudden alignment of two particles in the lf7/2 or- 
bit, which provides an intuitive physical explanation for 
the abrupt change in the MED (Bentley et al, 2000). 

At the end one can add some extras: The electromag- 
netic moments of the ground state are also known: Their 
values /Lte£cp=3.568(2)/.tAr and Qea;p=42(7) efm^ (Fire- 
stone, 1996) compare quite well with the calculated 
^th='^-S97^N and Qt;i=35 efm^. 
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TABLE IV Transitions 


in ^^Mn. 






Exp. 


Th. 


B{M1) 






7- 6- 
2 2 


U.zU ( (^04 ) 


n 1 '7'7 
U.i M 


9~ ^ I~ 

2 2 


0.16(5) 


0.116 


11- _^ 9- 

2 2 


0.662(215) 


0.421 


17 - ^ 15 - 
2 2 


0.00012(4) 


0.00003 


19- ^ il" 
2 2 


>0.572 


0.797 


B(£;2) 


(e' /m^) 


(e' /m*) 


7 — 5 — 
2 2 


528(146) 


305 


9- 6- 
2 ^2 


169(67) 


84 


9- 7- 
2 2 


303(112) 


204 


11 - 7- 

2 2 


236(67) 


154 


11- 9- 
2 2 


232(75) 


190 


17- 13- 
2 2 


1.236(337) 


2.215 



TABLE V Transitions in ^^Cr. (a) from (Brown et al., 1974). 



We complete this section with an even-even and an 
odd-odd nucleus. In figure 27 we show the yrast states of 
^^Cr, up to the band termination, calculated in the full 
p/-shell space. The agreement of the KB3G results with 
the experiment is again excellent. 
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FIG. 27 Yrast band of ^'^Cr; experiment vs shell model cal- 
culation in the full j3/-shell space. 

The electromagnetic moments of the first J = 2+ state 
are known: /iej.p=2.41(13) /ijy (Spcidcl et al, 2000) and 
Qexp = —8.2(16) e fm^ (Firestone, 1996) and they are 
in very good agreement with the theoretical predictions: 
/xt/i=2.496 yUjv and Qtft=-9.4 e fm^ In table V the experi- 
mental values for the electromagnetic transitions between 
the states of the yrast sequence arc given. The calculated 
values are in very good correspondence with the experi- 
ment. 

The yrast states of the odd-odd nucleus ^^Mn, calcu- 
lated in the full p/-shell, are displayed in fig. 28. No- 
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FIG. 28 Yrast band of ^'^Mn; experiment vs shell model cal- 
culation in the full p/-shell space. 



tice the perfect correspondence between theory and ex- 
periment for the states belonging to the multiplet below 
1 MeV. The experimental data for the spins beyond the 
band termination (11"*" to 16"*") come from a recent ex- 
periment (Axiotis, 2000). The agreement between theory 
and experiment is spectacular. There are only three ex- 
perimental states with spin assignment not drawn in the 
figure: a second J = (5+) at 1.42 MeV, a third J = (5+) 
at 1.68 MeV and a second J = (6+) at 1.96 MeV. The cal- 
culation places the second J = 5+ at 1.37 MeV, the third 
J = 5+ at 1.97 MeV and the second J = 6+ at 1.91 MeV. 
The experimental magnetic and quadrupole moments for 
the 6+ ground state are also known^ i-iexp — 

062(2) 

and Qexp=+50{7) c fm^ (Firestone, 1996). The calcu- 
lated values, /X(?(=2.952 /xjv and Qth=+50 e fm^ repro- 
duce them nicely. Some electromagnetic transitions along 
the yrast sequence have been measured. The calculations 
are in very good agreement with them (see table VI). 
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TABLE VI Transitions 


in ^^Mn. 






Exp. 


Th. 


B(M1) 
7+^6+ 
8+^7+ 
9+^8+ 


0.501(251) 

>0.015 
1.074lg;°« 


0.667 
0.405 
0.759 


B{E2) 
7+^6+ 
8+^6+ 
8+^7+ 
9+^7+ 
11+ ^ 9+ 


(eVm*) 
92l|?* 
>1.15 
>4.15 
104li°° 
54(6) 


(eVm^) 
126 
33 
126 
66 
53 



C. Gamow Teller and magnetic dipole strength. 

Out of the approximately 2500 known nuclei that are 
bound with respect to nucleon emission, only 253 are 
stable. The large majority of the rest decay by (3 emis- 
sion or electron capture, mediated by the weak interac- 
tion. When protons and neutrons occupy the same or- 
bits, as in our case, the dominant processes arc allowed 
Fermi and Gamow- Teller transitions. The information 
obtained from the weak decays has been complemented 
by the {p,n) and {n,p) reactions in forward kinematics, 
that make it possible to obtain total GT strengths and 
strength functions that cannot be accessed by the decay 
data because of the limitations due to the Qp windows. 
From a theoretical point of view, the comparison of calcu- 
lated and observed strength functions provide invaluable 
insight into the meaning of the valence space and the na- 
ture of the deep correlations detected by the "quenching" 
effect. 

The half- life for a transition between two nuclear states 
is given by (Behrens and Biihring, 1982; Schopper, 1966): 



ifA+nt 



6144.4 ± 1.6 



{fv/fA)B{F)+B{GTy 



(79) 



The value 6144.4 ± 1.6 is obtained from the nine best- 
known supcrallowcd beta decays (Towner and Hardy, 
2002) (see Wilkinson, 2002a, b, for an alternative study). 
fv and /a are the Fermi and Gamow- Teller phase-space 
factors, respectively (Chou et ai, 1993; Wilkinson and 
Macefield, 1974). is the phase space for electron cap- 
ture (Bambynek et ai, 1977), that is only present in /?+ 
decays. If ti/2 is the total lifetime, the partial lifetime of 
a level with branching ratio br is t = ti/2/br- 



B{F) and B{GT) are defined as 

"(/ II Efc4 II iy' 



B{F) 



B{GT) = 



V2Jr + 1 

9a\ if II Ek'^'ti II ^) 
9v 



V2jrn" 



(80a) 
(80b) 



Matrix elements are reduced (A3) with respect to spin 
only, ± refers to decay, a = 2S and {gA/gv) = 
-1.2720(18) (Hagiwara et ai, 2002) is the ratio of the 
weak interaction axial-vector and vector coupling con- 
stants. 

For states of good isospin the value of B{F) is fixed. It 
can be altered only by a small isospin-symmctry-breaking 
Sc correction (Towner and Hardy, 2002). Shell model 
estimates of this quantity can be also found in (Ormand 
and Brown, 1995). 



B{F) = [T(r+1)-T,.T,,] %(l-5c), 



(81) 



where Si / allows only transitions between isobaric analog 
states. Superallowed decays may shed light on the depar- 
tures from unitarity of the Cabibbo Kobayashi Maskawa 
matrix. 

The total strengths S± are related by the sum rules 

S- (F) -S+{F) = N- Z, (82a) 
S-{GT) - S+{GT) = 3{N - Z), (82b) 

where N and Z refer to the initial state and S±{GT) 
does not contain the Qa/qv factor. The comparison of 
the (p, n) and (n, p) data with the Gamow Teller sum rule 
(Ikeda et al., 1963) led to the long standing "quenching" 
problem; only approximately one half of the sum rule 
value was found in the experiments. 

The GT strength is not protected by a conservation 
principle and depends critically on the wavefunctions 
used. Full Ohu) calculations already show a large quench- 
ing with respect to the independent particle limit as seen 
in Table VII where the result of a full pf calculation is 
compared with that obtained with an uncorrelated Slater 
determinant having the same occupancies: 



{2ji + l){2jk + 1) 



(83) 



the sum on i runs over the proton (neutron) orbits in the 
valence space and k over the proton (neutron) orbits for 
5+ (-?-). and n'' denote the number of particles and 
holes, respectively. The determinantal state demands a 
quenching factor that is almost twice as large as the stan- 
dard quenching factor = (0.74)^) that brings the full 
calculation in line with experiment. 



1. The meaning of the valence space 

Before moving to the explanation of the quenching of 
the GT strength, it is convenient to recall the meaning 
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TABLE VII Comparison of GT+ strengths. For ^''Fe, ^^Mn, 
^*Ni and ^^Co the calculations are truncated to t = 8, t = 4, 
t = 6 and t = 4 respectively. The data are from (Alford et al, 
1993; El-Kateb et al., 1994; Vetterh et al, 1990; Williams 
et al., 1995). 



Nucleus Uncorrelated 



Correlated 



Expt. 



Unquenched Q = 0.74 



51y 


5.15 


2.42 


1.33 


1.2 ±0.1 


5%e 


10.19 


5.98 


3.27 


3.3 ±0.5 




7.96 


3.64 


1.99 


1.7 ±0.2 


56pe 


9.44 


4.38 


2.40 


2.8 ±0.3 


5«Ni 


11.9 


7.24 


3.97 


3.8 ±0.4 


5«Co 


8.52 


3.98 


2.18 


1.9 ±0.1 


62Ni 


7.83 


3.65 


2.00 


2.5 ±0.1 



of the valence space as discussed in Section II.A.l. The 
^*Ca(p, n)*^Sc reaction (Anderson et al., 1990) provides 
an excellent example. In Fig. 29, adapted from Caurier 
et al. (1995b), the experimental data are compared with 
the strength function produced by a calculation in the 
full p/-shell using the interaction KB3. The peaks have 
all J = 1, T = 3. In the pf shell there are 8590 of them 
and the calculation has been pushed to 700 iterations in 
the Lanczos strength function to ensure fully converged 
eigenstates below 11 MeV. Of these eigenstates 30 are 
below 8 MeV: They are at the right energy and have the 
right strength profile. At higher energies the peaks are 
much too narrow compared with experiment. It means 
that they may well be eigenstates of the effective Hamil- 
tonian in the pf shell, but not eigenstates of the full 
system. Therefore, they should be viewed as doorways, 
subject to further mixing with the background of intrud- 
ers which dominates the level density after 8 MeV, as 
corroborated by the experimental tail that contains only 
intruders and can be made to start naturally at that en- 
ergy. 

The KB3 effective interaction is doing a very good job, 
but it is certainly not decoupling 8590 pf states from 
the rest of the space. If the fact is not explicitly rec- 
ognized we end up with the often raised (Hjorth- Jensen 
et al., 1995) "intruder problem": Decoupling cannot be 
enforced perturbatively when intruders are energetically 
close to model states. Figure 29 indicates that although 
few eigenstates are well decoupled, it is possible to make 
sense of many others if one interprets them as doorways. 
The satisfactory description of the lowest states indicates 
that the model space makes sense. It means that it en- 
sures good decoupling at the 5*2 level, or simply in sec- 
ond order perturbation theory, which guarantee a state- 
independent interaction. 

Energetically we are in good shape, and we concentrate 
on the renormalization of the GT operator: the calcu- 
lated strength has been quenched by a factor « (0.74)^. 
Why? Why this factor ensures the right detailed strength 
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FIG. 29 Gamow Teller strength in ***Ca(p, n)^**Sc from (An- 
derson et al., 1990) — after elimination of the Fermi peak at 
around 6 MeV — compared with the calculated peaks (KB3 in- 
teraction) after 700 Lanczos iterations (Caurier et al., 1995b). 
The peaks have been smoothed by gaussians having the in- 
strumental width of the first measured level. 



for about 30 states below 8 MeV? 



2. Quenching 

To understand the quenching problem it is best to start 
by the (tentative) solution. The dressed states in Eq. (4) 
are normalized to unity in the model space. This trick is 
essential in the formulation of linked cluster or exp S the- 
ories (hence Eq. (63)). It makes possible the calculation 
of the energy — and some transitions, such as the E2 — 
without knowledge of the norm of the exact wavefunc- 
tion. In general though, we need an expectation value 
between exact, normalized states: (/ || T \\ i)'^ . If we 
write 



— a\Qfiw) 



(84) 



and a similar expression in a' , (3' for (/|, we find 

2 

(/||T|| IY= |a«'ro±^/3„/3;T„ 



(85) 



since the GT operator does not couple states with differ- 
ent number of hw excitations. If we make two assump- 
tions: (a) neglect the n ^ Q contributions; (b) a w a'; 
it follows that if the projection of the physical wavefunc- 
tion in the Ohuj space is Q ^ a^, its contribution to the 
transition will be quenched by . 

Exactly the same arguments apply to transfer 
reactions — for which T = as (or a\) — but are simpler 
because T„^o = 0. The transition strength is given by 
the spectroscopic factor, which can be identified with Q 
when one particle is removed and 1 — Q when it is added. 
The assumption that the model amplitudes in the ex- 
act wavefunctions are approximately constant is borne 
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out by systematic calculations of Q in the p shell (Chou 
et at, 1993, Q = 0.820(15)), the sd shell (Wildcnthal 
et al, 1983, Q = 0.77(2)) and the pf shell (illustrated in 
Fig. 30, Martmez-Pinedo et al., 1996b, Q = 0.744(15)). 




0.4 0.6 

T(GT) Th. 



1.0 



FIG. 30 Comparison of the experimental and theoretical val- 
ues of the quantity T{GT) in the pf shell (Martmez-Pinedo 
et al., 1996b). The x and y coordinates correspond to theo- 
retical and experimental values respectively. The dashed line 
shows the "best-fit" for Q — 0.744. The solid line shows the 
result obtained in the sd-shell nuclei (Wildenthal et al, 1983) 



These numbers square well with the existing informa- 
tion on spectroscopic factors from {d,p) (Void et at, 1978, 
Q w 0.7) and (e, e'p) (Cavedon et al, 1982, Q = 0.7) data 
(see also Pandharipande et al, 1997). This consistency is 
significant in that it backs assumption (a) above, which 
is trivially satisfied for spectroscopic factors. It opens 
the perspective of accepting the GT data as a measure 
of a very fundamental quantity that does not depend 
on particular processes. The proposed "solution" to the 
quenching problem amounts to reading data. 
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FIG. 31 Effective spin g-factor of the Ml operator deduced 
from the comparison of shell model calculations and data for 
the total B(M1) strengths in the stable even mass N=28 iso- 
tones (from von Neumann-Cosel et al. (1998)) 

Experimentally, the challenge is to locate all the 



strength, constrained by the Ikeda sum rule that relates 
the direct and inverse processes. The careful analysis 
of Anderson et al. (1985) suggests, but does not prove, 
that the experimental tail in Fig. 29 contains enough 
strength to satisfy approximately the sum rule. A similar 
result is obtained for ^^Fc{p,n) (Anderson et al, 1990). 
More recent experiments by the Tokyo group establish 
that the strength located at accessible energies exhausts 
90(5)% of the sum rule in the ^°Zr(n,p) and 84(5)% in 
27Al(p,n) (Wakasa et al, 1998, 1997). 

The theoretical problem is to calculate Q. It has been 
compounded by a sociological one: the full GT operator 
is {gA/ gv)o'T, where gA/gv ~ —1.27 is the ratio of weak 
axial and vector coupling constants. The hotly debated 
question was whether Q was due to non nucleonic renor- 
malization of gA or nuclear renormalization of ar (Arima, 
2001; Osterfeld, 1992). We have sketched above the nu- 
clear case, along the lines proposed by Caurier, Poves, 
and Zuker (1995b), but under a new guise that makes 
it easier to understand. The calculations of Bertsch and 
Hamamoto (1982), Drozdz et al. (1986) and Dang et al 
(1997) manage to place significant amounts of strength 
beyond the resonance region, but they arc based on 2p-2h 
doorways which fall somewhat short of giving a satisfac- 
tory view of the strength functions. No-core calculations 
are under way, that should be capable of clarifying the 
issue. 

The purely nuclear origin of quenching is borne out by 
(p.p'), (7,7') and (e,e') experiments that determine the 
spin and convection currents in Ml transitions, in which 
9A/gv play no role (see Richter, 1995, for a complete re- 
view). An analysis of the data available for the N = 28 
isotones in terms of full p/-shell calculations concluded 
that agreement with experiment was achieved by quench- 
ing the ar operator by a factor 0.75(2), fully consistent 
with the value that explains the Gamow- Teller data (von 
Neumann-Coscl et al, 1998) (see Fig. 31). These re- 
sults rule out the hypothesis of a renormalization of the 
axial- vector constant qa- it is the err operator that is 
quenched. 



VI. SPHERICAL SHELL MODEL DESCRIPTION OF 
NUCLEAR ROTATIONS 

Theoretical studies in the pf shell came in layers. The 
first, by McCuUen, Bayman, and Zamick (1964), re- 
stricted to the /7/2 space, was a success, but had some 
drawbacks: the spectra were not always symmetric by in- 
terchange of particles and holes, and the quadrupole mo- 
ments had systematically the wrong sign. The first diag- 
onalizations in the full shell [(Pasquini, 1976), (Pasquini 
and Zuker, 1978)] solved these problems to a large extent, 
but the very severe truncations necessary at the time 
made it impossible to treat on the same footing the pair- 
ing and quadrupole forces. The situation improved very 
much by dressing perturbatively the pure two-body f'j/2 
part of the Hamiltonian H2, with a three body term Hm, 
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mostly due to the quadrupolc force (Povcs and Zuker, 
1981b). The paper ended with the phrase: 

"It may well happen, that in some cases, not in the pf 
shell but elsewhere, Hm will overwhelm H2. Then, and 
we are only speculating, we shall speak, perhaps, of the 
rotational coupling scheme." 

Indeed, some nuclei were indicating a willingness to 
become rotational but could not quite make it, simply 
because the perturbative treatment was not enough for 
them. The authors missed prophecy by one extra condi- 
tional: "not in the pf shell should have been "not neces- 
sarily..." 

At the time it was thought impossible to describe ro- 
tational motion in a spherical shell-model context. The 
glorious exception, discovered by Elliott (1958a, b), was 
(apparently) associated to strict SU(3) symmetry, ap- 
proximately realized only near ^"^Ne and ^^Mg. 



A. Rotors in the pf shell 



TABLE VIII ''^Cr; quadrupole properties of the yrast band 
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The fourth layer was started when the ANTOINE code 
(Section III.D) came into operation: ^^Cr was definitely 
a well deformed rotor (Caurier et aZ., 1994), as can be 
surmised by comparing with the experimental spectrum 
from Lenzi et al. (1996) in Fig. 32 (see also Cameron 
et al., 1993) and the transition properties in Table VIII 
from Brandolini et al. (1998) where we have used 



(J + l)(2J + 3) 

B{E2, J ^ J - 2) = 



K ^ 1 (86) 



— e^\{JK20\J-2,K)\^Qo{tf K^l/2,1; (87) 

lOTT 

to establish the connection with the intrinsic frame de- 
scriptions: A good rotor must have a nearly constant Qo, 
which is the case up to J = 10, then ''^Cr backbends. 
This is perhaps the most striking result to come out of 
our pf calculations, because such a behaviour had been 
thought to occur only in much heavier nuclei. 

Fig. 33 compares the experimental patterns with those 
obtained with KB3, and with the Gogny force. The lat- 
ter, when diagonalized gives surprisingly good results. 



FIG. 32 ''^Cr level scheme; experiment vs. theory 



When treated in the cranked Hartree Fock Bogoliubov 
(CHFB) approximation, the results are not so good. The 
discrepancy is more apparent than real: The predictions 
for the observables are very much the same in both cases 
(see Caurier et al., 1995a, for the details). The reason 
is given in Fig 34, where exact KB3 diagonalizations are 
done, subtracting either of the two pairing contributions, 
JT = 01 and 10 (Poves and Martmez-Pinedo, 1998). It 
is apparent that the JT = 01-subtracted pattern is quite 
close to the CHFB one in Fig. 33, especially in the ro- 
tational regime before the backbend. The JT = 10 sub- 
traction goes in the same direction. The interpretation is 
transparent: CHFB does not "see" proton- neutron pair- 
ing at all, and it is not very efficient in the low pair- 
ing regime. As it does everything else very well, the in- 
evitable conclusion is that pairing can be treated in first 
order perturbation theory, i.e., the energies are very sen- 
sitive to it, but not the wavefunctions. Floods of ink have 
gone into "neutron-proton" pairing, which is a problem 
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FIG. 33 The yrast band of "^Cr; experiment vs. the sheU FIG. 35 The yrast band of ^°Cr; experiment vs. the shell 
model calculations with KB3 and the Gogny force and the model calculations with the KB3 interaction 
CHFB results also with the Gogny force. 
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FIG. 34 Gamma-ray energies along the yrast band of *®Cr (in 
MeV), full interaction (KB3), isoscalar pairing retired (KB3- 
POl) and isovector pairing retired (KB3-P01). 

for mean field theories but neither for the Gogny force 
nor for the shell model. Furthermore, the results show 
that ordinary pairing is also a mean field problem when 
nuclei are not superfluids. 

**Cr has become a standard benchmark for mod- 
els (Cranked Hartree Fock Bogoliuvov (Caurier et ai, 
1995a), Cranked Nilsson Strutinsky (Juodagalvis and 
Aberg, 1998), Projected Shell Model (Kara et at, 1999), 
Cluster Model (Descouvemont, 2002), etc.) Nuclei in the 
vicinity also have strong rotational features. The mirror 
pairs "'^V-^'^^Cr and ''^Cr-^^Mn closely follow the semiclas- 
sical picture of a particle or hole hole strongly coupled to 
a rotor (Martmez-Pinedo et ai, 1997) in full agreement 
with the experiments (Bentley et at, 1998; Cameron 
et ai, 1991, 1994; O'Leary et at, 1997; Tonev et ai, 
2002). '^"Cr was predicted to have a second backbending 
by Martmez-Pinedo et al. (1996a), confirmed experimen- 
tally (Lenzi et ai, 1997) (see figure 35). When more par- 
ticles or holes are added, the collective behaviour fades, 
though even for ^^Fe, a rotor like band appears at low 
spin with an yrast trap at J=12"'", both are accounted 



for by the shell model calculations (Poves and Zuker, 
1981b)(sic), (Ur et ai, 1998). For spectroscopic com- 
parisons with the odd-odd nuclei, see Brandolini et al. 
(2001); Lenzi et al. (1999) for ■**'V and Svensson et al. 
(1998) for ^"Mn. Recently, a highly deformed excited 
band has been discovered in ^^Ni (Rudolph et ai, 1999). 
It is dominated by the configuration (lf7/2)^^ (2p3/2, 
lf5/2, 2pi/2)''- The calculations reproduce the band, that 
starts at about 5 MeV excitation energy and has a defor- 
mation close to (3=0.A. 



B. Quasi-SU(3). 

To account for the appearance of backbending rotors, 
a theoretical framework was developed by Zuker et al. 
(1995), and made more precise by Martmez-Pinedo et al. 
(1997). Here we give a compact overview of the scheme 
that will be shown to apply even to the classic examples 
of rotors in the rare-earth region. 

Let us start by considering the quadrupole force alone, 
taken to act in the p-th oscillator shell. It will tend to 
maximize the quadrupole moment, which means filling 
the lowest orbits obtained by diagonalizing the operator 
Qo = 2(720 = — — y"^. Using the cartesian repre- 
sentation, 2^20 ~ 'i.riz — Tlx ~ Uy^ we find eigenvalues 2p, 
2p— 3,. . . , etc., as shown in the left panel of Fig. 36, where 
spin has been included. By filling the orbits orderly we 
obtain the intrinsic states for the lowest SU(3) represen- 
tations (Elliott, 1958a, b): (A, 0) if all states are occupied 
up to a given level and (A, /i) otherwise. For instance: 
putting two neutrons and two protons in the K = \/2 
level leads to the (4p,0) representation. For four neu- 
trons and four protons, the filling is not complete and we 
have the (triaxial) (8(p — 1),4) representation for which 
we expect a low lying 7 band. 

In jj coupling the angular part of the quadrupole op- 
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where 



FIG. 36 Nilsson orbits for SU(3) (fc = -2p) and quasi-SU(3) 
(fe = -2p + l/2)). 



erator q^^ = r^C"^^ has matrix elements 

3[(j + 3/2)2 -m^] 



{jm\C^\j + 2m)^ 
{jm\C^\j + lm) = - 



2{2j + 3)2 ' 
3m[(j + l)2-m2]V2 
2.7(2.7 + 2)(2j+ 4) 



(88) 
(89) 



The Aj = 2 numbers in Eq. (88) are — within the ap- 
proximation made — identical to those in LS scheme, ob- 
tained by replacing ) by I. The A) = 1 matrix elements 
in Eq. (88) are small, both for large and small m, corre- 
sponding to the lowest oblate and prolate deformed orbits 
respectively. If the spherical )-orbits arc degenerate, the 
Aj — 1 couplings, though small, will mix strongly the 
two Aj = 2 sequences (e.g., (/7/2P3/2) and (/5/2Pi/2))- 
The spin-orbit splittings will break the degeneracies and 
favour the decoupling of the two sequences. Hence the 
idea (Zuker et ai, 1995) of neglecting the Aj = 1 matrix 
elements and exploit the correspondence 



I 



j = 1 + 1/2 m — > m + 1/2 x sign(m). 



which is one-to-one except for m = 0. The resulting 
"quasi SU(3)" quadrupole operator respects SU(3) rela- 
tionships, except for m = 0, where the correspondence 
breaks down. The resulting spectrum for quasi-2(720 is 
shown in the right panel of Fig 36. The result is not 
exact for the K = 1/2 orbits but a very good approxi- 
mation. 

To check the validity of the decoupling, a Hartree cal- 
culation was done for H = eHgp + Ha, where Hgp is 
the observed single particle spectrum in ^^Ca (essentially 
equidistant orbits with 2MeV spacings) and Hq is the 
quadrupole force in Eq. (33) with a properly renormal- 
ized coupling. The result is exactly a Nilsson (1955) cal- 
culation (Martmez-Pinedo et ai, 1997), 



H 



Nilsson 



^2920 



(90) 



1 (2g20) 

4 (r2) 



(2(72 



(p + 3/2) 



4 ■ 



(91) 



In the right panel of Fig. 37 the results are given in the 
usual form. 
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FIG. 37 Nilsson diagrams in the pf shell. Energy vs. sin- 
gle particle splitting e (left panel), energy vs. deformation 5 
(right panel). 

In the left panel we have turned the representation 
around: since we are interested in rotors, we start from 
perfect ones (SU(3)) and let e increase. At a value of 
~ 0.8 the four lowest orbits are in the same sequence as 
the right side of fig. 36 (Remember here that the real 
situation corresponds to e « 1.0). The agreement even 
extends to the next group, although now there is an in- 
truder (1/2 [3 10] orbit). The suggestion is confirmed by 
an analysis of the wavefunctions: For the lowest two or- 
bits, the overlaps between the pure quasi-SU(3) wave- 
functions calculated in the restricted A) = 2 space (fp 
from now on) and the ones in the full pf shell exceeds 
0.95 throughout the interval 0.5 < e < 1. More interest- 
ing still: the contributions to the quadrupole moments 
from these two orbits vary very little, and remain close 
to the values obtained at e = i.e., from fig. 36). 

We have learned that -for the rotational features- cal- 
culations in the restricted {fp)" spaces account remark- 
ably well for the results in the full major shell (pf)" (see 
last column of Table VIII). Let us move now to larger 
spaces. In Fig. 38 we have yrast transition energies for 
different configurations of 8 particles in A) = 2 spaces. 
The force is KLS, hui = 9 MeV, the single particle split- 
tings uniform at e = 1 MeV, and gds, say, is the lower 
sequence in the sdg, p = 4, shell. Rotational behav- 
ior is fair to excellent at low J. As expected from the 
normalization property of the realistic quadrupole force 
[Eq. (35)] the moments of inertia in the rotational region 
go as {p + 3/2)2 _^ 3/2)2, ^_g_ jf multiply all the 
Ej values by this factor the lines become parallel. The 
intrinsic quadrupole moment Qo [Eq. (86)] remains con- 
stant to within 5% up to a critical J value at which the 
bands backbend. 

Why and how do the bands backbend? We have no 
simple answer, but Fig. 39 from (Velazquez and Zuker, 
2002) shows the behavior of the (gds)^ space under the 
influence of a symmetric random interaction, gradually 
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FIG. 38 Yrast transition energies = E{J + 2) - E{J) for 
different configurations, KLS interaction. 
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FIG. 39 Backbending patterns in {gdsf T = 0, with the KLS 
interaction and with a random interaction plus a constant 

Wri^tu = -a. 



made more attractive by an amount a. There is no need 

to stress the similarity with Fig. 38. The appearance of 
backbending rotors seems to be a general result of the 
competition between deformation and alignment, char- 
acteristic of nuclear processes. 

The group theoretical aspects of quasi-SU3 have been 
recently discussed and applied to the description of the 
sd-shell nuclei by Vargas, Hirsch, and Draayer (2001, 
2002a). 



C. Heavier nuclei: quasi+pseudo SU(3) 

We have seen that quasi-SU(3) is a variant of SU(3) 
that obtains for moderate spin-orbit splittings. For 

other forms of single particle spacings, the pseudo-SU(3) 
scheme (Arima et al. (1969); Hecht and Adler (1969); 




FIG. 40 Schematic single particle spectrum above ^^^Sn. rp 
is the set of orbits in shell p excluding the largest. For the 
upper shells the label / is used for j = / + 1/2 



Raju et al. (1973), see also Vargas et al. (2002b) for 
more recent applications) will be favored (in which case 
we have to use the left panel of Fig. 36, with pseudo- 
p=p—l). Other variants of SU(3) may be possible and 
are well worth exploring. In cases of truly large deforma- 
tion SU(3) itself may be valid in some blocks. 

To see how this works, consider Fig. 40 giving a 
schematic view of the single particle energies in the space 
of two contiguous major shells — in protons (tt) and neu- 
trons {v) — adequate for a SM description of the rare 
earth region. 

We want to estimate the quadrupolc moments for nu- 
clei at the onset of deformation. We shall assume quasi- 
SU(3) operates in the upper shells, and pseudo-SU(3) in 
the lower ones. The number of particles in each shell 
for which the energy will be lowest will depend on a bal- 
ance of monopole and quadrupole effects, but Nilsson 
diagrams suggest that when nuclei acquire stable defor- 
mation, two orbits K=l/2 and 3/2 originating in the 
upper shells of Fig. 36 become occupied, i.e., the up- 
per blocks are precisely the 8-particle configurations we 
have studied at length. Their contribution to the electric 
quadrupole moment is then 



Qo = 8[e„(p7r - 1) + e^iPu - 1)], 



(92) 



with = 5, p„ = 6; e^r and e„ are the effective charges. 

Consider even-even nuclei with Z=60-66 and N=92- 
98, corresponding to 6 to 10 protons with pseudo-p = 3, 
and 6 to 10 neutrons with pseudo-p = 4 in the lower 
shells. From the left part of Fig. 36 we obtain easily 
their contribution to Qo, which added to that of Eq. (92) 
yields a total 



Qo = 56e^ -t- (76 -|- 4n)ei,, 



(93) 
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TABLE IX B{E2) T in e^6^ compared with experiment (Ra- 
man et al., 1989) 



IVT 

IN 


iNcl 


om 




Dy 


92 


4.47 


4.51 


4.55 


4.58 




2.6(7) 


4.36(5) 


4.64(5) 


4.66(5) 


94 


4.68 


4.72 


4.76 
5.02(5) 


4.80 
5.06(4) 


96 


4.90 


4.95 


4.99 
5.25(6) 


5.03 

5.28(15) 


98 


5.13 


5.18 


5.22 


5.26 
5.60(5) 



for 152+2nNj^ 154+2ngj^^^ 156+2nQ^ 158+2nj3y 

spectively. At fixed n, the value is constant in the 
four cases because the orbits of the triplet K=l/2, 
3/2, 5/2 in Fig. 36 have zero contribution for p=3. 
Qo (given in dimensionless oscillator coordinates, i.e., 
r — > r/b with 6^ w fw?), is related to the 

E2 transition probability from the ground state by 
B{E2) t= 10-^A2/3q2^ The results, using effective 
charges of e^r = 1.4, e„ = 0.6 calculated in (Dufour and 
Zuker, 1996) are compared in table IX with the available 
experimental values . The agreement is quite remarkable 
and no free parameters are involved. Note in particu- 
lar the quality of the prediction of constancy (or rather 
A^/^ dependence) at fixed n, which does not depend on 
the choice of effective charges. The discrepancy in ^^^Nd 
is likely to be of experimental origin, since systematics 
indicate, with no exception, much larger rates for a 2+ 
state at such low energy (72.6 keV). It is seen that by 
careful analysis of exact results one may come to very 
simple computational strategies. In the last example on 
B(E2) rates, the simplicity is such that the computation 
reduces to a couple of sums. 



D. The Ar and Ca super-deformed bands 

The arguments sketched above apply to the region 
around ^^O, where a famous 4-particles 4-holes (4p-4h) 
band starting at 6.05 MeV was identified by Carter 
et al. (1964), followed by the 8p-8h band starting at 
16.75 MeV (ChevaUier et al, 1967). SheU model cal- 
culations in a very small space, ^1/2^5/2^1/21 could ac- 
count for the spectroscopy in ^®0, including the 4p-4h 
band (Zuker et at, 1968, ZBM), but the 8p-8h one needs 
at least three major shells and was tackled by an a-cluster 
model (Abgrall et ai, 1967a, b). It is probably the first 
superdeformed band detected and explained. 

In '*°Ca, the first excited 0+ state is the 4p-4h band- 
head. It is only recently that another low-lying highly 
deformed band has been found (Ideguchi et ai, 2001), fol- 
lowing the discovery of a similar structure in '^^Ar (Svens- 
son et al, 2000). 




Ql , \ , \ , 1 , \ , 1 

1000 2000 3000 4000 5000 

Ey (keV) 

FIG. 41 The superdeformed band of '^^Ar, experiment vs. SM 



By applying the quasi-|-pseudo-SU(3) recipes of Sec- 
tion VI. C we find that the maximum deformations at- 
tainable are of 8p-8h character in '^''Ca and 4p-4h in '^^Ar 
with Qo ~ 180 efm^ and Qq = 136 efm^ respectively. 

The natural generalization of the ZBM space consists 
in the d3/2Si/2/p orbits (remember: fp = /7/2P3/2), 
which is relatively free of center of mass spuriousncss. 
And indeed, it describes well the rotational regime of the 
observed bands. However, to track them beyond back- 
bend, it is convenient to increase the space to ^3/251/2^/- 
The adopted interaction, sdfp.sm is the one originally 
constructed by Retamosa et al. (1997), and used in (Cau- 
rier et al, 1998) (USD, KB3 for intra shell and KLS for 
cross shell matrix elements), with minor monopole ad- 
justments dictated by new data on the single particle 
structure of ^'^Si (Nummela et ai, 2001a). 

The calculations are conducted in spaces of fixed num- 
ber of particles and holes. In figure 41 the calculated 
energy levels in '^^Ar are compared to the data. The 
agreement is excellent, except at J = 12 where the data 
show a clear backbending while the calculation produces 
a much smoother upbending pattern. 

In table X the calculated spectroscopic quadrupole 
moments (Qs) and the B(E2)'s are used to compute 
the intrinsic Q as in Eqs. (86,87), using standard ef- 
fective charges SqTr=Sq^=0.5. As expected both Qo(s) 
and Qo(t) are nearly equal and constant — and close to 
the quasi-l-pseudo SU(3) estimate — up to the backbend. 
The calculated B(E2)'s agree well with the experimental 
ones (Svensson et ai, 2001). The value of Qq corresponds 
to a deformation /3 « 0.5. 

Now we examine the 8p-8h band in ^"Ca (Ideguchi 
et ai, 2001). The valence space adopted for '^^Ar is trun- 
cated by limiting the maximum number of particles in 
the lf5/2 and 2pi/2 orbits to two. 

The experimental (Ideguchi et al, 2001) and calcu- 
lated yrast gamma-ray energies are compared in fig. 42. 
The patterns agree reasonably well but the change of 
slope at J = 10 — where the backbend in ''^Cr starts 
(Fig. 33) — is missed by the calculation, which only back- 
bends at J = 20, the band termination for the configura- 
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TABLE X Quadrupole properties of the 4p-4h configuration's 
yrast-band in ^®Ar (in e^fm* and efm'^) 



B(E2)(J->J-2) 
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Qo(s) 


Qo(t) 


2 




315 


-36.0 


126 


126 
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372(59) 
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-45.9 


126 
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-50.7 
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120 
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-52.8 
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10 
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-52.7 


121 


104 


12 
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315 
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96 
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-54.3 


120 


82 
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FIG. 42 The superdeformed band in ''"Ca; exp. vs. 8p-8h 
calculation 

tion /|/2(^3/2'Si/2)~*- The extra collectivity induced by 
the presence of sd particles in pseudo-SU(3) orbitals is 
responsible or the delay in the alignment, but apparently 
too strong to allow for the change in slope at J = 10. 

The experimental Qo(t)=180l29 obtained from the 
fractional Doppler shifts corresponds to a deformation 
P ~ 0.6 (Ideguchi et ai, 2001). It is extracted from an 
overall fit that assumes constancy for all measured values, 
and corresponds exactly to the quasi-|-pseudo SU(3) es- 
timate. It also squares well with the calculated 172 efm^ 
in Table XI where the steady decrease in collectivity re- 
mains consistent with experiment within the quoted un- 
certainties. A reanalysis of the experimental lifetimes in 
the superdeformed band of '^"Ca (Chiara et ai, 2003) sug- 
gest that for low spins the deformation would be smaller 
due to the mixing with less deformed states of lower np- 
nh rank. 

This calculation demonstrates that a detailed descrip- 
tion of very deformed bands is within reach of the shell 
model. The next step consists in remembering that states 
at fixed number of particles are doorways that will frag- 
ment in an exact calculation that remains to be done. 

The band in '^^Ar has been also described by the pro- 
jected shell model in (Long and Sun, 2001). Beyond mean 



TABLE XI Quadrupole properties of the 8p-8h configura- 
tion's yrast-band in ^''Ca (in e^fm* and efm'^), calculated in 
the sdpf valence space 
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B(E2)(J J-2) 
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111 


162 


20 
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-85.1 






22 
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-79.1 






24 
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-81.5 







field methods using the Skyrme interaction have been re- 
cently applied to both the '^^Ar and ^°Ca superdeformed 
bands by Bender, Flocard, and Heenen (2003a). 



E. Rotational bands of unnatural parity 

The occurrence of low-lying bands of opposite parity to 
the ground state band is very frequent in pf shell nuclei. 
Their simplest characterization is as particle-hole bands 
with the hole in the ld3/2 orbit. In a nucleus whose 
ground state is described by the configurations (p/)" the 
opposite parity intruders will be (ld3/2)~^ (p/)"^^- The 
promotion of a particle from the sd to the pf shell costs 
an energy equivalent to the local value of the gap, that in 
this region is about 7 MeV. On the other side, the pres- 
ence of one extra particle in the p/-shell may produce 
an important increase of the correlation energy, that can 
compensate the energy lost by the particle hole jump. 
For instance, the very low-lying positive parity band of 
''''V, can be interpreted as (ld3/2)~^ (a proton hole) cou- 
pled to (pf^) T=0. Indeed, the correlation energy of this 
pseudo-'^^Cr is larger than the correlation energy of the 
ground state of ^^V and even larger than the correlation 
energy of the real ^*Cr, explaining why the band starts 
at only 260 keV of excitation energy (Poves and Sanchez 
Solano, 1998). The most extreme case is ''^Sc, where the 
intruder band based in the configuration (Id3y2)^^ cou- 
pled to ''^Ti, barely miss (just by 12 keV) to become the 
ground state. Many bands of this type have been ex- 
perimentally studied in recent years, mainly at the Gasp 
and Euroball detectors (Brandolini et ai, 1999), and ex- 
plained by shell model calculations. 

We present now some recent results in (Brandolini 
et ai, 2002a). The positive parity levels of this odd- 
odd nucleus are very well reproduced by the shell model 
calculation using the interactions KB3 (or KB3G) as can 
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FIG. 43 The ground state (K=4+) and the side band (K=l+) 
of **V, experiment vs. theory. 



be seen in Fig. 43, another example of good quality "pure 
spectroscopy". The levels are grouped in two bands, one 
connected to the 4+ ground state, that would correspond 
to K=4+ in a Nilsson context (the aligned coupling of 
the last unpaired proton K=3/2~ and the last neutron 
K=5/2~) and a second one linked to the 1+ member of 
the ground state multiplet that results of the anti-aligned 
coupling. 

The negative parity states can be interpreted as the 
result of the coupling of a hole in the ld3/2 orbit to ^^Cr. 
As we have mentioned earlier, the ground state band of 
■^^Cr can be understood in a particle plus rotor picture, 
and assigned K=5/2~. Now we have to couple the ld3/2 
hole to our basic rotor ^^Cr. Therefore the structure of 
the negative parity states of ''^V is one particle and one 
hole coupled to the rotor core of **Cr. Hence we expect 
that the lowest lying bands would have K=l~ and K=4~ 
(aligned and anti-aligned coupling of the quasiparticle 
and the quasihole). These two bands have been found 
experimentally at about 500 kcV of excitation energy. 
Their structure is in very good agreement with the SM 
predictions (see Fig. 44) except for a small shift of the 
K=4~ and K=8~ bandheads relative to the K=l~ state. 

A more stringent test of the theoretical predictions is 
provided by the electromagnetic properties. In this par- 
ticular nucleus the experimental information is very rich 
and a detailed comparison can be made for both the pos- 
itive and negative bands and for E2 and Ml transitions. 
The results for the ground state band are collected in 



FIG. 44 The negative parity K=l- and K=4- bands of **V, 
experiment vs. theory. 




FIG. 45 Theoretical and experimental B(M1) and B(E2) 
transition probabilities in the ground state band of 



Fig. 45. Notice that the calculation reproduces even the 
tiniest details of the experimental results, as, for instance 

the odd-even staggering of the B(Ml)'s. The B(E2)'s cor- 
respond to a deformation /3=0.2 that is clearly smaller 
than that of *^Cr. 

The transition probabilities in the negative parity 
bands are presented in Fig. 46. Now the experimen- 
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FIG. 46 Theoretical and experimental B(M1) and B(E2) 
transition probabilities in the negative parity bands of 



this monopolc drift provides direct evidence for the need 
of thrcc-body mechanisms. 

For the p shell, the situation is similar: The imposing 
11.5 MeV ^1/2-^5/2 ga-P in is down to some 3 MeV 
in ^^C. The monopole drift of the 51/2-^5/2 gap brings 
it from about 8 MeV in ^^Si to nearly -1 MeV in ^^C. 
According to Otsuka et al. (2001a) the drift may well 
continue: in ^He, the S1/2-P1/2 bare gap is estimated at 

0. 8 MeV. 

All the numbers above (except the last) are experi- 
mental. The bare monopole values are smaller because 
correlations increase substantially the value of the gaps, 
but they do not change qualitatively the strong monopole 
drifts. Their main consequence is that "normal" states, 

1. e., those described by Qhw p or sd calculations often 
"coexist" with "intruders" that involve promotion to the 
next oscillator shell. Let us examine how this happens. 

A. N=8; "Li: halos 



tal results have larger uncertainties. Despite that, the 

agreement in absolute values and trends is very good for 
the lowest band (K=l^). The quantitative agreement is 
worse for the other two bands, nevertheless, the main fea- 
tures of the data are well reproduced. The deformation 
of the K=l~ band extracted from the B(E2) properties 
turns out to be larger than the deformation of the ground 
state band and close to that of ^^Cr. 



VII. DESCRIPTION OF VERY NEUTRON RICH NUCLEI 

The study of nuclei lying far from the valley of stabil- 
ity is one of the most active fields in today's experimen- 
tal nuclear physics. What it spec;ific; about these nuclei, 
from a shell model point of view? As everywhere else: 
the model space and the monopole behavior of the in- 
teraction. They go together because the effective single 
particle energies (ESPE, see Section V) depend on occu- 
pancies, which in turn depend on the model space. 

In broad terms, the specificity of light and medium 
neutron rich nuclei is that the EI closures (correspond- 
ing to the filling of the p -I- 1/2 orbit of each oscillator 
shell) take over as boundaries of the model spaces. As 
discussed in Section II.B.2, the oscillator closures (HO) 
may be quite solid for double magic nuclei, but they be- 
come vulnerable in the semi-magic cases. 

For instance the d3/2-,/7/2 neutron gap in *°Ca of ~ 
7 MeV goes down to « 2.5 MeV around ^^O. Since the 
orbit is well below its sd partners, now quite close 
to /7/2, the natural model space is no longer the sd shell, 
but the EI space bounded by the = 14 and 28 closures, 
supplemented by the P3/2 subshell whenever the P3/2-/7/2 
gap becomes small: from « 6 MeV in ^^Ni, it drops to « 
4.5 MeV in "SQa, then « 2 MeV in ^"Ca, and finally w 
MeV in ^^Si. As explained at the end of Section II.B.2 



The 1/2+ ground state in ^^Be provided one of the first 
examples of intrusion. The expected Ohiu normal state 
lies 300 keV higher. The explanation of this behavior 
has varied with time (Aumann et al., 2001; Sagawa et al, 
1993; Suzuki and Otsuka, 1994; Talmi and Unna, 1960), 
see also (Brown, 2001) for a recent review and (Suzuki 
et al., 2003) for a new multi-?KJ calculation, but the idea 
has remained unchanged. The normal state corresponds 
to a hole on the N = 8 closure. The monopole loss of 
promoting a particle to the sd shell is compensated by a 
pairing gain for the p~^^ holes. A quadrupole gain due to 

the interaction of the sd neutron with the protons is 
plausible, but becomes questionable when we note that 
the same phenomenon occurs in ^He which has no 
particles, suggesting the need for a further reduction of 
the monopole loss (Otsuka et al., 2001a). 

The interest in ^^Be — as a "halo" nucleus — was revived 
by the discovery of the remarkable properties of ^^Li, 
which sits at the drip line (» 200 keV two-neutron sepa- 
ration energy) and has a very large spatial extension, due 
to a neutron halo (Hansen and Jonson, 1987; Tanihata 
et al., 1985). 

The shell model has no particular problem with halo 
nuclei, whose large size is readily attributed to the large 
size of the si/2 orbit. As shown by (Kahana et al., 1969a) 
the use of Woods Saxon wavefunctions affect the matrix 
elements involving this orbit, but the uncertainties in- 
volved are easily absorbed by the monopole field. As a 
consequence, the whole issue hinges on the Si/2 contribu- 
tion to the wavefunctions, which is sensitively detected 
by the (3 decay to the first excited l/2~ state in ^^Be. In 
Borge et al. (1997) it was shown that calculations produc- 
ing a 50% split between the neutron closed shell and the 
■^1/2^1/2 configuration lead to the right lifetime. This 
result was confirmed by Simon et al. (1999). Further 
confirmation comes from Navin et al. (2000), with solid 
indications that the supposedly semi-magic ^^Be ground 
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state is dominated by the same s^j^ P1/2 configuration. 



B. N=20; ^^Mg, deformed intruders 

In the mid 1970's it was the sd shell that attracted the 
most attention. Nobody seemed to remember ^^Be, and 
everybody (including the authors of this review active 
at the time) were enormously surprised when a classic 
experiment Thibault et al. (1975) established that the 
mass and /3-decay properties of the ^^Na ground state — 
expected to be semimagic at = 20 — could not possi- 
bly be that of a normal state (Wildenthal and Chung, 
1979). The next example of a frustrated semi- magic was 
•^^Mg (Detraz et al, 1979). Early mean field calcula- 
tions had interpreted the discrepancies as due to defor- 
mation (Campi et al., 1975) but the experimental confir- 
mation took some time (Guillemaud-Mueller et al., 1984; 
Klotz et al, 1993; Motobayashi et al, 1995). 

Exploratory shell model calculations by Storm et al. 
(1983), including the lf7/2 orbit in the valence space, 
was able to improve the mass predictions, however, de- 
formation was still absent. To obtain deformed solutions 
demanded the inclusion of the 2p3/2 orbit as demon- 
strated by Poves and Retamosa (1987). These calcu- 
lations were followed by many others (Fukunishi et al., 
1992; Heyde and Woods, 1991; Otsuka and Fukunishi, 
1996; Poves and Retamosa, 1994; Siiskonen et al, 1999; 
Warburton et al, 1990), that mapped an "island of inver- 
sion" , i. e. the region where the intruder configurations 
are dominant in the ground states. The detailed con- 
tour of this "island" depends strongly on the behavior 
of the effective single particle energies (ESPE), in turn 
dictated by the monopole hamiltonian. Let's recall, that 
according to what we have learned in Section VLB, the 
configurations {d^/2Si/2)l~^ (/7/2P3/2)^~^,corresponding 
to iV = 20, Z = 10-12, have a "quasi-SU(3)" quadrupole 
coherence close to that of SU(3) {i.e., maximal). 

The ESPE in Fig. 47 represent Hm for the sdfp.sm and 
Tokyo group interactions (Utsuno et al, 1999). Within 
details they are quite close, and such as efhciently to 
favor deformation. Both interactions lead to an "island 
of inversion" for Z=10, 11 and 12; 7V=19, 20 and 21. 

Consider now some detailed information obtained with 
sdpf.sm by Caurier et al (2001c). The iS'2Ar values of 
Fig. 48 locate the neutron drip line, consistent with what 
is known for oxygen and fluorine, where the last bound 
isotopes are ^''O and ^^F (Sakurai et al, 1999). Note 
the kink due to deformed correlations in the latter. For 
the other chains, the behavior is smoother and the last 
predicted bound isotopes are ■^'^Ne, ■^^Na and ^'^Mg. 

Some results for the even Mg isotopes (N=18, N=20 
and N=:22) are gathered in table XII. In '^"Mg the normal 
configuration is the one that agrees with the existing ex- 
perimental data (Pritychenko et al, 1999). In '^^Mg the 
situation is the opposite as the experimental data (the 2+ 
excitation energy (Guillemaud-Mueller et al, 1984) and 
the 0+ 2+ B(E2) (Motobayashi et al, 1995)) clearly 




FIG. 47 Effective single particle energies (in MeV) at N=20 
with the sdpf.sm interaction (upper panel) and the Tokyo 
group interaction (bottom panel) 
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FIG. 48 Two neutron separation energies (in MeV) calculated 
with the sdpf.sm interaction. 



prefer the intruder. A preliminary measure of the 4+ 
excitation energy reported by Azaiez (1999) goes in the 
same direction. Data and calculations suggest prolate 
deformation with f3 « 0.5. In "^'^Mg the normal config- 
uration that contains two pf neutrons, is already quite 
collective. It can be seen in the table that it resembles 
the '^^Mg ground state. The 4p-2h intruder is even more 
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TABLE XII Properties of the even magnesium isotopes. N stands for normal and I for intruder. Energies in MeV, B(E2)'s in 
e^fm'' and Q's in efm^ 







Mg 






Mg 






Mg 






N 


I 


EXP 


N 


I 


EXP 


N 


I 


EXP 


AE(0+) 




+3.1 






-1.4 






+ 1.1 




0+ 


0.0 


0.0 




0.0 


0.0 




0.0 


0.0 




2+ 


1.69 


0.88 


1.48 


1.69 


0.93 


0.89 


1.09 


0.66 


0.67 


4+ 


4.01 


2.27 




2.93 


2.33 


(2.29) 


2.41 


1.86 


2.13 


6+ 


6.82 


3.75 




9.98 


3.81 




3.52 


3.50 




B(E2) 




















2+ ^0+ 


53 


112 


59(5) 


36 


98 


90(16) 


75 


131 


126(25) 


4+ ^ 2+ 


35 


144 




17 


123 




88 


175 




6+^4+ 


23 


140 




2 


115 




76 


176 




Q spec (2"'") 


-12.4 


-19.9 




-11.4 


-18.1 




-15.4 


-22.7 





deformed (/? « 0.6) and a better rotor. Results from the 
Riken experiments of Yoneda et al. (2000) and Iwasaki 
et al. (2001) seem to favor the intruder option. In the 
QMCD calculations of Utsuno et al. (1999), the groimd 
state band is dominantly 4p-2h; the 2+ comes at the right 
place, but the 4+ is too high. Clearly, ^^Mg is at the edge 
of the "island of inversion" . Another manifestation of the 
intruder presence in the region has been found at Isolde 
(Nummela et al., 2001b): The decay of '^'^Na indicates 
that the ground state of ^^Mg has J'^=3/2+ instead of 
the expected J'^=3/2~ or J'^=7/2~. This inversion is 
nicely reproduced by the sdpf.sm calculation. 



C. N=28: Vulnerability 

Let us return for a while on Fig. 47. The scale does 
not do justice to a fiuidamental feature: the drift of the 
P3/2 — fiji gaP) which decreases as protons are removed. 
As explained at the end of Section II. B. 2, this behavior 
is contrary to a very general trend in heavier nuclei, and 
demands a three-body mechanism to resolve the contra- 
diction. Next we remember that the oscillator closures 
arc quite vulnerable even at the strict monopolc level. As 
we have seen, quadrupole coherence takes full advantage 
of this vulnerability. On the contrary, the EI closures 
are very robust. However, because of the drift of the 
P3/2 — /7/2 gap even the iV = 28 closure becomes vul- 
nerable. The sdpf.sm interaction leads to a remarkable 
result summarized in Table XIII: the decrease of the gap 
combined with the gain of correlation energy of the 2p-2h 
lead to a breakdown of the = 28 closure for '*^S and 
40Mg. The double magic 42sij.ggig^s^ Towards N=Z, ^^Cr 

and ^^Fe exhibit large correlation energies associated to 
the prolate deformed character of their 2p-2h neutron 
configurations (/? ~ 0.3). 

When the 2p-2h bandheads of *°Mg and ^^Cr are al- 
lowed to mix in the full Ohiu space, using them as pivots in 
the LSF procedure, they keep their identity to a large ex- 



tent. This can be seen in Fig. 49, where we have plotted 
the strength functions of the 2p-2h 0+ states in the full 
space. In ^°Mg the 2p-2h state represents the 60% of the 
ground state while in ''^Cr it is the dominant component 
(70%) of the first excited 0+. This is a very interesting 
illustration of the mechanism of intrusion; the intruder 
state is present in both nuclei, but it is only in the very 
neutron rich one that it becomes the ground state. 
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FIG. 49 LSF of the 2p-2h 0+ bandheads of ''"Mg (upper 
panel) and ^^Cr (bottom panel) in the full OHlj space 

Table XIV gives an idea of the properties of the iso- 
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TABLE XIII N=28 isotones: quasiparticle neutron gaps, dif- 
ference in correlation energies between the 2p-2h and the Op- 
Oh configurations and their relative position 





«Mg 


«Si 


44s 


4«Ar 


4«Ca 




52 Cr 


s^Fe 


56Ni 


gap 


3.35 


3.50 


3.23 


3.84 


4.73 


5.33 


5.92 


6.40 


7.12 


AEcorr 


8.45 


6.0 


6.66 


5.98 


4.08 


7.59 


10.34 


10.41 


6.19 


^2p-2h 


-1.75 


1.0 


-0.2 


1.7 


5.38 


3.07 


1.50 


2.39 


8.05 



tones in which configuration mixing is appreciable, but no 
clearcut deviation from N=28 magicity can be detected. 
In Fig. 50 we have plotted the low energy spectra of the 
heaviest known sulphur isotopes. The agreement with 
the accumulated experimental results (sec Glasmacher, 
1998, for a recent review) is excellent and extends to the 
new data of Sohler et al. (2002). Analyzing their proton 
occupancies we conclude that the rise in collectivity along 
the chain is correlated with the equal filling of the 
and Si/2 orbitals, (the ^5/2 orbital remains always nearly 
closed). The maximum proton collectivity is achieved 
when both orbitals are degenerate, which corresponds to 
the Pseudo-SU(3) limit. For the neutrons, the maximum 
collectivity occurs at A'' = 24, the /7/2 mid-shell. Ac- 
cording to the calculation, '*^S is a prolate rotor, with 
an incipient 7 band. In ^^S the spherical and deformed 
configurations mix equally. 

The N=27 isotones also reflect the regular transition 
from sphericity to deformation in their low-lying spec- 
trum. The excitation energy of the 3/2~ state should 
be sensitive to the correlations and to the neutron gap. 
While in *'^Ca, it lies quite high (at around 2 MeV) 
due to the strong 77/2 closure, in ^^Ar it appears at 
about 0.4 MeV. Concerning "^^S, the information comes 
from two recent experiments: the mass measurements at 
Ganil by Sarazin et al. (2000) that observed a low-lying 
isomer around 400 keV excitation energy and the MSU 
Coulex experiment of Ibbotson et al. (1999) that detected 
a strong E2 transition from the groimd state to an excited 
state around 940 keV. According to our calculations the 
ground state corresponds to the deformed configuration 
and has spin 3/2~. The spherical single hole state 7/2~ 
would be the first excited state and his lifetime is con- 
sistent with that of the experimental isomer. The third 
known state is short lived and it should correspond to 
the 7/2~ member of the ground state band. 



D. IM=40; from "magic" ^*Ni to deformed ^''Cr 

Systematics is a most useful guide, but it has to be 
used with care. The p and sd shells can be said to be 
"full" Qhhj spaces, in the sense that the region bound- 
aries are well defined by the iV = 4, 8 and 20 closures, 
with the exception for the very neutron-rich halo nuclei 
and the island of inversion discussed previously. In the 
p/-shell it is already at = Z « 36 that the Qhuj model 



TABLE XIV N=28 isotones: Spectra, quadrupole properties 
and occupancies 





«Mg 


*2Si 


44g 


4«Ar 


E*(2+)(MeV) 


0.81 


1.49 


1.22 


1.51 


E*(4+) 


2.17 


2.68 


2.25 


3.46 


E*(0+) 


1.83 


1.57 


1.26 


2.93 


Q(2+)(c fm^) 


-21 


16 


-17 


20 


B(E2)(e2 fm'') 


108 


71 


93 


93 


{^7/2) 


5.54 


6.16 


6.16 


6.91 


{f7/2f % 


3 


28 


24 


45 



TABLE XV 2f energies and gg/2 intruder occupation in the 

Nickel isotopic chain, from (Sorlin et al., 2002) 

«2Ni 64Ni ee^j es^j ro^j 72-^- 74^^ 

E(2+)caic. 1-11 1-24 1.49 1.73 1.50 1.42 1.33 

E(2+)ea,p. 1.173 1.346 1.425 2.033 1.259 

B(E2t)caic. 775 755 520 265 410 505 690 

(rig/a) 0.24 0.43 0.67 1.07 0.84 0.55 0.45 



space collapses under the invasion of deformed intruders. 

Naturally, we expect something similar in the neighbor- 
hood of iV = 40 isotones. Recent experiments involving 
the Coulomb excitation of ^e-esj^j (Sorlin et al., 2002), 
the /? decay of 60-63 V (SorHn et al, 2003), and the de- 
cay and the spectroscopy of the 67nipg (Sawicka et al, 
2003), support this hypothesis. In addition, the (3 decays 
of the neutron-rich isotopes 6'*Mn and 66 Mn, investigated 
at Isolde (Hannawald et al, 1999), that show a sudden 
drop of the 2+ energies in the daughter nuclei ^^Fe and 
^^Fe, point also in the same direction. 

In defining the new model spaces, the full pf results are 
quite useful: "''^Ni turns out to be an extraordinarily good 
pivot (refer to Fig. 15), and as soon as a few particles 
are added it becomes a good core. Therefore, for the 
first time we are faced with a genuine EI space 
(remember ra is the "rest" of the p — 'i shell). Its validity 
is restricted to nearly spherical states. Deformation will 
demand the addition of the c?5/2 subshell. 

The effective interaction is based on three blocks; i) 
the TBME from KB3G effective interaction (Poves et al, 
2001), ii) the G-matrix of ref . (Hjorth-Jensen et al., 1995) 
with the modifications of ref. (Nowacki, 1996) and iii) the 
Kahana, Lee and Scott G-matrix (Kahana et al., 1969b) 
for the remaining matrix elements. We have computed 
all the nickel isotopes and followed their behavior at and 
beyond N=40. Let's first focus on 6SNi. The calculated 
low energy spectrum, an excited 0"*" at 1.85 MeV, 2^ at 
2.15 MeV, 5- at 2.72 MeV and 4+ at 3.10 MeV, fit re- 
markably well with the experimental results. The ground 
state wave function is 50% closed shell. This number 
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FIG. 50 Predicted level schemes of the heavy sulphur isotopes, compared with the experiment. Energies in keV, B(E2)'s in 

e^ fm* 



is very sensitive to modifications of the pf-g gap, and 
smaller values of the closed shell probability can be ob- 
tained without altering drastically the rest of the prop- 
erties. The SM results have been compared with those of 
other methods in (Langanke et al., 2003b). In table XV 



TABLE XVI Spectroscopic properties of ^o-e^Cr in the fpgd 
valence space 





80 Cr 


82Cr 


s^Cr 


E*(2+) (MeV) 


0.67 


0.65 


0.51 


Qs(e.fm2) 


-23 


-27 


-31 


B(E2)i(e2.fm4) 


288 


302 


318 


Qi(e.fm^) from Qs 


82 


76 


109 


Qi(e.fm^) from B(E2) 


101 


103 


106 


E*(4+) (MeV) 


1.43 


1.35 


1.15 


Q.(e.fm2) 


-37 


-30 


-43 


B(E2)i(e2.fm'') 


426 


428 


471 


Qi(e.fm^) from Qs 


102 


84 


119 


Qi(e.fm^) from B(E2) 


117 


117 


123 



we compare the excitation energy of the 2+ states in the 
Nickel chain with the available experimental data (in- 
cluding the very recent ^"Ni value (Sorlin et al., 2002)). 
The agreement is quite good, although the experiment 
gives a larger peak at N=40. Similarly the B(E2)'s in 
Fig. 51 follow the experimental trends including the drop 
at ^*Ni recently measured (Sorlin et al., 2002). Note that 
for a strict closure the transition would vanish. 

The /3 decays of ^"V and ^^V indicate very low en- 
ergies for the 2^ states in ^"Cr and ^^Cr (Sorlin et al.. 
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FIG. 51 B(E2) (0+ 2+) in e^ fm" for the Nickel isotopes, 
from (SorUn et al., 2002) 

2000, 2003). If N=40 were a magic number, one would 

expect higher 2+ energies as the shell closure approaches, 
particularly in '^"'Cr. As seen in Fig. 52, exactly the oppo- 
site seems to happen experimentally. The three possible 
model spaces tell an interesting story: pf = rs is accept- 
able at A'' = 32, 34. The addition of 99/2 does more harm 
than good because the gaps have been arbitrarily reduced 
to reproduce the A" = 36 experimental point by allowing 
2p-2h jumps including ^5/2- But then, in A/' = 38 the 
2+ is too high, and the experimental trend promises no 
improvement at A^ = 40. 

Properties of the 2+ energies given by the pfgd calcu- 
lation are collected in table XVI. They point to prolate 
structures with deformation /3 « 0.3. 
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FIG. 52 Experimental and calculated 2"*" energies in the pf, 
pfg and pfgd spaces. 

VIII. OTHER REGIONS AND THEMES 
A. Astrophysical applications 

Astrophysical environments involve conditions of tem- 
perature and densities that are normally not accessible in 
laboratory experiments. The description of the different 
nuclear processes requires then theoretical estimates. As 
discussed in previous sections, shell-model calculations 
are able to satisfactorily reproduce many experimental 
results so that it should be possible to obtain reliable 
predictions for nuclei and/or conditions not yet accessi- 
ble experimentally. 

As a first example, consider the decay properties of 
nuclei totally stripped of the atomic electrons at typi- 
cal cosmic rays energies of 300 MeV/A. A nucleus such 
as ^^Mn, unstable in normal conditions, becomes stable. 
Another isotope, ^'*Mn, could be used as a chronome- 
ter to study the propagation of iron group nuclei on 
cosmic-rays (Duvernois, 1997) provided its half-life un- 
der cosmic-rays conditions is known. The necessary cal- 
culation involves decay by second-forbidden unique tran- 
sitions to the ground states of ^^Cr and ^''Fe. Due to 
phase space arguments the (3~ decay to ^^Fe is expected 
to dominate. In two difficult and elegant experiments 
the very small branching ratio for the /S"*" decay to the 
ground state of ^'^Cr has been measured: (1.8±0.8) x 10~^ 
by Zaerpoor et al. (1997) and (1.20±0.26) x 10"^ by Wu- 
osmaa et al. (1998). Taking the weighted mean of this 
values, (1.26 ± 0.25) x 10^^, and knowing that of ^-^Mn, 
312.3(4) d, a partial /?+ half-life of (6.8 ± 1.3) x 10^ yr 
is obtained compared with 5.6 x 10^ yr from the shell 
model calculation of Martmez-Pinedo and Vogel (1998), 
which yields 5.0 x 10^ yr for the dominant /3~ branch. 
Both theoretical results are sensitive to uncertainties in 
the renormalization of the unique second-forbidden oper- 
ators that should be removed taking the ratio of the /3~ 
and half-lives. Multiplying the theoretical ratio by 
the experimental value of the /3+ half-life yields a value 
of (6.0 ± 1.2) X 10^ yr for the (3~ branch. Using a similar 



argument Wuosmaa et al. (1998) estimate the partial f3 
half-life to be (6.0± 1.3[stat] ± l.l[theor]) x 10^. The in- 
fluence of this half-life on the age of galactic cosmic-rays 
has been discussed recently by Yanasak et al. (2001). An- 
other possible cosmic ray chronometer, ^^Ni, could mea- 
sure the time between production of iron group nuclei in 
supernovae and the accelaration of part of this matterial 
to form cosmic rays (Fisker et ai, 1999). Before acceler- 
ation, the decay of ^^Ni proceeds by electron capture to 
the 1+ state in ^^Co with a half-live of 6.075(20) days. 
After acceleration ^^Ni is stripped of its electrons, the 
transition to the 1+ state is no longer energetically al- 
lowed and the decay proceeds to the 3+ state at 158 keV 
via a second forbidden unique transition. Currently only 
a lower limit for the half-life of totally ionized ^^Ni could 
be estabhshed (2.9 x 10"* y) (Sur et ai, 1990). A recent 
shell-model calculation by (Fisker et ai, 1999) predict a 
half-live of 4 x 10"* y that is too short for ^^Ni to serve as 
a cosmic ray chronometer. 

Nuclear beta-decay and electron capture are important 
during the late stages of stellar evolution (see Langanke 
and Martmez-Pinedo, 2003, for a recent review). At the 
relevant conditions in the star electron capture and P de- 
cay are dominated by Gamow- Teller (and Fermi) transi- 
tions. Earlier determinations of the appropriate weak in- 
teraction rates were based in the phenomenological work 
of Fuller, Fowler, and Newman (1980, 1982a,b, 1985). 
The shell model makes it possible to refine these esti- 
mates. For the sd shell nuclei, important in stellar oxygen 
and silicon burning, we refer to Oda et al. (1994). More 
recently, it has been possible to extend these studies to pf 
shell nuclei relevant for the pre-supernova evolution and 
collapse (Caurier et ai, 1999a; Langanke and Martmez- 
Pinedo, 2000; Langanke and Martmez-Pinedo, 2001). 
The astrophysical impact of the shell-model based weak 
interaction rates have been recently studied by Heger 
et al. (2001a,b) 

The basic ingredient in the calculation of the differ- 
ent weak interaction rates is the Gamow- Teller strength 
distribution. The GT^ sector directly determines the 
electron capture rate and also contributes to the beta- 
decay rate through the thermal population of excited 
states (Fuller et al, 1982a). The GT_ strength con- 
tributes to the determination of the /3-decay rate. To 
be applicable to calculating stellar weak interaction rates 
the shell-model calculations should reproduce the avail- 
able GT+ (measured by (n,p)-type reactions) and GT_ 
(measured in (p, n)-type reactions). Figure 53 com- 
pares the shell-model GT-f distributions with the pi- 
oneering measurements performed at TRIUMF. These 
measurements had a typical energy resolution of sa 
1 MeV. Recently developed techniques, involving (d, ^He) 
charge-exchange reactions at intermediate energies (Rak- 
ers et ai, 2002), have improved the energy resolution 
by an order of magnitude or more. Figure 54 compares 
the shell-model GT-|_ distribution computed using the 
KB3G interaction (Poves et al., 2001) with a recent ex- 
perimental measurement of the ^^Y{d, ^He) performed at 
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FIG. 53 Comparison of shell model GT+ distributions with 
experimental data (Alford et al, 1993; El-Kateb et ai, 1994; 
Ronnqvist et al, 1993) for selected nuclei. The shell model 
results (discrete lines) have been folded with the experimen- 
tal resolution (histograms) and include a quenching factor of 
(0.74)^ (adapted from Caurier et al, 1999a). 



KVI (Baumer et al, 2003). 

Figure 55 compares the shell-model GT_ distributions 
with the data obtained in (p, n) charge-exchange reaction 
measurements for 54,56pg g^j^j 58,60]\jj (Anderson et ai, 
1990; Rapaport et ai, 1983). The GT_ operator act- 
ing on a nucleus with neutron excess and ground state 
isospin T can lead to states in the daughter nucleus with 
three different isospin values (T— 1, T, T+1). As a conse- 
quence, the GT_ strength distributions have significantly 
more structure and extend over a larger excitation energy 
interval than the GT-|_ distributions, making their theo- 
retical reproduction more challenging. Nevertheless, the 
agreement with the experimental data is quite satisfac- 
tory. The shell-model results for ^*Ni have been recently 
compared with high-resolution data (50 keV) obtained 
using the (■^He,t) reaction (Fujita et ai, 2002). 

Shell-model diagonalization techniques have been used 
to determine astrophysically relevant weak interaction 
rates for nuclei with ^ < 65. Nuclei with higher masses 
are relevant to study the collapse phase of core-collapse 
supernovae (Langanke and Martmez-Pinedo, 2003). The 
calculation of the relevant electron-capture rates is cur- 
rently beyond the possibilities of shell-model diagonal- 
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ization calculations due to the enormous dimensions of 
the valence space. However, this dimensionality prob- 
lem does not apply to Shell-Model Monte-Carlo meth- 
ods (SMMC, see section IV.B.3). Moreover, the high 
temperatures present in the astrophysical environment 
makes necessary a finite temperature treatment of the 
nucleus, this makes SMMG methods the natural choice 
for this type of calculations. Initial studies by Langanke 
et al. (2001) showed that the combined effect of nuclear 
correlations and finite temperature was rather efficient 
in unblocking Gamow- Teller transitions on neutron rich 
germanium isotopes. More recently this calculations have 
been extended to cover all the relevant nuclei in the range 
A = 65-112 by Langanke et al. (2003a). The resulting 
electron-capture rates have a very strong influence in the 
collapse (Langanke et ai, 2003a) and post-bounce (Hix 
et ai, 2003) 

The astrophysical r-process is responsible for the syn- 
thesis of at least half of the elements heavier than A « 
60 (Wallerstein et al., 1997). Simulations of the r-process 
require the knowledge of nuclear properties far from the 
valley of stability (Kratz et ai, 1998; Pfeiffer et ai, 2001). 
As the relevant nuclei are not experimentally accessible, 
theoretical predictions for the relevant quantities (i.e., 
neutron separation energies and half-lives) are needed. 
The calculation of (3 decay half-lives usually requires two 
ingredients: the Gamow- Teller strength distribution in 
the daughter nucleus and the relative energy scale be- 
tween parent and daughter (i.e. the Q/3 value). Due 
to the huge number of nuclei relevant for the r pro- 



49 




5 10 15 

E (MeV) 



5 10 15 

E (MeV) 



FIG. 55 Comparison of {p, n) L — forward angle cross sec- 
tion data (Rapaport et al, 1983) (upper panels) with the cal- 
culated GT_ strength distributions (lower panels). For ^''Fe 
the solid curve in the lower panel shows the experimental 
GT_ data from (Anderson et al., 1990), while the shell model 
results are given by the dashed line. The Fermi transition to 
the isobaric analog state (IAS) is included in the data (upper 
panels) but not in the calculations. A quenching factor of 
(0.74)'^ is included in the calculations, (from Caurier et al., 
1999a). 



cess, the estimates of the half-lives are so far based 
on a combination of global mass models and the quasi 
particle random-phase approximation (see Langanke and 
Marti'nez-Pincdo, 2003, for a description of the different 
models). However, recently shell- model calculations have 
become available for some key nuclei with a magic neu- 
tron number TV = 50 (Langanke and Marti'nez-Pinedo, 
2003), = 82 (Brown et al, 2003; Marti'ncz-Pincdo and 
Langanke, 1999), and N = 126 (Martmez-Pinedo, 2001). 
All this calculations suffer from the lack of spectroscopic 
information on the regions of interest that is necessary 
to fine tune the effective interactions. This situation is 
improving at least for = 82 thanks to the recent spec- 
troscopic data on ^^°Cd (Dillmann et al., 2003). 

Nuclear reaction rates are the key input data for sim- 
ulations of stellar burning processes. Experiment-based 
reaction rates for the simulation of explosive processes 
such as novae, supernovae, x-ray bursts, x-ray pulsars, 
and merging neutron stars are scarce because of the ex- 
perimental difficulties associated with radioactive beam 



measurements (Kappeler et al., 1998). Most of the re- 
action rate tables are therefore based on global model 
predictions. The most frequently used model is the sta- 
tistical Hauser Feshbach approach (Rauscher and Thiele- 
mann, 2000). For nuclei near the drip-lines or near closed 
shell configurations, the density of levels is not high 
enough for the Hauser Feshbach approach to be appli- 
cable. For these cases alternative theoretical approaches 
such as the nuclear shell model need to be applied. Shell- 
model calculations were used for the determination of the 
relevant proton capture reaction rates for sd-shell nu- 
clei necessary for rp process studies (Herndl et al., 1995). 
This calculations have been recently extended to include 
p/-shell nuclei (Fisker et al, 2001). 

Knowledge of neutrino nucleus reactions is necessary 
for many applications, e.g. the neutrino oscillation stud- 
ies, detection of supernova neutrinos, description of the 
neutrino transport in supernovae and nucleosynthesis 
studies. Most of the relevant neutrino reactions have 
not been studied experimentally so far and their cross 
sections are typically based on nuclear theory (see Kolbe 
et al., 2003, for a recent review). The model of choice 
for the theoretical description of neutrino reactions de- 
pends of the energy of the neutrinos that participate in 
the reaction. 

For low neutrino energies, comparable to the nuclear 
excitation energy, neutrino-nucleus reactions are very 
sensitive to the appropriate description of the nuclear 
response that is very sensitive to correlations among nu- 
cleons. The model of choice is then the nuclear shell- 
model. Qhid calculations have been used for the calcula- 
tion of neutrino absorption cross sections (Sampaio et al., 
2001) and scattering cross sections (Sampaio et al., 2002) 
for selected pf shell nuclei relevant for supernovae evo- 
lution. For lighter nuclei complete diagonalizations can 
be performed in larger model spaces, e.g. Afiio calcula- 
tions for ^^O (Haxton, 1987; Haxton and Johnson, 1990) 
and Qfiio calculations for ^^C (Hayes and Towner, 2000; 
Volpe et al., 2000). Other examples of shell- model cal- 
culations of neutrino cross sections are the neutrino ab- 
sorption cross sections on ^°Ar of Ormand et al. (1995) 
for solar neutrinos (see Bhattacharya et al., 1998, for 
an experimental evaluation of the same cross section), 
this cross section have been recently evaluated by (Kolbe 
et al., 2003) for supernova neutrinos. And the evaluation 
by Haxton (1998) of the solar neutrino absorption cross 
section on ''^Ga relevant for the GALLEX and SAGE 
solar neutrino experiments. 

For higher neutrino energies the standard method of 
choice is the random phase approximation as the neutrino 
reactions are sensitive mainly to the total strength and 
energy centroids of the different multipoles contributing 
to the cross section. In some selected cases, the Fermi and 
Gamow- Teller contribution to the cross section could be 
determined from shell-model calculation that is supple- 
mented by RPA calculations for higher multipoles. This 
type mix calculation has been carried out for several iron 
isotopes (Kolbe et al., 1999; Toivanen et al., 2001) and 
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for 20Ne (Heger et al., 2003). 



B. /3/3-decays 

The double beta decay is the rarest nuclear weak pro- 
cess. It takes place between two even-even isobars, when 
the decay to the intermediate nucleus is energetically for- 
bidden or hindered by the large spin difference between 
the parent ground state and the available states in the 
intermediate nuclei. It comes in three forms: The two- 
neutrino decay /3(i2v- 



is just a second order process mediated by the Standard 
model weak interaction It conserves the lepton number 
and has been already observed in a few nuclei. 
The second mode, the neutrinoless decay fi^ov'- 



^Z+2 ^N- 



+ 61 +eo 



needs an extension of the standard model of electroweak 
interactions as it violates lepton number. A third mode, 
PPov,x is also possible 



N 



*Z+2 



Xjv-2 + ei +62 +x 



in some extensions of the standard model and proceeds 
via emission of a light neutral boson, a Majoron x- The 
last two modes, not yet experimentally observed, require 
massive neutrinos an issue already settled by the recent 
measures by Super-Kamiokande (Fukuda et al., 1998), 
SNO (Ahmad et al, 2002) and KamLAND (Eguchi et al, 
2003) . Interestingly, the double beta decay without emis- 
sion of neutrinos would be the only way to sign the Majo- 
rana character of the neutrino and to distinguish between 
the different scenarios for the neutrino mass differences. 
Experimentally, the three modes show different electron 
energy spectra ((see figure 2 in Zdesenko, 2002)), the 
f}/32v and PPov,x are characterized by a continuous spec- 
trum ending at the maximum available energy Q/3/3, while 
the (3f3oiy spectrum consists in a sharp peak at the end 
of the Qfjjs spectrum. This should, in principle, make 
it easier the signature of this mode. In what follows we 
shall concentrate in the /3/32;y and /3/3o;y modes. 

The theoretical expression of the half-life of the 2v 
mode can be written as: 



1%] ^ = <^2i/|AfQ^^|^, 



with 



(J+||^t-||l+)(l+||^t-||0+) 

Em + Eq{J) 



(94) 



(95) 



(there is an implicit sum over all the nucleons). G2V 
contains the phase space factors and the axial coupling 

constant qa- The calculation of M'qj, requires the precise 
knowledge of the ground state of the parent nuclei and 



the ground state and occasionally a few excited states of 
the grand daughter, both even-even. Besides, it is nec- 
essary to have a good description of the Gamow- Teller 
strength functions of both of them, which implies a de- 
tailed description of the odd-odd intermediate nucleus. 
This is why this calculation is a challenge for the nuclear 
models, and why agreement with the experiment in this 
channel is taken as a quality factor to be applied to the 
predictions of the models for the neutrinoless mode. 

It is also a show-case of the use of the Lanczos strength 
function (LSF) method. It works as follows: Once the 
relevant wave functions of parent and grand daughter 
are obtained, it is straightforward to build the doorway 
states at-\Qj^^^^J and Bt+\Jf^^J. In a second step, 
one of them is fragmented using LSF, producing, at it- 
eration N, N 1"*" states in the intermediate nucleus, with 
excitation energies E^. Overlapping these vectors with 
the other doorway, entering the appropriate energy de- 
nominators and adding up the N contributions gives an 
approximation to the exact value of M'qj, (N=1 is just 
the closure approximation). Finally, the number of iter- 
ations is increased until full convergence is reached. The 
method is very efficient, for instance, in the A=48 case, 
20 iterations suffice largely. The contributions of the dif- 
ferent intermediate states to the final matrix element are 
plotted in Fig. 56. 



0.01 



5x10 




Excitation Energy(MeV) 

FIG. 56 LSF for **Ca **Ti 2v decay. Each bar corresponds 
to a contribution to the matrix element. Notice the interfering 
positive and negative contributions. 



In the p/-shell there is only one double beta emit- 
ter, namely "^^Ca. For many years, the experimental 
information was limited to a lower limit on the 2f-/3/3 
half-life T2^2> 3.6x1019 yr (Bardin et a/., 1970). The 
calculation of the ***Ca half-life was one of the first 
published results stemming from the full ^)/-shell cal- 
culations using the code ANTOINE. The resulting ma- 
trix elements are: M§^(0+)=0.083 and M^^(2+)=0.051. 



51 



For the ground state to ground state decay ^Ca^^Ti, 
G2^=l.lxl0-" yr-^ (Tsuboi et ai, 1984). The 
phase space factor hinders the transition to the 2+ 
that represents only about 3% of the total probabil- 
ity. Putting everything together, and using the Gamow- 
Tcllcr qiicnching factor, ahoady discussed, the resulting 
half-life is ^1'/^= 3.7x101^ yr, (Caurier et ai, 1990) 
(see also erratum in Caurier et at, 1994). The pre- 
diction was a success, because a later measure gave 
Tl'/2= 4.3li;i[stat]±1.4[syst] xlO^^ yr (Balysh et al, 
1996). 

Among the other (Sfi emitters in nature (around 30), 
only a few are potentially interesting for experiment since 
they have a Q/3/3 value, sufficiently large (> 2.5MeV) for 
the Oiy signal not to be drowned in the surrounding natu- 
ral radioactivity. With the exception of ^^"Nd, all of them 
can be described within a shell model approach. The re- 
sults for the 2iy mode of the lightest emitters '***Ca, ^^Ge 
and ®^Se, for which full space calculations are doable, 
are gathered in table XVII (Caurier et al., 1996). The 
SMMC method has also been applied to the calculation 
of 21^ double beta decays in (Radha et ai, 1996). 



TABLE XVII Calculated Tf/^ half-lives for several nuclei and 
0"*" — »■ 0"*" transitions 



Parent 


4«Ca 




825e 


T?/2 th.(y) 


3.7 X 10^" 


2.6 X 10^^ 


3.7 X 10^" 


^1% exp.(y) 


4.3 X 10^^ 


1.8 X 10^^ 


8.0 X 10^^ 



The expression for the neutrinoless beta decay half-life, 
in the 0+ — > 0+ case, can be brought to the following form 
(Doi et al, 1985; Takasugi, 1981): 

[t(;,^)(o+->o+]-^ = 

where {m,y) is the effective neutrino mass, Gqi, the kine- 
matic space factor and M^^^ and Mp'^^ the following 
matrix elements (m and n sum over nucleons): 

M^GT = (O/ II Yl Hr){^n ■ &m)tntm. ||0 + ) (96) 

where, due to the presence of the neutrino propagator, 
the "neutrino potential" h{r) is introduced. In this case, 

the matrix elements arc just expectation values of two 
body operators, without a sum over intermediate states. 



This was believed to make the results less dependent of 
the nuclear model employed to obtain the wave func- 
tion. An assumption that has not survived to the actual 
calculations. Full details of the calculations as well as 
predictions for other Ov and decays in heavier double 
beta emitters can be found in Retamosa et al. (1995) and 
Caurier et al. (1996) (see also Suhonen and Civitarese, 
1998, for a recent and very comprehensive review of the 
nuclear aspects of the double beta decay). The upper 
bounds on the neutrino mass resulting of our SM calcu- 
lations, assuming a reference half-life Tj^2 ^ 10^^ V- 
collected in table XVIII. 



TABLE XVlll Ov matrix elements and upper bounds on the 
neutrino mass for Tj/j ^ 10^^ y. (mi/) in eV. 



Parent 










0.63 


1.58 


1.97 




-0.09 


0.19 


-0.22 




0.94 


1.33 


0.49 



For the heavier emitters, some truncation scheme has 
to be employed and the seniority truncation seem to be 
the best, since the dimensions are strongly reduced in 
particular for 0"*" states. An interesting feature in the 
calculation of the double beta decay matrix elements is 
shown in Fig. 57 for the ^^Gc case : the convergence of 
the Mq^ matrix element is displayed as a function of the 
truncation of the valence space, either with the seniority 
V or with the configurations t (t is the maximum number 
of particles in the g9 orbital). The t = 16 and v = 16 
values correspond to the full space calculation and are 
consequently equal. The two truncation schemes show 
very distinct patterns, with the seniority truncation be- 
ing more efficient. Such patterns have been also observed 
in the Tellurium and Xenon isotopes. This seems to in- 
dicate that the most favorable nuclei for the theoretical 
calculation of the Oiy mode would be the spherical emit- 
ters, where seniority is an efficient truncation scheme. 



C. Radii isotope sliifts in the Calcium isotopes 

We have already seen in section VI. E that deformed 
np-nh configurations can appear at very low excitation 
energy around shell closures and even to become yrast 
in the case of neutron rich nuclei. This situation simply 
reflects the limitation of the spherical mean field descrip- 
tion of the nucleus, and shows that even in magic cases, 
the correlations produce a sizeable erosion of the Fermi 
surface. This is the case in the Calcium isotopic chain. 
Experiments based on optical isotope shifts or muonic 
atom data, reveal that the nuclear charge radii (r^) fol- 
low a characteristic parabolic shape with a pronounced 
odd-even staggering (Fricke et al., 1995; Palmer, 1984). 
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FIG. 57 Variation of the (m^) limit as a function of senior- 
ity truncations (circles) and configuration truncations (dia- 
monds) 



For the description of the nuclei around N=Z=20, a 
model space comprising the orbits Ida , 2s i , If? and 2p3 
is a judicious choice (sec Cauricr et at, 2001a). The inter- 
action is the same used to describe the neutron rich nuclei 
in the sd-pf valence space, called sdpf.sm in section VII, 
with modified single particle energies to reproduce the 
spectrum of ^^Si. 

Several important features of the nuclei at the sd/pf 
interface are reproduced by the calculation, among oth- 
ers, the excitation energies of the intruder O2 states in 

the Calcium isotopes as well as the location of the | 

states in the Scandium isotopes (sec Fig. 58), the excita- 
tion energies of the 2+ and 3^ and the B(E2)'s between 
the 2f and the states in the Calcium isotopes. 



functions, can be expressed as: 

1 



6ri{A) = -n;f{A)b^ 



ipf\^ju (98) 

where Z=20 for the calcium chain, b is the oscillator pa- 
rameter, and is extracted from the calculated wave 
functions. The charge radii isotope shifts of the Calcium 
isotopes relative to ''"Ca arc shown in Fig. 59, together 
with the experimental values. The global trends are very 
well reproduced, although the calculated shifts are a bit 
smaller than the experimental ones. This is probably due 
to the limitations in the valence space, that excludes the 
Ids , the If 5 and the 2p 1 orbits. 
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FIG. 59 Isotopes shifts in Calcium chain: experimental (cir- 
cles) versus shell model calculations (stars) 



4 - 



3-^ 



-1 



exp ^ 
/ 

/ 



+ G- 

3/2 inSc 



0* in Ca 

/ 

/ 

' k' 

sm 



J3- 



. . . • exp 



40 



42 



44 



46 



48 



FIG. 58 Comparison between calculated and experimental in- 
truder excitation energies in the Calcium and Scandium iso- 
topic chains 

Due to the cross shell pairing interaction, protons and 

neutrons are lifted from the sd to the /p orbitals. The 
former produce an increase of (r^) that, using HO wave 



D. Shell model calculations In heavier nuclei 

There are some regions of heavy nuclei in which phys- 
ically sound valence spaces can be designed that are at 
the same time tractable. A good example is the space 
comprising the neutron orbits between N=50 and N=82 
for the tin isotopes (Hjorth-Jensen et al., 1995; Nowacki, 
1996). When protons are allowed the dimensions grow 
rapidly and the calculations have been limited until now 
to nuclei with few particle or holes on the top of the 
closed shells (see Covello et al., 1997, for a review of the 
work of the Napoli group). 

The SMMC, that can overcome these limitations, has 
been applied in this region to ^^*Te and ^^*Xe that are 
candidates to 7 soft nuclei by Alhassid et al. (1996) 
and also to several Dysprosium isotopes A=152-162 in 
the Kumar-Baranger space by White et al. (2000). The 
QMCD method has been also applied to the study of 
the spherical to deformed transition in the even Barium 
isotopes with A=138-150 by Shimizu et al. (2001). 

Hints on the location of the hypothetical islands of su- 
per heavy elements arc usually inspired in mean-field cal- 
culations of the single particle structure. The predictions 
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for nuclei far from stability can be tested also in nuclei 
much closer to stability, where shell model calculations 
arc now feasible. In particular, recent systematic mean 
field calculations suggest a substantial gap for Z = 92 
and N = 126 (^^*U) corresponding to the Ihg shell clo- 
sure. On the other hand, the single quasiparticle energies 
extrapolated from lighter N=126 isotones up to ^^^Ac, do 
not support this conclusion. To shed light on these con- 
troversial predictions, shell model calculations (as well as 
experimental spectroscopic studies of ^^^Th (Hauschild 
et al, 2001)) were undertaken in the N = 126 isotones up 
to ^^*^U. Shell model calculations were performed in the 
lh|, 2f7, li^, 3p|, 2f5, 3pi proton valence space, using 
the realistic Kuo-Herling interaction (Kuo and Herling, 
1971) as modified by Brown and Warburton (Warburton 
and Brown, 1991). The calculation reproduces nicely the 
ground state energies, the 2+ energy systematics as well 
as the high-spin trends. The cases of ^^^Ra, ^^^Th and 
^^*U are shown in figure 60. The only deviations between 
theory and experiment are in the 3~ energy and reflect 
the particular nature of theses states which are known 
to be very collective and corresponding to particle hole 
excitations of the ^°®Pb core. 

No shell gap for Z=92, corresponding to the Ihg clo- 
sure, is predicted. On the contrary, the ground state 
of the N=126 isotones is characteristic of a superfluid 
regime with seniority zero components representing more 
than 95% of the wave functions for all nuclei. 




FIG. 60 Experimental vs theoretical spectrum of ^^''Ra, 



E. Random Hamiltonians 

The study of random Hamiltonians is a vast interdis- 
ciplinary subject that falls outside the scope of this re- 
view (see Porter, 1965, for a collection of the pioneering 
papers). Therefore here we shall only give a bibliograph- 
ical guide to work that has recently attracted wide at- 
tention and may have consequences in future shell model 
studies. Johnson et al. (1998, 1999) noticed that random 
interactions had a strong tendency to produce J = 
ground states. This is an empirical fact that was hith- 
erto attributed to the pairing force. Bijker and Frank 
(2000a, b); Bijker et al. (1999) showed that in an inter- 
acting boson context this also occurred, and was asso- 
ciated to the typical forms of collectivity found in the 
IBM. For the fcrmion problem no collectivity occurs for 
purely random interactions (Horoi et al., 2001). The ori- 
gin of J = ground state-dominance was attributed to 
"geometric chaoticity" (Mulhall et al., 2000). The ge- 
ometric aspects of the phenomenon were investigated in 
some simple cases by Zhao and Arima (2001); Zhao et al. 
(2001) and Chau Huu-Tai et al. (2003). Velazquez and 
Zuker (2002) argued that the general cause of J = 
ground state-dominance was to be found in time-reversal 
invariance, and showed that when the random matrix ele- 
ments were displaced to have a negative centroid, well de- 
veloped rotational motion appeared in the valence spaces 
where the realistic interactions would also produce it (as 



in Fig. 39). This is very much in line with what was found 
by Cortes et al. (1982), and the interesting point is that 
the collective ingredient induced by the displacement of 
the matrix elements is the quadrupolc force. 

It is unlikely that we shall learn much more from purely 
random Hamiltonians. However, the interplay of random 
and collective interactions may deserve further study. In 
particular: We know that a monopole plus pairing Hamil- 
tonian is a good approximation. Would it be a good idea 
to replace the rest of the interaction by a random one, 
instead of neglecting it? 

IX. CONCLUSION 

Nuclei are idiosyncratic, especially the lighter ones, ac- 
cessible to shell model treatment. Energy scales that are 
well separated in other systems (as vibrational and ro- 
tational states in molecules) do overlap here, leading to 
strong interplay of collective and single particle modes. 
Nonetheless, secular behavior — both in the masses and 
the spectra — eventually emerges and the pf shell is the 
boundary region were rapid variation from micleus to nu- 
cleus is replaced by smoother trends. As a consequence, 
larger calculations become associated with more trans- 
parent physics, and give hints on how to extend the shell 
model philosophy into heavier regions where exact diag- 
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onalizations become prohibitive. In this review we have 
not hesitated to advanee some ideas on how this could 
be achieved, by suggesting some final solution to the 
monopole problem, and exploiting the formal properties 
of the Lanczos construction. 

The shell model has been craft and science: one in- 
vented model spaces and interactions and forced them 
on the spectra. Sometimes it worked very well. Then 
one wondered why such a phenomenology succeeded, to 
discover that there was not so much phenomenology after 
all. It is to be hoped that this state of affairs will persist. 
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APPENDIX A: Basic definitions and results 

• p is the principal oscillator quantum number 

• Dp = {p + l){p + 2) is the (single fluid) degeneracy 
of shell p 

• Orbits are called r, s, etc. Dr = 2jV -|- 1 

• 6rs = ^jr-j, but r ^ s, and both have the same 
parity. 

• nir is the number of particles in orbit r, Tr is used 

for both the isospin and the isospin operator. In 
neutron-proton (np) scheme, nirx specifies the fluid 
X. Alternatively we use rir and Zr- 

• ^rstu ^rstu ^wo body matrix elements. W^(„ 
is used after the monopole part has been sub- 
tracted. 

A few equations have to exhibit explicitly angular mo- 
mentum (J), and isospin (T) conservation. We use Bruce 
French's notations (French, 1966): F stands for JT. 
Then (-1)'^ = (-1)-^+^, [F] = (2J + 1)(2T + 1), and 
in general F{T) = F{J)F{T). Also(-l)'' = (-1)>+V2, 
[r] = 2{2 jr + !)• Expressions carry to neutron-proton 
formalism simply by dropping the isospin factor. 



The one particle creation and anhilation operators 
can be coupled to quadratics in A and B, 



Xl^ (rs) = {ArAs)l^, Xrr.(rs) = (B.B.)F^, 



(A2) 



Zpp^(rs) = (1 -|-(5rs)~^^^^rr^(^*) the normalized pair 
operator and Ztt^ {ts) its Hermitian conjugate. 

For reduced matrix elements we use Racah's definition 



<aa,\Pl\pl3, >= 

y-Oiz Iz Pz 



(A3) 



The normal and multipole representations of H are 
obtained through the basic recoupling 



{Xl{rs)Xr{tu)f = -{-lT+'-^ - 



1/2 



^[r7]V2(-l)«+*-7-r Y sT\ ^^^^ 



whose inverse is 

(57,S7jo = (-1)"-*+-^ 

^[F7]^/^(-l)^+*---^ 




(A5) 



(4(rs)Xr(to))°. 



r * 7 

The "normal" representation of V is then 

r<st<u,r 

- E irsitunl^Vl,,^{xl^X^^^)\ (A6) 

(rstu)r 



where we have used 



(l + (5^,)-i/2 ifr<s, 
(1 -I- (5rs)^^^/2 if no restriction. 



(A7) 



so as to have complete fiexibility in the sums. Accord- 
ing to (A4), V can be transformed into the "multipole" 
representation 

(rstu)7 

+<5.t^™W'/'a;0„,,50j , (A8) 
where (sum only over Pauli-allowed F) 



55 



It follows that the form of HmT must be 



CO, 



u t ^ 



+,_,_r r s r I ^^^^^ 



Eq.(A5) suggests an alternative to (A8) 



(All) 



(57,si)°-(-ir+'- 



1/2 



where each term is associated with a pure two body 
operator. 



APPENDIX B: Full form of Hm 

If we were interested only in extracting the diagonal 
parts of the monopole Hamiltonian in jt scheme, HmT, 
the solution would consist in calculating traces of 7i, and 
showing that they can be written solely in terms of num- 
ber and isospin operators. Before we describe the tech- 
nique for dealing with the non diagonal parts of Ti.mT 
i. e. for generalizing to them the notion of centroid, some 
remarks on the interest of such an operation may be in 
order. 

The full Hm contains all that is required for Hartree 
Fock (HF) variation, but it goes beyond. Minimizing the 
energy with respect to a determinantal state will invari- 
ably lead to an isospin violation because neutron and 
proton radii tend to equalize (Duflo and Zuker, 2002), 
which demands different orbits of neutrons and protons. 
Therefore, to assess accurately the amount of isospin vio- 
lation in the presence of isospin-breaking forces we must 
cnsiire its conservation in their absence. More generally, 
a fidl diagonalization of Hm is of great intrinsic interest. 
Now, the technical details. 

Define the generalized number and isospin operators 

Srs = Srs[r]'^^S',° , Trs = ^^r«M'/'5°i (Bl) 

which, for Srs = 1, reduce to Srr = n-r, Trr — Tr- By 

definition Hm contains the two body quadratic forms in 
Sr<i and o i 



Jrtsu 



■ CrsCtui^rt^su ~ ^st^ru), (B2a) 



HmT — IC + ^^^{ttrtsu^rtsu ~ brtsuTrtsu)^rt^suj (BS) 



all 



where the sum is over all possible contributions. This is a 
special Hamiltonian containing only A = terms. Trans- 
forming to the normal representation through Eq. (AID) 

a;Oi = <5;,o^,i^V/l = 3M-V2,v/;,„ = -[rs]-V2 
and therefore 

Srtsu = ^ -^rsr ■ ^tuF = (B4a) 

r 



-Trtsu - 4 X] ^rsJO ' ^tuM ~ ^ X] -^IsJl ' ^tuJl (B4b) 

and inverting 



rtsu) 



(B5a) 



^IsJl ' ^tuJ'i- — -^{'^Srtsu + 4T'rts„). (B5b) 

J 

When jr ~ is = jt = ju and r ^ s and t ^ u, both 
Srust, Trust and Srtsu, Trtsu are present. They can be 
calculated from (B5a) and (B5b) by exchanging t and u, 

= -^{Srust — ^Trust) (B6a) 
= ^{^Srust + 4:Trust) (B6b) 

and by combining (B5a), (B5b), (B6a) and (B6b) 

2^ ^rsjo ■ ^tujo ^ (B7a) 



Trtsu = CrsCtu{TrfTsu-^SstSru), (B2b) Zlj, ■ ZtuJl^^^^^^ = (B7b) 



2 

J 



which in turn become for 5rt = f^su = 1) n^rs and Trs in 1 

Eqs. (10a) and (10b) ^((35'rtsu + 4:Trtsu) ± (JiSmst + -^Trust)) 
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To write HmT we introduce the notations 



then 



$(P) = 1 - (1 - 5rs){l - 5tu), (B8) 
^{e) = {5rs-5rs){6tu-5tu), (B9) 

So for r = s OT t = u, '^{P) = 1 and for = js = jt = ju 
and $(P) = 0, $(e) = 1. Then 



-,PT ^pT 
rstu 



r<s 
t<u 
T,p=± 



p = Sign (-1)^ f2^,j„ = ^ zIj^ ■ ZtujT 



a 



±T 
rstu 



JT ' ^tuJT 



(1± (-!)■') 



(BIO) 



W™t = /C + 5^(1 - $(e))(a„5„ + 



+<^(e)(a'^Sn 


+ b'^T„ 




' a. Ct) 


, with 






= 


4 


+ v°) 






hoc 


= ^(V^ 
4 


-v°) 


ao = 




+ 3V-1 




- V-") 


«a = 




- 3V-^ 




-v-°) 


b'^ = 
"a 


i(3V^'+3V-' 








5(3V« 


- 


+ v„+" - 


- V-") 



The values of the generahzed centroids V^stu and VpT are 
determined by demanding that H — HmT = contain 
no contributions with A = 0. In other words 



— 5rt5s 



' ' rstu ^ rstu 

(1 - $(e))vL„ + HeWrlu 



(Bll) 



must be such that cOrtsu — Oi from Eq. (AlO) 

(J) (J) (J) 

Applying this prescription to all the terms leads to (ob- 
viously SrtSsu = 1 in all cases) 

vL. = E^"^«M/Ew 
(J) (J) 

1 Dr{Ds + 2^P){-ir) 
4 1 + $(P) 

vf'L = Y.^^Lm ± (-!)■')/ E["^](i ± (-1)') 

Y^[J]{l±{-iy) = lD,{D^T2) 

(B13) 



Ei^] = z 



J 

D, 



$(e) = 1, for vfl . 



Through cqs.(B5a), (B5b), (B7a), and (B7b) we can 
obtain the form of TimT in terms of the monopole opera- 
tors by regrouping the coefficients affecting each of them. 
To simplify the presentation we adopt the following con- 
vention 

\ a = rstu r < s,t < w, SriSgu = 1, BUT 

I Sec — Srtsu^ Sqc = Srust-, Pa — Prtsu^ Pot — Prust^ 



1. Separation of Hmnp and Timo 

In j formalism Timnp is TimT under another guise: neu- 
tron and proton shells are differentiated and the oper- 
ators Prs and Srs are written in terms of four scalars 
Sr^sy-, x,y = n or p. 

We may also be interested in extracting only the purely 
isoscalar contribution to HmT, which we call Hmo- The 
power of French's product notation becomes particularly 
evident here, because the form of both terms is identical. 
It demands some algebraic manipulation to find 



Hn 



. or Hmo = ^ + E 



+ 2(C^^ + '^a)Sc + (V- - V^)5a)$(e) 



(B15) 



with 



v.«t„ = E^-t«[r]/E[r] 

(r) (r) 

J^[r] = Dr{D,-^{P))/{l + ^P)) 

(r) 

r 

E[r](l ± (-1)^+^'-) = Dr{Dr T i-lfl (B16) 



Of course we must remember that for H 



mnp-i 



Dr = 2jr + l, {-ly = -l,r = J, etc., while for H 



mO 



Dr = 2(2> + 1), (-1)^'' = -Fl, r = JT; etc. 

It should be noted that Hma is not obtained by simply 
discarding the b coefficients in eqs.(B14), because we can 
extract some 7 =00 contribution from the operators. 
The point will become quite clear when considering the 
diagonal contributions. 
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2. Diagonal forms of Hm 

We are going to specialize; to tiie diagonal terms of 
HmT and Hmnp^ which involve only V^^^^ matrix ele- 
ments whose centroids will be called simply Vrs and 

T 

(the overline in V^g^y^ was meant to avoid confusion with 

possible matrix elements V^^f^ or V^gj„, it can be safely 
dropped now). Then (B14) becomes Eq. (11) and (B15) 
becomes 



Knp or = }C + Y,Vrs nr{ns-Srs)/il + Srs) (B17) 



A coordinates and momenta while only A—loi them can 
be linearly independent since the solutions cannot de- 
pend on the center of mass coordinate R = Q2i ri)/VA 
or momentum P = (^iPi)/\^- The way out is to im- 
pose a factorization of the wavefunctions into relative and 
CM parts: 3>(ri r2 ... ta) = $rei 4>cm- The potential en- 
ergy is naturally given in terms of relative values, and for 
the kinetic energy we should do the same by referring to 
the CM momentum, 



Pi)', (Cl) 



We rewrite the relevant centroids incorporating explic- 
itly the Pauli restrictions 



Vr 



Dr{Ds — Srs) 



(B18a) 
(B18b) 



and change accordingly /Cy in Eq. (2). As we are only 
interested in wavefunctions in which the center of mass 
is at rest (or in its lowest possible state), we can add a 
CM operator to the Hamiltonian n^n + X{R^ + P^), 
and calling = r, — Vj we have 



i?2 



-Pi,), (C2) 



brs = - V°, (B18c) 

The relationship between and Vrs makes it possible 
to combine eqs. (11) and (B17) in a single form 



^1 = ^ + 5^ 7^ 



(1 + 5rs) 

4(£>.-l)' 



Vrs nr{ns - 5rs) + 



(B19) 



in which now the brs term can be dropped to obtain 

In the np scheme each orbit r goes into two r„ and Vp 
and the centroids can be obtained through {x,y = n or 
P, xy^y) 



Vr 



_ 1 

~ 2 
= V^, 



Vl 1 



2Sr. 
Dr 



VI 1 



25r_ 
ITr 



(B20) 



Note that the diagonal terms depend on the represen- 
tation: H!^np 7^ ^mT ill general. 



APPENDIX C: The center of mass problem 

What do you do about center of mass? is probably the 
standard question most SM practitioners prefer to ignore 
or dismiss. Even if we may be tempted to do so, there is 
no excuse for ignoring what the problem is, and here we 
shall explain it in sufficient detail, so as to dispel some 
common misconceptions. 

The center of mass (CM) problem arises because in a 
many body treatment it is most convenient to work with 




so that upon diagonalization the eigenvalues will be of 

the form E = Erei + M^CM + 3/2), and it only re- 
mains to select the states with Ncm=0- We are tak- 
ing for granted separation of CM and relative coordi- 
nates. Unfortunately, this happens only for spaces that 
arc closed under the CM operator (C2): for Ncm to be a 
good quantum number the space must include all possi- 
ble states with Ncm oscillator quanta. The problem was 
raised by Elliott and Skyrme (1955) who initiated the 
study of "particle-hole" (ph) excitations on closed shells. 
They noted that acting with R — iP on the IPM ground 
state of ^^O (|0)) leads to 



V2p3Si - VEpid3+p3d3 + md5)\0), (C3) 



where P2j removes a particle and si or d2j add a particle 
on |0). As R-iP has tensorial rank J'^T = I'D, Eq. (C3) 
is telling us that out of five possible 1~0 excitations, one 
is "spurious" and has to be discarded. 

Assume now that we are interested in 2p2p excita- 
tions. They involve jumps of two oscillator quanta (2/k<; ) 
and the CM eigenstates {R — iP)^|0) involve the op- 
erator in Eq. (C3), but also jumps to other shells, of 
the type spp{sd) = s (sd) or (sd) {pf)p (sd) = p{pf). 
Therefore, as anticipated, relative-CM factorization can 
be achieved only by including all states involving a given 
number of oscillator quanta. The clean way to pro- 
ceed is through complete iVfiw spaces, discussed in Sec- 
tions II. A. 1 and III.F. The 0/kjand EI spaces are also 
free of problems (the latter because no Ihio 1~0 states 
exist). It remains to analyze the EEI valence spaces, 
where CM spuriousness is always present but strongly 
suppressed because the main contributors to R — iP -of 
the type pj =^ p + Ij ± 1, with the largest j- are al- 
ways excluded. Consider EEI(l)=pi/2, ^5/2; Si/2- The 
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only possible lhu!l~0 state is PiSi, which according to 
Eq. (C3) accounts for (1/18)% = 5.6% of the spurious 
state. This apparently minor problem was unduly trans- 
formed into a serious one through a proposal by Gloeck- 
ner and Lawson (1974). It amounts to project CM spuri- 
ousness through the H TC+X{R^+P^) prescription in 
the EEI(l) space by identifying R — iP = pisi. The pro- 
cedure is manifestly incorrect, as was repeatedly pointed 
out (sec for instance Whitehead et ai, 1977) but the mis- 
conception persists. An interesting and viable alternative 
was put forward by Dean et al. (1999) in calculations with 
two contiguous major shells. Further work on the subject 
would be welcome. 
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